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792 HADUPTAUFSÄTZE I. 
Ein Pohlhausenverfahren zur Berechnung laminarer 
© kompressibler Grenzschichten an einer geheizten Wand 
a Ei | Von J. Ginzel in Göttingen PR akt, SBIHL UN, 
Das für inkompressible Strömung bekannte Pohlhausenverfahren wird durch Hinzunahme einer weileren m on 
gewöhnlichen Differentialgleichung (Energie-Integralsatz) auf den kompressiblen Fall erweitert und en AR 
Lösungsverfahren für das gekoppelte System zweier linearer Differentialgleichungen angegeben. Für wi 
fest gewählte Wandtemperaturen werden Rechengrundlagen für die nunmehr sehr bequeme numerische Bi- 


spielrechnung bereitgestellt. Zum Schluß werden Wandschubspannung und Wärmeübergangszahl in einem _ ne e) 3 
Beispiel diskutiert. a ; i v N 


The known Pohlhausen method for ihe caleulation of the laminar boundary layer of an incompressible Bi 
flow is extended to compressible flow by a further ordinary differential equation (energy integral iheorem). N 
. A method of solution for the two differential equations one depending on the other is given. There are auxi- } 
liary curves for the numerical solution of an example in two fixed cases of wall temperature. Finally, skin 0... 
[rietion and heat transfer .coeffieient are discussed for an example. ® ; 


‚La. methode bien connue de Pohlhausen pour la calculation d’une couche limite laminaire d’un fluide 
= incompressible est &iendue au cas compressible en regardant une deuxieme £quation differentielle ordinaire { 
Be (theoreme integrale de l’energie). Ume methode de solution est donnee pour le systeme simultane de deux N 
 equations differentielles lineaires. Des courbes auxiliaires pour ld calculation des exemples sont preparees 1 
pour deuz temperatures de la surface. Enfin le frottement a la surface et le co&fficivent du transport de la: 

chaleur sont discul&s pour un ewemple. 

Meron HlonpraysenHa, IpuMeHseMmbIä IIPu MCCAIeNOBAHHN IIOTOK HECiKHMAMITeÄCH YKHUAKOCTH, 
PACHPOCTpaHAeTcaıyTeM BEIIWOYeHHA BTOPOTO UPOCTOrO AuhepeHnmanbHoro YpaBHeHua (Teo- 
PeME HUHTerpana 3Hepruu) Ha CITYAY CHKUMAIMINETOCH IIOTOKA. aerca MeTon pemeHun 9ToA 

CHCTEMBI ABYX JIMHeÜHBIX AuhepeHHHanbHsIx ypaspeHnä. I'nıAa NByX NOCTOAHHEIX TeMIECHATYP 
CTEHBL, HAITCA BCHOMOTATENBHLIE AUHATPAMMBI, SHAYUTENIBHO OÖAETYAWINHEe YHCHEHHOE HSCHE-. d 
MOBAHHE YACTHBIX CKyyaes. HaxoHen o6cy3kNamTcaA HA IPHMepe TPaHnyHLe HAaupmkeHuna Van Br 
. cABura (HA CTeHe) uU KO9bHNHeHT Hepexona Temna. | | S 


1. Bisheriger Stand der Arbeiten über kompressible laminare Grenzschichten Be 
Der heutige Stand der Forschung in Deutschland auf dem Gebiet der kompressiblen } “er 
Grenzschichten ist dargestellt in dem Artikel B7 von M. Schäfer der Göttinger Mono- ar. 
graphien [1]. Da die Monographien nicht überall sehr leicht zugänglich sein dürften, müssen 2 
wir ein kurzes Referat vorausschicken. | SR 
Zuerst ist die strenge Lösung des Differentialgleichungssystems (s. Gl. 1 bis 3 des nächsten 
Abschnittes) für den Fall der ebenen Platte zu nennen, die vnHantzscheundWendt 
angegeben worden ist [2], [3]. Die Lösungsmethode deckt sich teilweise mit einem vonGrocco 
[4] angegebenen Verfahren. Es handelt sich um die Bestimmung der Geschwindigkeitsver- 
teilung u (x, y) und der Temperaturverteilung 7 (z, y) in der Grenzschicht, wo x die Koordinate 
in Plattenrichtung, y die Koordinate senkrecht zur Platte ist. Für den Fall der Platte zeigt 


es sich nun, daß die Lösungen einander ähnlich sind, d.h. als Funktionen einer Variablen yly& 
geschrieben werden können. Das System der partiellen Differentialgleichungen geht dann in 
ein System gewöhnlicher Differentialgleichungen über. 

_— — -In zwei Spezialfällen werden die gewöhnlichen Differentialgleichungen voneinander un- 
abhängig und lassen sich nacheinander (z. T. sogar geschlossen) integrieren. Der erste Fall 
betrifft einen speziellen Zahlenwert der als Materialkonstante angesehenen Prandtlschen 
Zahlo = c, u[A (u = Zähigkeitskoeffizient, c, = spez. Wärme bei konstantem Druck, A = Wär- 
meleitfähigkeitskoeffizient). Für o = 1 werden die Gleichungen voneinander unabhängig, und 
das Integral der zweiten läßt sich sofort hinschreiben. (Die Prandtlsche Zahl für Luft 
beträgt o = 0,7). Eine ähnliche ausgezeichnete Rolle spielt noch ein anderer Zahlenwert. Wie 
wir im nächsten Abschnitt sehen werden, gehört zu dem vollständigen ‚Gleichungssystem außer 
den drei Differentialgleichungen und außer der Zustandsgleichung des Gases noch ein Gesetz 
über die Abhängigkeit des Zähigkeitskoeffizienten von der Temperatur. Diesem Gesetz kann 
man näherungsweise die Form eines Potenzgesetzes geben (s. Gl. (7) im nächsten Abschnitt). 
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u FA ar erfahren zur r 
"Der konstante Exponent der Potenz heiße ®, sein Wert muß dem jeweilige: 
‘peraturen ® = 0,7 bis © = 0,5. Fürro=1 zerfällt das System der gewöhnlichen 


"werte auf die Ergebnisse studiert. Außerdem liegt speziell für das sogenannte Thermometer- 
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angepaßt werden. Für tiefe Temperaturen ist der beste Näherungswert o= 1, für hohe Tem- 


gleichungen wieder in zwei voneinander unabhängig Gleichungen. Hantzsche undWend t 
haben nun Beispiele gerechnet einerseits für e = 1 und verschiedene Werte von ® und anderer- 
seits ® = 1 und verschiedene Werte von o und auf diese Weise den Einfluß der beiden Zahlen- 


problem (der nichtgeheizten bzw. ungekühlten Platte), das durch die Bedingung verschwinden- 
den Wärmeüberganges an der Wand charakterisiert ist, eine numerische Lösung des simultanen 
Systems gewöhnlicher Differentialgleichungen für o = 0,7 und eine Reihe von wa-Werten vor, 
die von Grocco [4] stammt. } | 

Es hat sich gezeigt, daß für das Thermometerproblem der Einfluß des Exponenten & 
auf die Lösung gering ist, die Geschwindigkeits- und Temperaturprofile werden durch den 
leicht zu integrierenden Fall ® = 1 gut wiedergegeben. Dagegen wird der Reibungswiderstand 
vom Exponenten & des Zähigkeitsgesetzes beeinflußt. DiePrandtlsche Zahl beeinflußt für 
höhere Machsche Zahlen die Geschwindigkeitsprofile beim Thermometerproblem, wogegen 
der Reibungswiderstand von der Prandtlschen Zahl nur unbedeutend abhängt. Die Ab- 
hängigkeit der Machschen Zahl von der Temperatur (für den Fall der geheizten und ge- 
kühlten Platte) kann man aus den Kurvenblättern in [2] und [3] entnehmen. 

Die Resultate vonHantzscheund Wendt werden gestützt und ergänzt durch die 
Ergebnisse einer neueren Arbeit von Crocco [5], die sich wieder mit der laminaren Grenz- 
schicht an der ebenen Platte befaßt. Neu an dieser Arbeit ist eine ausführliche Diskussion des 
Zähigkeitsgesetzes, und zwar hat Crocco eine Menge Fälle durchgerechnet erstens für ver- 
schiedene Exponenten des Potenzgesetzes und zweitens für verschiedene Parameterwerte der 
vollen Sutherlandschen Formel (s. Gl. (6) des nächsten Abschnitts). Dabei zeigt sich, 
daß die Enthalpieverteilung vom Zähigkeitsgesetz praktisch unabhängig ist, von der Heizung 
Gplda (ip Enthalpiewert an der Wand, ?, Enthalpiewert in der äußeren Strömung) und vom Qua- 
drat derMach schen Zahl linear abhängt. Die Unabhängigkeit der Temperaturverteilung vom 
Zähigkeitsgesetz wurde auchvonHantzscheundWendtifestgestellt (s. oben). Die lineare 
Abhängigkeit vom Quadrat der Mach -Zahl wird auch vnGopeundHartree[$]benutzt 
(s. weiter unten). Das für den Spezialfall@ = 1 gewonnene Integral der Energiegleichung (s. 
Gl. (3)) gilt also praktisch auch für Fälle #1. Für diese Fälle muß nur die Grenzschicht- 
gleichung (1) neu integriert werden. Der Reibungsbeiwert c, ist vom Zähigkeitsgesetz merklich 
abhängig, wenn der Temperaturbereich bzw. Machzahlbereich groß ist, wie Croccoin vielen 
Abbildungen zeigt. In den Abbildungen ist c, YRe aufgetragen, wo c, der Reibungsbeiwert und 
Re die mit der Plattenlänge gebildete Reynoldssche Zahl der Außenströmung ist. Die im 
Inkompressiblen auftretende Abhängigkeit von der Reynoldsschen Zahl ist in den Abbil- 
dungen also nicht enthalten, denn im Inkompressiblen ist ce, Yke eine Konstante. 

Das Problem der laminaren Grenzschicht an der Platte (verschwindender Druckgradient) 
kann also als gelöst gelten. Für das Thermometerproblem verschwindenden Wärmeüberganges 
an der Wand ist auch noch der Fall des konstanten Druckgradienten vollständig behandelt 
worden von Gopeund Hartre[8]. Die Verfasser benutzen für den Energiesatz (3) die für 
die Prandtlsche Zahl o=1 geltende Lösung als Näherungslösung (s. Gl.(11)). Mit dieser 
Näherungsannahme ist die weitere verknüpft, daß Geschwindigkeits- und Temperaturgrenz- 
schichtdicke gleich groß sind. Die Verfasser haben es dann nur noch mit der Grenzschicht- 
gleichung (1) zu tun und schließen ihr Verfahren an die vonHowarth[9]für den inkom- 
pressiblen Fall und für den Fall des konstanten Geschwindigkeitsgradienten entwickelte Me- 
thode der Reihendarstellung an. Die numerische Berechnung wurde mit Hilfe der Eniae durch- 
geführt, die die Anwendung eines Iterationsschemas erlaubte. Auf: diese Weise wurden die 
Koelfizientenfunktionen mit großer Genauigkeit für einen großen Mach zahlbereich 
(0 <M <A) berechnet. 

Die nächste Aufgabe besteht in der Lösung für den Fall einer beliebigen Kontur bzw. be- 
liebigen Geschwindigkeitsverteilung in der Außenströmung. Es ist bekannt, welche Schwierig- 
keiten die exakte Behandlung der Grenzschichtgleichung schon im inkompressiblen Fall bietet, 
Es scheint hoffnungslos, das Problem für den kompressiblen Fall angreifen zu wollen. Nun 
haben sich in praktischen Fällen im Inkompressiblen die sogenannten Mittelwertsverfahren oder 
Impulsverfahren, die an den Namen von Pohlhausen geknüpft sind, sehr bewährt. Man 
erfüllt die Grenzschichtgleichung nicht exakt, sondern nur im Mittel, Das erreicht man durch 
Integration der partiellen Differentialgleichung über die eine Variable (senkrecht zur Kontur). 
Die partielle Differentialgleichung wird damit zu einer gewöhnlichen, die man mit irgend- 
welchen bekannten Näherungsmethoden lösen kann. Es liegt nahe, dieses Verfahren auch im 


..  zialfall Prandtlsche Zahl o =1. In diesem Fall ist die Abhän gigkeit der Temperaturverteilu 

bzw. Zähigkeitsverteilung von. der Geschwindigkeit ( Lösung der Gl. (3)) bekannt (s. Gl. (11 

so daß das ganze Verfahren auf die Integration der gemittelten Grenzschichtgleichung hinaus- 
läuft. Die Rechnung verläuft wie im Inkompressiblen, nachdem als Hilfsmittel der IMogTatinaR 
‚ eine Kurvenschar mit der Machschen Zahl als Parameter bereitgestellt wurde. 


ebenen Platte und das Thermometerproblem: für den Fall des konstanten Druckgradienten Br 
en  Näherungslösung für Prandtlsche Zahl o = 1 mit Hilfe eines Impulsintegralverfahrens (Pohl- 
hausen - Verfahren) vor. Die vorliegende Arbeit stellt ein Pohlhausen verfahren Turaaı 


‚eine beliebige geheizte oder gekühlte Wand dar, wobei für die Prandtlsche Zahl und dene 
BE  üaeneen En Zähigkeitsgesetzes keine speziellen Zahlenwerte vorausgesetzt Nach s 


j Grenzschichtgleichung, Kontinuitätsgleichung ‚und Energiegleichung bezeichnet. Sie lauten 


Wert 113. Für kleine a 


N >r ers :h in dieser 
]. Os wa " its ch beschränkt sich auf ds Trdmomeis oben und Behandeit‘ nur den Sp 


Zusammenfassend stellen wir fest, daß das Problem der laminaren Grenzschicht a an der 


vollständig gelöst sind. Für das Thermometerproblem einer beliebigen Strömung liegt eine 


2. Die Ausgangsgleichungen 
Die Grundlage der Rechnung bilden drei partielle Daerentialeterunn en die man als 
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Dabei bedeuten v und v die Geschwindigkeitskomponenten in Wandrichtung «© und senkrecht 
zur Wand (y-Richtung), p den Druck, u den Zähigkeitskoeffizienten. © heißt Wärmeinhalt 


oder Enthalpie, und o ist die Prandtlsche Zahl 
_ pU | 2 

| EN RE (4), % 
wo c, die spezifische Wärme bei konstantem Druck und A der Koeffizient der Wärmeleit- B, 
ee ist. Wir setzen im folgenden konstante Prandtlsche Zahl o und konstante Wärme ba 
„ der Luft voraus. Dann ist für ideale Gase Be; 
RE ER NEAN (5), By 

wo T die absolute Temperatur ist, die nicht in Graden gemessen ist, sondern mit der Grad- Far “ 


messung T° durch die Formel 7 = T° J.(J mechanisches Wärmeäquivalent) zusammenhängt. 
Wir brauchen nun noch einen Zusammenhang zwischen dem Zähigkeitskoeffizienten und der 
Temperatur. Dieser Zusammenhang wird gut durch die Sutherlandsche Formel be- | 


schrieben \ 
5 wi N ES ARE y, 
wo T, die Ruhetemperatur oder Kesseltemperatur ist. Für Luft hat die Konstante ( den 


u ‚ also nicht zu große Spannen zwischen Temperatur und Ruhe- 


2 

temperatur kann man das Gesetz entwickeln. 
I RN ee ! 7 
a nr en, 


Für den Koeffizienten » kann man bei tiefen Temperaturen den Wert 1, bei hohen Tempe- 
“ raturen den Wert 1/2 einsetzen. 


Während in den bisher besprochenen Arbeiten (außer in [6]) das Problem der Lösung 
der Differentialgleichungen (1) bis (3) direkt in Angriff genommen wird, wollen wir ein ein- 
faches Mittelwertverfahren entwickeln, d.h. wir wollen die Gleichungen (1) und (3) unter Be- 
rücksichtigung der Gl. (4) bis (7) in y-Richtung nicht exakt, sondern nur im Mittel erfüllen. Wir 
integrieren für diesen Zweck die Gleichungen (1) und (3) über y unter Berücksichtigung der 
Randbedingungen an der Wand und in der Außenströmung. An der Wand gilt, wenn wir die 
geheizte oder gekühlte Kontur betrachten, 

y=V0: v=v=0 d= ip = konst 
21* 
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Br Bedingungen an der Wand wäre das v 3 ] Y en le? a: 2 WERL 4 nd a7 

ae in r In ara er ee w ın 1. 3 bedeutet ae en 
Be en h. die y-Entfe Benieiiperath 
ist; ö4 die Temperaturgrenzschichtdicke, An die g-Entlernäig; far e die A ® Ta Ri 

. ‚erreicht ist. Nur nr 
Wir betrachten Überschallströmungen und stellen jetzt einige gun auf, die 

wir immer y wieder brauchen. Für die Schallgeschwindigkeit c gilt bei adiabatischer Ström 
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SE | A An gi 
wo “ De instant | 
En R=J (9%) r alas 
ist. c, ist die spezifische Wärme bei konstantem Volumen, x ist das Verhältnis der spezifischen BL: 
Wärmen Se * 
R ä , “= Co[% 07 ie ee ae EL Ka re i un 9). \ > 
Dann gilt Zr | er‘ 
BE ”»—1 ; ’ | 
erg) 3 U: 
” ” ” ” 1 ” 
PORN Ra al | c® 1 *#%—G)To *%—1 up 
m m U2 2 
Bi ‚Nun ist nach den Energiesatz für die AUDERSEUDTUDEN die wir als adiabatisch Vor | 
ER‘ haben, 
Br: ..,w = Ö k { nr +%(7T—T)=0 
E Das 0er 2 re Eh PETE ET. 11) 
F i) 3 4 Gg \ )- 
N = 7, konst 
N fr | 
l Setzt man ‚diese Beziehung in (10) ein, so erhält man nach einiger Rechnung 
B: Me—)- ET: In De 2 ante u RC (12) 
k Buy unter Benutzung dieser Beziehung ae aus (10) 
| ia ul —1) - ie ( —1 ) 
| ah RT 1 + 1 RE (13). 


Die hiermit eingeführte Zahl A, hängt von der Außenströmung ab. Aus der Definition der 
Schallgeschwindigkeit 


» __ d0a 
ce dp» 
und aus der Bernoullischen Gleichung leiten wir ferner her 
U 0a 
Mia Te A RE 14), 
SF (14) 


.wo die Striche Ableitungen in Strömungsrichtung bedeuten. 

Wir kehren zurück zu den Gleichungen (1) bis (3) und wiederholen die Herleitung des 
Impulsintegralsatzes, den wir in der vonHolsteinundBohlen[7J, Oswatitsch [6] 
und W.alz [11] angegebenen Form benutzen wollen. Zunächst machen wir die in der Grenz- 
schichttheorie übliche Voraussetzung, daß der Druck quer zur Wand (in y-Richtung) konstant 


ist und erhalten aus der Bernoullischen Gleichung für den Druckanstieg in PLÖN: 
richtung 


dp 2 ar der 

Mn am al 0} U U U = EN . 
Dieser Wert soll unabhängig von y in der ganzen Grenzschicht gelten. Die Zustandsgleichung 
des Gases heißt 


DaßRol....ı, in 7. ....(18) 
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Der erste Wert heißt Verdrängungsdicke, weil er ein Maß für die Strecke ist, um die die äuße 


Strömung durch die Grenzschicht vom umströmten Körper abgedrängt wird. Die zweite Größe 


 mißt den Impulsverlust, den die Strömung unter der Wirkung der Grenzschicht erieideh. a 
een daher SHIT SUSUUCHS- ‚Mit diesen Größen heißt die letzte Gleichung _ % 
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oder nach leichter Uimformung und Muliplikation mit 


“N 


Das ist eine lineare Differentialgleichung, auf deren Integration wir eingehen, sobald wir aus 
der Ausgangsgleichung (3) die entsprechende Gleichung abgeleitet haben. 
Wir addieren zu Gleichung (3) die mit @ multiplizierte Kontinuitätsgleichung (2) 


doui en ; ch a en 
0% + 70 an Aere 


und erhalten unter nochmaliger en der Kontinuitätsgleichung und Integration über 4 


"unter Berücksiehtigung ‚der Randbedingungen 
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Diese Gleichung ist nicht neu, sie findet sich in etwas anderer Schreibung bei Goldstein 
[12] (Bd. 2 Gl. (39) S. 615). Unter Benutzung der Gleichung (11) kann man schreiben 


d 1 ee ) 
4 m A ee dy .. (21). 
Efeu -iway+ vu [ua—ody -—e), +[ul%)as . eu 
Das letzte Glied heißt Dissipationsfunktion. Wir definieren jetzt eine Enthalpieverlustdicke 
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ir ersetzen i„ durch den aus Gl. (12) folgenden sr dividieren a Gihung durch 0 
und aa nach A. des ersten Gliedes 
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‘Während wir für Satz eo) = Wandabstand mit der sen Gimenatidalog ab, 
macht haben, werden wir ihn beim Temperaturprofil mit der Enthalpieverlustdicke dueBTE SE 
Sale machen und erhalten nach Umformung des ersten Gliedes ’ 
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Die Querstriche unter dem Logarithmus sollen andeuten, daß man die dimensionsbehafteten 
Größen durch eine für die betreffende Strömung feste Größe derselben Dimension (etwa An- 
strömgeschwindigkeit und Ruhedichte) dividieren soll, damit die Gleichung dimensionsrichtig 
bleibt. In der Form (26) ist der Energieintegralsatz für die Integration besonders geeignet, 
denn die Gl. 129) ist wieder eine N Differentialgleichung. Sie hat die Gestalt 


TE +epl@) +rla)=0, 


wobei in unserem: Falle 


ua 
p(x) =z 7, Me 
Fe RE 0. (27) 
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Größen pun AFrastachtlens im wesentlichen vo ht: von n2: abhän 

a all rdings erst die numerische Rechnung gezeigt. Es wurd n einige Beispielrechnungen & 

. zwei Weisen durchgeführt. Das eine Mal wurde der Energieintegralsatz i in der Form (28) be- 
nutzt, das andere Mal in einer anderen umständlicheren (hier nicht angegebenen) Form, die 


nicht die Gestalt einer linearen Differentialgleichung hat. Die numerischen Resultate stimmten 
im Rahmen der Rechengenauigkeit überein. Wendet man auf (27) die bekannte Integrations- 
formel für lineare Er a an, so sieht man, daß sich eine Integration ausführen 
‚läßt, "und man erhält, wenn man Zur Baimenls der. rn einen Bea A 
wert 0 von re lenmtet R ; 


Betrachtung zeigt, daß (20) außer von ö (bzw. von z) noch von Ile und ebenso (28) AUBEEN 
von 6, (bzw. von z,) noch von d,/ö abhängt. 
fihe wir nun im nächsten Kapitel zur praktischen Integration der Gl. (20) ee (28) kom- 
men, muß einiges über den physikalischen Sinn des Energieintegralsatzes gesagt werden. 
Zunächst ist festzustellen, daß sich die Aussage unserer Gl. (24) bzw. (26) nicht mit dem 
Energiesatz deckt, den Wieghardt [13] zur Berechnung laminarer Grenzschichten ein- 


‚geführt hat und den auch Walz[1 1] benutzt. Wenn man von der Kompressibilität absieht, _ 


‚liefert j ja die Gl. (3) keine Aussage. Wann man die Grenzschichtgleichung (1) mit der Geschwin- 
digkeit « multipliziert, haben ihre Glieder die Dimension eines Energietransportes, durch ‚Inte- 
gration in y-Richtung bekommt man eine Aussage über den Zusammenhang der kinetischen 
Energie der Strömung mit der Dissipation, und diesen Zusammenhang hat Wieghardt 
für den inkompressiblen Fall aufgestellt und bei seinem Verfahren berücksichtigt. Wir wollen 
die entsprechende Beziehung für den kompressiblen Fall kurz herleiten. 

Wir addieren die mit w multiplizierte Grenzschichtgleichung zu der mit u?/2 multipli- 
zierten Kontinuitätsgleichung. Unter neuerlicher Benutzung der Kontinuitätsgleichung ergibt 
sich dann 

2 
‚oou ? 


eier ae : 0 ou a u\? 
= dy —guUU= 2 (un ie — u ie) : 
0%. ey 2) 0% oy °y °Y 


Die Integration dieser Gleichung über y bis zur Außenströmung ergibt unter Beachtung der 


Randbedingungen nach Umformung des zweiten Gliedes 
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Im inkompressiblen Fall fällt das zweite Glied weg, und man kann das-Ergebnis dahin deuten, 
daß der Verlust an kinetischer Energie sich in Reibungsenergie umgesetzt hat. Im kom- 
pressiblen Fall tritt das zweite Glied, das von der Veränderlichkeit der Dichte abhängt, dazu. 
Gehen wir zurück bis zur Gl. (21), so sehen wir, daß dieses Glied auch dort auftritt. Nun muß 
bei Berücksichtigung der Kompressibilität in Gl. (29) auch die Gl. (21) erfüllt sein. Wir können 
also den Energiesatz auch so formulieren, daß wir die beiden Gleichungen addieren und erhalten 


d DW. ) 3 = ( = 
an ea) 2 ee ren an) ERS 


Hier sieht man, daß die kinetische Energie und die Wärmeenergie in der Strömung mit dem 
Wärmeübergang von der Wand in die Strömung im Gleichgewicht sein müssen. Die Einfachheit 
der Beziehung (30) legt den Gedanken nahe, sie der weiteren Rechnung zugrunde zu legen. 
Das bewährt sich aber nicht, weil wegen (12) die beiden Bestandteile des ersten Gliedes von 
derselben Größenordnung sind. Mindestens für die geheizte Kontur, für die der Wärmeinhalt 
in der Grenzschicht größer ist als in der Außenströmung, steht dann auf der linken Seite eine 
kleine Differenz zweier fast gleicher Größen. Das ist aber rechentechnisch für ein Näherungs- 
verfahren sehr ungünstig. Die Gl..(21) bzw. (28) weist diese rechentechnischen Schwierig- 
keiten nicht auf und soll deshalb dem nun zu beschreibenden numerischen Verfahren zugrunde 
gelegt werden. 


Die Integration des Bereits ist also auf eine Quadratur zürickgefthrt. Es = vs 5 
aber zu bemerken, daß die Gl. (20) und (28) gekoppelt sind, denn in r kommt außer 9, auch 9, 
vor; ©, und 9, hängen aber über o nach (16) von 4 und damit von 6, oder 9, ab. Dienähere 


RD N REN wie Bin 
' daß das Pohlhausen verfahren im Inkomp: erheblich 
- außer dem Impulsintegralsatz (20) noch den Energieintegralsatz (29) (zweites ( 
kompressiblen) erfüllt. Man erzielt auch noch eine aeg dem alten Pohlhau r 
ren erhöhte Genauigkeit, wenn man auf die Erfüllung der iet 
‘statt dessen die Bedingung (29) befriedigt, weil diese Aue Mittelwert über die ee AT 
schicht darstellt. Nach diesen Erfahrungen wäre es also auch in unserem Falle das 
außer den Gl. (20) und (28) noch die Gl. (29) zu erfüllen. Wir haben aber die 
Schwierigkeiten bei der Erfüllung dreier simultaner Differentialgleichungen gescheut und 
zu den Gl. (20) und (28) wieder die alten Grenzschichtbindungen dazu genommen. Das Er Ä 
fahren kann also keine größere Genauigkeit als das alte Pohlhausen verfahren (im 
nach Holstein und Bohlen) beanspruchen. 3 


3. Bestimmung der gewählten Profilklasse aus den Randbedingungen 


Außer der schon aus al. (13) definierten Größe A, benutzen wir den Pohlhaus en- a KT | 
Parameter A, rg 


Bi 2 1-8. NER, Se Ur Be 0 a 


u 2 


Dr 4 Neben die Geschwindigkeitsgrenzschichtdicke ö tritt jetzt die Tenperahrehenai ö, 

: als die Entfernung von der Wand, in der die Temperatur und der Wärmeinhalt ihren der Außen- 

strömung entsprechenden Wert 7, bzw. i, angenommen haben. 

Bi} Wir machen die Geschwindigkeit und die Enthalpie mit ihren Werten der Außenströmung 

Fe dimensionslos 


DR) dag 
Den Wandabstand machen wir für das Geschwindigkeitsprofil mit 6, für das Temperatur- 
profil mit ö, dimensionslos 


und erinnern an die Gl. (7) und (16). 
Als Profilklasse wählen wir die von Mangler [14] angegebene Form der Pohl- 
“hausen profile 


i | a | BTRON, 9 (33). 
t=1— (1)? (b, +5, nı) 


Durch diesen Ansatz sind die Randbedingungen am äußeren Ende der Grenzschicht schon 
erfüllt, denn dort ist 


07 a? n? Ri Mm 7 
Am inneren Rande muß % = 0 sein, was auch erfüllt ist. Die Konstante b, ergibt sich aus der 
Forderung %n=0)— —e = /, u 
& 
ea 1 N ERTRRe (34). 


Die noch freibleibenden Konstanten «a, und b, ergeben sich aus den Grenzschichtbindungen, 
d.h. aus den Gl. (1) und (3) für n = 0 bzw. 7, = 0. Diese Grenzschichtbindungen heißen 
R Pu du ou 

or ayoy 


0, due “ Tou\® 
a ayay en er 
Führt man in diese Gleichungen die Beziehungen (7), (1 6) und (32) ein, so hat man zunächst 
aus der ersten 


— 9, UU' 
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Im Inkompressiblen heißt diese Beziehung einfach 3 = — 4. In (36) sind nun für die Funk- 


tionswerte ihre Werte an der Wand einzusetzen. Das ergibt nach (33) 


61,430 281345) Ad | 
I RR N A URN (37). 


61402 (31,34 b,) 
1 


7 Das ist der Zusammenhang zwischen a, und b.. 

2 Behandeln wir die zweite Gl. (35) in derselben Weise, so erhalten wir unter Beachtung 
5 von (12) | 

R- A LE, +) LaB 3, re N Behr. en (38) 
5 mit 


Ö 2 f 
oma —1)2, BR RE EHEN, (39). 
Unter Beachtung der Gl. (37) kann man jetzt die folgende Beziehung für b, aufstellen; 
124, + 4,478 2 


64, +02 (3A; —3-+b,) 
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Ss 
x 
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Diese Gleichung für b, muß in jedem Fall numerisch gelöst werden. Da man den Bereich 
kennt, in dem 5, sinnvollerweise liegen muß, kann man die Gleichung ziemlich schnell durch 
Probieren lösen. 


4. Herstellung zweier. Kurvenblätter als Grundlage des Verfahrens 

Nachdem die Profilklasse festliegt, kann ein Teil der Arbeit, die zur Auflösung der Diffe- 
rentialgleichungen (20) und (26) bzw. (28) geleistet werden muß, unabhängig vom besonderen 
Beispiel (d.h. unabhängig von der Form der Wand oder von der vorgegebenen Außenströmung) 
ein für allemal ausgeführt werden. Dabei ist eine Kleinigkeit vorweg zu bemerken. Unsere 
Gl. (20) und (26) bevorzugen die Impulsverlustdicke 9, und die Enthalpieverlustdicke ®;, 
während in unserer Profilklasse die Grenzschichtdicken ö und ö, bevorzugt werden. Zwischen 
diesen beiden Größen besteht Proportionalität. Es ist nach Gl. (17), (18), (22) und (23) 

1 


9= 5 | (I) am 8,7, 08, 24 84/) 


0] 
1 
9,= 0, |" 1 an = 81J2(M, 1, 0,/0) 
() 
0 
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Wir haben in jedem Fall n, als unabhängige Variable gewählt, weil 1.a. ö, > ö ist; der Inte- 
grand verschwindet dann evtl. vor Erreichung des Wertes 7, — 1. Die Größen J,J2 J3 Ja 
müssen als Funktionen der angegebenen Argumente A,, M und ö,/ö ausgerechnet werden. Es 
empfiehlt sich, diese Werte ebenso wie die Parameter A, und A, zunächst in die Gl. (20) und (28) 
einzuführen. Es ergibt sich für Gl. (20) 


da RR Da (2%) ' = 
aa 2: Ad 2J;, +J,8—M?)] +2 al ENG: (42). 


At an der Gl. 08) brauchen wir die Größe en, dis sich in den 6 
Ä yes M folgendermaßen U 
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Die in dieser An mitgeteilte ein für allemal a Vorarbeit zur Integration derDitter 


Be B) e ae 
Be rentialgleichungen (20) und (28) besteht in der Berechnung der Größen 2( (&). und‘ £ 
rk 


En e (Kurvenblätter 3 bis 6) und der kuryenmäßigen Darstellung der Größen Kr und 
rn  » Leider ist die Zahl der eingehenden Parameter bei unserem Problem weitaus größer als 
im inkompressiblen Fall und zwingt zu einer gewissen Auswahl unter den Fällen, deren Be- 


we .: rechnung an und für sich wünschenswert wäre. Wir haben uns für den Fall interessiert, daß 
Bu die Wand durch Heizung oder Kühlung auf einer festen Temperatur gehalten wird, und zwar 

Be: haben wir für die numerische Behandlung zwei Fälle gewählt. Der erste Fall soll etwa bei sonst Re: 
+ gewöhnlichen Luftverhältnissen einer Aufheizung der Wand auf.100° C entsprechen, der zweite Br: 
f- | 

er ne [3 w=05 

“ o 
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as fie Machzahl RE so müssen wir im reinen u Überschallhereich bleiben, ı und habeı 2 
s im ı ee auf ‚den Bereich > A EN AR N be 


ET f BA N 1S<M<2 ER RR ö 
er beschränkt mit einigen ana für ‚größere N ne Ma c en sehen uhr - 
3 = Ferner haben wir gesehen, daß sowohl in den Gleichungen für die J als auch in ‘den 
(42) und (44) die Größe E es Parameter auftritt. pa —. sich gezeigt, daß es Be 
den Bereich Kr T Se; BE e 
AN TS, 3% 3 EEE IEN et ri BB. 
zu nBerdbcken mit einigen Stichangaben® außerhalb des Bereiches. Der Po 0 " l he ausen- 
parameter überdeckt den Bereich von A, — 0 bis zuf Ablösung, also etwa bis er — B* 20 im 
Balls I Ye a = — - im ale I. Aue 
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m=b„-m„ (M- 15) 
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Die Größe dz/dz ist 
in den Bildern 1, — in 
denBildern2, J,undJ, 
in den Bildern 3 und 4 


und 2(&*) und >) 
ön/w On 


in den Bildern 5 und 6 
niedergelegt. Die in 
der Darstellung be- 
nützte Linearisierung 
gilt im obengenann- 
ten Parameterbereich 
—20<A,<0, 15<M 
<2 und 1< ö,/ö< 1,5. 
Eine einzige Größe 


(&) für denFautt) 
on ” 


ließ sich nicht in der 
angegebenen Weise li- 
nearisieren. Für sie 
werden deshalb nicht 
die partiellen Ableitun- 
gen, sondern die Größe 
selbst angegeben. ‚Alle 
Größen sind für die 
Fälle I (?„/i,=1,3) und 
IL (&o/t, =.0,85) _dar- 
gestellt. Für den Ex- 
ponenten des Zähig- 
keitsgesetzes (7) wurde 
im Falle Io= 0,5, im 
Falle II = 1gewählt. 
Die Prandtlsche 
Zahl wurde immer mit 
co —= (0,7 Angesetzt. Um 
den Einfluß des Zä- 


und © = 0,5 (für dieselbe Heizung Gy/%o = 1,3) zusammenfallen, womit die auch von andere 


 (s. auch die Bemerkung zu [5] in der Einleitung). 
Zeichen für die Zweckmäßigkeit der Fassung (28) für den Energie-Integralsatz. (In einer zu- 
nächst ausprobierten Fassung fielen die Kürven für = 0,5 nicht mit denen-für o=1 zu- 

. sammen). Die Bilder 1, die zum Impulsintegralsatz gehören, hängen hingegen vom PREIS 
gesetz ab; wie es nach den Ergebnissen von [5] sein muß. r 
Pe Mit Hilfe der Bilder kann nun die Integration sehr EA geleistet. erden Die) 
Gl. (20) wird als Gleichung für die unbekannte Funktion 2 oder auch ö und die Größe ö,/ö zu- 
- nächst als Parameter aufgefaßt. Ebenso wird die Gl. (28) als Gleichung für 2, oder ö, und 6,/6 
zunächst als Parameter aufgefaßt. Da man auf jeden Fall schrittweise integrieren muß, ist ' 
es zweckmäßig, abwechselnd einen Schritt im Lösungsverfahren der Gl. (20) und einen indem 


Schritt. 


Ve ild) h: ur m 7 (& Be 
wiedergegeben) becedttmeh. Dabei hat sich gezeigt, daß für de Größe —r die a), gie oo 


Autoren behauptete Unabhängigkeit des Energiesatzes vom Zähigkeitsgesetz bestätigt wird 


der Gl. (28) zu machen und das berechnete Verhältnis ö,/6 dem nächsten Schritt zugrunde 


zu legen. Dabei hat sich ergeben, daß das Verfahren nicht konvergiert, wenn man den Schritt k 


Diese Unabhängigkeit ist umgekehrt ein 


n-+ 1 mit dem beim Schritt n gewonnenen 6,/ö berechnet, vielmehr muß man dem Schrittn+1 


den Mittelwert aus den beiden 6,/ö zugrunde legen, die beim Schritt n einerseits zugrunde 


gelegt und andererseits herausgekommen sind. Eine Iteration war dann nicht mehr nötig. Es | 
ist bemerkenswert, daß auch das Crocc.osche Iterationsverfahren [5] zur Bereehnung der 


Schubspannung nur konvergiert, wenn beim zweiten Schritt nicht das Ergebnis des ersten 


Schrittes, sondern das Mittel aus dem Eingangswert und dem Ergebnis des ersten Schrittes als 


neuer Eingangswert für die Schubspannung genommen wird. 


Im einzelnen sind wir so vorgegangen. Die Werte der Außenströmung wurden alle über 


der Machzahl der Außenströmung aufgetragen, das erleichtert die Interpolation. Man beginnt 


mit einem Anfangswert von 2 (etwa = 0 an der Stelle x,), liest dz/dx ab, gewinnt durch Multi-. 


plikationen mit Az ein z an der Stelle %, berechnet daraus A und ö. Dabei braucht man außer 
den Kurvenblättern 1 noch die Werte J,. Zur Gewinnung von 2, liest man —r ab, multipliziert 
mit dem zweiten Faktor unter dem Integral und integriert graphisch. Aus 2, ergibt sich ö,. 
Hierbei braucht man J;. Mit dem schon genannten Mittelwert von 6,/ö, dem gewonnenen 
Parameterwert A, und der Machzahl des nächsten Schrittes berechnet man den nächsten Schritt 
wie den vorhergehenden. Dabei wächst die Kurve des Integranden in (28) von Schritt zu 


5. Übertragung des Verfahrens auf den rotationssymmetrischen Fall 
WieMangler[1ö]gezeigt hat, kann man die Grenzschichtrechnung an einem Rotations- 


körper auf die Grenzschichtrechnung an einem zweidimensionalen Profil eines in Spannweiten- 


richtung unendlich breiten Körpers zurückführen. Man hat dazu nur die folgenden Transfor- 
mationen auszuführen, wo sich die Größen mit Querstrich auf das ebene Ersatzprofil beziehen. 
Es ist 


P(s)= pe) : 

WE) -ile) -_ [ r(8) 443 
e(s)=eol) nl L? Di en: 
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sist die Bogenlänge des Meridianschnittes des Rotationskörpers, r,ist der Abstand eines Punk- 
tes der Oberfläche von der Achse, so daß die Kontur durch r = r, (s) gegeben ist. L ist irgend- 
eine Bezugslänge. s ist die Bogenlänge der ebenen Ersatzkontur. Dem Punkte s wird dann 
der Druck bzw. die Geschwindigkeit im Punkte s zugeordnet. Wenn also die äußere Strömung 
im rotationssymmetrischen Fall bekannt ist, dann ist durch (45) auch die äußere Strömung 
im zugehörigen ebenen Fall bekannt, und die Grenzschichtrechnung funktioniert wie bisher. 

Aus den Eigenschaften der Transformation (45) kann man nun schon einiges ablesen. 
Wir wollen uns überlegen, wie die zu einem Meridianschnitt r, (s) gehörige ebene Kontur aus- 
sehen wird. Für die Grenzschichtrechnung interessant sind Konturen von der Art, wie eine 
in Bild Ioben links angedeutet ist. Die Strömung um solche Konturen muß im Überschall 
einen Druckanstieg überwinden. Der Druckanstieg wird um so kräftiger sein, je dicker das 


- Profil ist. Nehmen wir an, daß er kräftig genug ist, um in ebener Strömung Ablösung zu er- 


zwingen. Jetzt fassen wir dieselbe Kontur als Meridianschnitt eines Rotationskörpers auf. Dann 
ist die Kontur nach (45) zu transformieren. Dabei bleibt der Integrand (r,) über ein sehr großes 
Stück der Achse sehr klein und fängt dann allmählich an zu wachsen. Das zum gegebenen 
Meridianschnitt gehörige ebene Profil ist also viel dünner als das Profil der Abb. I. Es ist also 
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Be iehkezn erwarten, daß am Rotationskörper mit dem ] 


auftritt. Er ae yes i Ä serie i £' he fi ” b :., 
Die Schwierigkeit liegt für diese Fälle nicht in der Grenzschichtrechnung, so n 
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Berechnung der Potentialströmung. Eine solche Rechnung führt man nach | harakte- 
ristikenverfahren durch, entweder nach dem sogenannten Felderverfahren, von denen das be 


kannteste das von Prandtlund Busemann ist, oder nach dem Gi 


aufwand bleibt erträglich. ’ 
6. Die Beispiele h h 
Da es bei der Auswahl der Beispiele auf mögliche Anwendungen des Verfahrens nicht 


ankam, haben wir Beispiele gesucht, die für die Berechnung der äußeren Strömung möglichst — 


günstig lagen, die einen vernünftigen Druckanstieg bis in die Nähe der Ablösung liefern und 
dabeietwa den Machzahlbereich 22 M> 1,5 durchlaufen, Mit Mach zahlen größer als 2 zu 
beginnen, war im Hinblick auf etwaige Vergleiche mit Messungen unzweckmäßig. Mit Rücksicht 
auf die Berechnung der äußeren Strömung waren Verdichtungsstöße in der ankommenden 
Strömung vor ihrem Auftreffen auf den Körper zu vermeiden. Wir haben also den Körper so 
gewählt, daß die Strömung tangential zur Wand ankommt. Für die ebene Strömung kann 
man dann nach den Untersuchungen über die Strömung an einer konvexen Ecke annehmen, 
daß Druck, Geschwindigkeit und Temperatur durch den Neigungswinkel gegeben sind. Wir 
haben also eine Kontur 
y=a(l—cos x) 
gewählt. Dann ist der Neigungswinkel y gegeben durch 
tgy=asinz 
und wir können zu diesem Neigungswinkel einfach aus einer Tabelle (etwa in [18] S. 57) Druck, 
Geschwindigkeit, M ach zahlund Temperatur ablesen und die Dichte berechnen. Dabei können 
wir den Nullpunkt der y-Zählung so wählen, daß wir die Machzahl M= 2 dem Winkely=0 
zuordnen und haben es nun durch Wahl der Konstanten a in der Hand, wie wir den Bereich 
2> Mz 1,5 über das Profil verteilen wollen, d.h. wie dick wir das Profil wählen. Wir zeigen 
hier ein Profil, für das die Ablösung kurz vor dem Wendepunkt liegt (s. Bild I oben). Für 
diese Profilkontur wurden nun die Fälle I und II durchgerechnet. Im Falle I ist also das 
Profil etwa auf 100° C geheizt, im Falle II auf —25° C abgekühlt. 
In den Abbildungen ist die dimensionslose Wandschubspannung 
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aufgetragen. Um den Vergleich mit den Rechnungen an der Platte zu ermöglichen, haben wir 
c noch mit der Wurzel aus einer Reynoldsschen Zahl multipliziert. Für die ‚Platte in 
inkompressibler Strömung ist bekanntlich & YUxo/u konstant. Wir haben aufgetragen 


ER /7 Re IRr =x 
Rn =) U00d _ jo es a (50). 
eg en ww Ykes | 


Weiter definieren wir eine Wärmeübergangszahl & durch die an der Wand übertragene Wärme- 
stromdichte gq | 


CG= 


oT 
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die wir so definiert haben, daß sie mit der vonHantzscheund Wendt [3] eingeführten 

Definition vergleichbar wird (dort tritt allerdings die Differenz aus Thermometertemperatur 

und Wandtemperatur auf). Weiter definieren wir einen dimensionslosen Beiwert 6% durch 
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 ponenttow=1 gewählt wird. Für andere Zähigkeitsgesetze © 1 besteht eine geringe Ab 


gigkeit ( der Größe 90) von Machzahl und Temperaturverhältnis. Jedenfalls ist c, YRe wiede 
eine Größe, deren Verhalten an der ‚Platte einen reinen Kompressibilitätseffekt darstellt. Wen 


von 9(0) vergleichen, so haben wir einigermaßen den Einfluß des Drnekanstieges in unserem E 
i Beispiel erfaßt. WAT tragen also auf x 
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g(0) ist unabhängig von der M ach schen Zahl und vom Teinperatehverhältnie, eng. d der 


wir Ca Y Res, in unserem Beispiel mit den von HantzscheundWend D gegebenen Kurve, 
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Wir betrachten zunächst Bild I, das sich auf die geheizte Wand der oben im Bild’ 
angegebenen Form bezieht. Man sieht, daß die ebene Strömung noch vor dem Wende- 
punkt ablöst, denn die Wandschubspannung vom Zähigkeitsgesetz ist sehr gering, äußert sich 


eigentlich nur darin, daß für die größeren Machzahlen die aufgetragene Größe im rechten Bild | “ $ 
etwas kleiner ist. Da das dem rechten Bilde zugrundeliegende Zähigkeitsgesetz dem Tempe- 27 
raturbereich besser entspricht als die einfache Annahme u proportional 7 (linkes Bild), so kann Br 


man sagen, daß das einfache Gesetz etwas zu große Wandschubspannungen liefert. Nach den " 
in [3] mitgeteilten Kurven g(O) zeigt sich dieser Einfluß des Zähigkeitsgesetzes besonders bei a 
1-0 a 
den größeren Mach zahlen.. Auch der Faktor A wirkt für © < 1 verkleinernd auf «Ye. = 
Für den Einfluß des Zähigkeitsgesetzes auf die Wärmeübergangszahl gilt dasselbe. Daß die rl 


-dimensionslose Wärmeübergangszahl c, längs der Wand abnimmt, liegt einfach am Anwachsen 


der Wärmegrenzschichtdicke ö,. Wenn die Temperaturgrenzschicht dicker wird, muß die 
Wandneigung des Temperaturprofiles flacher werden, d.h. die Größe (97/d y), muß ab- 


nehmen, während ( (ailamı)o nahezu konstant bleibt, wie es in unserem Falle I auch zutrifft 


(die Differenz 4 1= = — 1 nimmt längs der Wand ab, weil die Mach zahl abnimmt). 


Wir kommen zu Bild II, das die Entwicklung der 'Grenzschicht an demselben Profil 
zeigt, wenn dieses Profil auf —25° C abgekühlt wird. Hier ergibt. sich, daß die Schubspannung 
durch Heizung wenig zu beeinflussen ist. Nach den Ergebnissen der Rechnung an der Platte 


mußte das erwartet werden, weil der Faktor g(0) in der Formel für « YRe, wenig von der Hei- 


1-0 
zung abhängt und der zweite Faktor 4, 2 für die üblichen ®-Werte auch sehr klein ist. Die 
Ablösungsstelle hängt kaum von der Heizung ab. Die aus der Gl. (37) en Ablösungs- 


bedingung A, < — 1240 
liefert zwar für die geheizte Wand viel größere |A,|-Werte als für die gekühlte Wand. Das 
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Sen ERLTE A EN 
" kommt daher, daß für die geheizte Wand die Grenzschichtdicke 6 größer ist als für die 
und sagt über die Ablösungsstelle nichts aus. Daß die Grenzschichtdicken durch Heizung 
wachsen, zeigen auch die Ergebnisse der Rechnungen an der Platte. Das Bild II scheint dafü r 
zu sprechen, daß die Ablösung durch die Abkühlung etwas hinausgezögert wird: Esel: 


Wr 


ö,/ö unseren vorbereiteten Bereich, worunter die Rechen- 
genauigkeit gelitten hat. Grundlegend verändert ist. durch 
die Abkühlung das Verhalten der Wärmeübergangszahl. Ob- 
wohl die Wandtemperatur kleiner ist als die Ruhetemperatur, 
ist sie doch noch größer als die bei den hohen Mach zahlen, 
sehr niedrige Außentemperatur. Die Wand gibt also Wärme 
an die Strömung ab. Infolge des starken Anwachsens der 


Temperaturgrenzschichtdicke wird wie im Falle I c, / Res, zu- 
nächst kleiner. Man sollte erwarten, daß diese Abnahme 
in größerer Entfernung von der Profilvorderkante immer 
deutlicher wird, weil die Temperaturgrenzschichtdicke dort 
stark anwächst, und weil für große Verhältnisse ö,/ö die Nei- 


gung der Temperaturprofile'(ö:/2n,), sehr klein und unter 
Umständen sogar negativ wird. (Es treten dann Werte 
Blaım iin > iolia auf, solche Profile siehe in [3].) Wie das Bild II 
N zeigt, steigt aber die Zahl Res, nach dem anfänglichen 
Absinken wieder an. Das ist ein Einfluß des mit den kleineren Machzahlen abnehmenden 
Unterschiedes von Wandtemperatur und Strömungstemperatur, aber auch ein Einfluß der in- 
folge des Druckgradienten (trotz ansteigender Temperaturgrenzschichtdicke) steiler abfallender 


dimensionslosen Temperaturprofile (|97/ön,|» wird tatsächlich größer). Dazu geben wir fol- 
gendes zu bedenken. Allgemein ist einzusehen, daß nach der Ablösung hin, wo die v-Kompo- 
B nente der Geschwindigkeit (insbesondere 9v/dy) wächst, der Wärmeaustausch an der Wand 
RB größer werden wird. Bei ohnehin sehr schwachen Wandneigungen des Temperaturprofiles 
{ (gekühlte Wand und große Temperaturgrenzschichtdicke) ist die prozentuale Änderung infolge 
| des Druckanstiegs besonders groß. Es tritt also nicht so leicht der Fall ein (in unserem Beispiel. 
Br, jedenfalls nicht), daß die Wärmeübergangszahl ihr Vorzeichen wechselt, daß also die Tempe- 
ratur in Wandnähe höher ist als an der Wand, wenn die Außentemperatur niedriger ist als die 
Wandtemperatur (solche Fälle treten in unseren vorbereitenden Rechnungen im Fall der ge- 
kühlten Wand bei kleinem A, und großen ö,/ö auf). Beisehr starker Kühlung und großen Mach- 
zahlen würden sie auch bei größerem |A,| auftreten. Es ist aber darauf aufmerksam zu machen, 
daß unser Ansatz für das Temperaturprofil zur Erfassung solcher Fälle nicht gut geeignet ist. 

Die Überlegungen des vorigen Abschnittes lassen nun schon erwarten, daß wir von der 
Ablösung weit entfernt sein werden, wenn wir das Profil des Bildes I als Meridianschnitt eines 
Rotationskörpers wählen, wie es die Rechnung auch zeigt. Wir haben das Profil dann auf das 
Doppelte verdickt, und noch immer ist an Ablösung nicht zu denken. Das Profil ist aber 
schon so dick, daß man nicht hoffen kann, durch weitere Verdickung noch etwas Vernünftiges 
zu erreichen, insbesondere ist es nicht möglich, für ein noch dickeres Profil mit der vorliegenden 
Methode [16] die Außenströmung zu rechnen. Schon die Strömung um das auf das Doppelte 
verdickte Profil zeigte einen Verdichtungsstoß, dessen Behandlung im rotationssymmetrischen 
Fall einige Arbeit gemacht hat. Natürlich sind Körper denkbar, an denen Ablösung auftreten 
wird, ihre Entwicklung auf dem von SchäferundTollmien[16]gewiesenen Wege würde 
wohl möglich sein, Fonnte aber nicht das Ziel dieser Arbeit sein, die sich mit Grenzschicht- 
rechnungen befassen wollte. Deshalb wurde die Untersuchung hier abgebrochen und die offen 
gebliebenen Fragen einer weiteren Untersuchung vorbehalten. 


D 


7. Zusammenfassung 


Zur Berechnung der Grenzschicht in kompressibler Strömung ist außer der als Grenz- 
schichtgleichung bekannten Impulsgleichung noch eine Energiegleichung zu berücksichtigen 
Durch Integration dieser Gleichung über die Grenzschicht gewinnt man eine gewöhnliche Diffe- 
rentialgleichung, die man als Energie-Integralsatz bezeichnen kann. Diese Gleichung tritt zu der 
integrierten Grenzschichtgleichung (Impulsintegralsatz) hinzu. Ein Pohlhausen verfahren 
im Überschall läuft also auf die Integration zweier gekoppelter gewöhnlicher Differentialglei- 
chungen. hinaus, die man beide in die Form von linearen Differentialgleichungen bringen kann 
Das Integrationsverfahren, das schrittweise abwechselnd die eine und die andere Gleichun 
befriedigt, wird beschrieben. Kurvenblätter für die beiden Fälle einer geheizten und She 


\ 


Aussage wäre voreilig, weil die Rechnungen im Fall II in der Nähe der Ablösung unsicher sin. 
Dort überschreitet das Verhältnis der Grenzschichtdicken 
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ZurFehlerschätzung bei gewöhnlichen Differentialgleichungen 


erster Ordnung 
Von @eorg Hamel in Berlin 


Verf. gibt eine Rekursionsformel für die Koeffizienten einer majoranten Reihe der Lösung einer gewöhn- 
lichen Differentialgleichung erster Ordnung, die zur Abschätzung des Fehlers beim Abbrechen der Reihe dienen 
kann, und einige damit zusammenhängende Formeln. Die Formeln wurden durch Zurückgehen auf eine 
partielle Differentialgleichung gewonnen. 


The author gives a formula of regression for Ihe coefficients of a series that is always major than the so- 
lution of an ordinary differential equalion of the first order, and in connewion with it, some further formulae. 
The formula that serves for the valuation of the error, üf Ihe series is cut off, is oblained by regression to a 
partial differenlial equation. 

L’auteur donne une formule regressive pour les coefficients d’une serie, qui est la majorante dc la solution 
d’une öquation differentielle ordinaire du premier ordre et qui peut servir pour l’evaluation du erreur, quand 
le developpement est rompu, en outre il donne quelques formules rapportantes. Toutes. Toutes les formules 
sont derivees d’une Equation differentielle ‚partielle. 


Asrop naeT pekypCcuBHy® dopMmyıy AA KOabmmmeHToB PANa, MASKOPAHTHOTO OTHOCH- 
TeI5HO PeINeHHA IPOCToroO AuhepeHNHansHoro YPaBHeHus epBoro HOopAnka. YTa Popmyıra 
CHY:KHT IA OMEHKH IOTPEINHOCTH CYMMEL KOHEYHOTO YMCHA WIeHOB PANA, KPOMe TOrO 
Aal0Tca HEKOTOPHE APpyrue oTHocAmmecH cıoya hopmyısı. Bce dopMmyısı Ösimm MONyYeHBL 
HpH moMoImm nubepeHlumanbHoro YPABHeHHA B YACTHEIX IIPOHSBONHEIX. 


Konvergenzbeweise sind beim Praktiker höchst unbeliebt, sie wären es nicht, wenn Bücher 
und Professoren mehr Nachdruck auf die so notwendige Fehlerschätzung legen würden als auf die 
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Diese Gleichung hat die für &= 0 verschwindende Lösung 
n=1-—Y1+2a log 1—& mit 


Und nun begnügt man sich gewöhnlich mit der Feststellung, daß diese Lösung im Kreise 


5? 
lii<1-e ” 


regulär ist, die wahre Lösung um so mehr. Aber viel wichtiger wäre die Berechnung der majoranten 
Reihe 


Denn mit ihr hat man die Kikenntnig) daß die N RL der wahren Reihe kleiner als die C, 


sind, Das ist aber das, was die gewonnene Potenzreihe überhaupt erst brauchbar macht. Denn 


damit übersieht man den Fehler, den man macht, wenn man die Reihe abbricht. 

Nun scheint es kein Problem, das den reinen "Mathematiker reizen könne, die Reihe für das 
majorante n aufzustellen; er wird sagen, das seitrivial. Aber vielleicht hat es doch einiges Interesse, 
zu erfahren, daß man eine sehr bequeme Rekursionsforme] angeben kann, die gestattet, die Koeffi- 
zienten fast so schnell zu berechnen wie die Binomialkoeffizienten. Auch kann man aus ihr eine 
sehr einfache Abschätzung der C, finden (Formel 9). 


Die Formellautetso: 
Man erkennt leicht, daß C', die Form hat 


C N? < Kna x SM 

iur, nal eh Er Er 

und nun gilt h 
Kn,u = (n — 1) Ko-ı.a +(24 —3) Kar-ı za DEREN TIITE (5), 
Kr, =( für A<n, Kıı =1. 

Als Sonderfälle erkennt man sofort 


Kun=m-ll,. Enno (n>2) 
Man berechnet schnell die Tafel’ 


1 1 = _ _ _ _ = 
2 1 1 = = dit = = 
3 2 3 3 = 3 = = 
4 6 11 18 15 ae = = 
5 24 50 105 150 105 a “= 
6 120 274 675 1275 1575. 945 N 
7 720 1 764 4 872 11.025 18 375 19 845 10 395 
8 | 5040 13 068 39 396 101535 |. 205 800 347 040 291 060 


z.B. berechnen sich die K,a aus Kr =6Rgı +(2A — 3) KRsı-ı 
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n ) a ie Sn Herr Dir Dr.-Ing. 
se Koeffizienten berechnen und fand empirisch die oben angegebene Rekursionsfor 
Sie auf diesem Wege beweisen zu wollen, scheint ‚kaum möglich. Der Eingang in die- 
gleichung selbst führte etwas weiter, nämlich bis zum Beweis der, Formellür Kun 
Nun gelingt aber der Beweis der Krackeschen Formel ( 5) durch Rückgang : auf par- 
A Ditferentialgleichung, der n als Funktion von hi und Runen Es ee der . der für ähn. 

e ‚liche Aufgaben Bedeutung haben dürfte: a 


Satz: Die Funktion tr er JR a 
f Ren se[itefzoe] $% a a 


Br Ba = 7 der Quotient zweier Potenzreihen mit 10) +0. sei 


und m eine Set Zahl, genügt der, linear homogenen partiel- 
len ee Ehe we Ton 
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ES : ? \ ' Hal \ a Er 
a + Bla er AR ARE - N, 
mit der ed N a u 2 Sind A und Baer ratio- 
na]l,"so 'hat.dTe ‚Differentialgleichung ganze rationale Koeitase 
en und es gibt eine ähnliche Rekursionsformel wie dieoben 
im Sonderfall angegebene Formel (5). Unsere Sonderaufgabe (2) ist enthalten, 


ae =— og setzt, Hier heißt die partielle Differential- 


© wenn mn1—7=3,2=&m 22: 


% gleichung ! 
Br: j on { y . 
| —_ 2 — = a att ea kat art, SE Ma ee 
Be za. SA—2Rz tan ma N 
| VI ist als 1 zu nehmen. Se 
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> 2 m—1 R 
=? % _m|i+af Ra)da] Ras ie, 
Er R ; 1rajR) de 
3 2 ® AN mz|R(a) da 
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Ir (®) )da, so bekommt man sofort die partielle Differentialgleichung (6). In die Gleichung (7) 


Gehen wir nun mit den Ansatz (3) hinein und WR - = 
BR, tee + &n Der=a. 


Für N folgt 
N NE OS ee 
- Vergleich von &*-1 aber gibt für n>1 
ne 
nn (n—]1) On — a0ıı +20 I. x 
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4 = n—1 einsetzen und erhält 


Ks 3 29) Ki Pr 


RN R,=1, K,=3 ist, folgt K, <3-2- Sand aa 
ae! \ 
sand. 
Daher ist der Koetfizient selbst WR 
und 

Er On< ar falls a+list 

| are in ne a ; 
Die wahre Reihe konvergiert also sicher besser als i 

| .n v2 Saal s 
+ a A 


womit eine Fehlerschätzung von selbst gegeben ist. 
‚Im allgemeinen Fall der partiellen Differentialgleichung (6) wird man ebenso vorgehen. 


SeiA= Saw mit ,+=0, B= 3 bux", so ergibt der Ansatz 
‚0 - - 


n>0 
dagegen für A=1 
a, Kn,ı & bu Kn,o 


nen +l (nl) 
n>0 Nn>0 


= mby. N =0, 2883 


2=14+20ne a ER I (11) 
(wir können ohne Beschränkung der Allgemeinheit 1” — 1 setzen) 
ao’ ZnOı ar +&bum da” Fe ar ma ci) = ma I bua* 
0 
und daraus folgt durch Vi ergleich gleich a Potenzen von x i 
dc 
> nyOat+ 3% | . a — ma0,| = maby, N=0,)1,2%: 
n+v=N+1 n+u=N da Ri 
n>0 n>0 
Insbesondere folgt für N = 0 
also 
Q, = ul a 
Setzt man nun ’ 
4 < Ka 2 
so bekommt man weiter durch Vergleich der Potenzen von «a? 
Qy. En. Kn, = 
+ b Al Se ü a 
Bzeae (n—1) DR eh m) =D für Az . (12), 


Dirads ee man alle er 1 N Die allgemeine Formel (2) n N > i gestaltet dann 
Berechnung der Kn,ı für a Ve 


gr Sind A, B ganz rational, so bekommt man Bekursiorsfonmein a: ie Gliederzahl, 

Be eu wieder einfach, wenn A und B nur wenige net nl SER 

RATE Beispiel: Soll ES PER OR LIE fa RR 
“ j en + aarctg z)m SR RE 2 


Bonickett werden, so setze man A=1+ 2, Br —1 und hat dann _ Sormel) 
| dl | Re 
d+9) Srommre Pr ma Zoat -ma Kar 


und somit für » = 0 | EL a Ma 


‚sonst dc F | | Ber 
| “ at On u: (m Ge y Cn+1+ a? Ey „=0, n>0 (G,=P0). 
, ; Kıa Bu H \ = 
Nir.C, = a folgt 
wenn a 
n+1,A eo A ee ; NER ; ih 


Für A=0 ergibt sich daraus 3 RER BR 
Karı, o=— Kn-1,0 nn — 1) » 


was mit Koo=0 und XK,,o = alle Kt —() ergibt in Übereinstimmung mit einer obigen 2 
Bemerkung. Allgemein folgt a 
Kari = — nn — 1) Kn,a—(A—1—m) Knı-ı | 
als Gegenstück zur Krackeschen Rekursionsiormel. 
Eingegangen am 28.2. 49. 
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Beiträgezurnumerischen Lösunglinearer Eigenwertprobleme. I’ 
Von N. Joachim Lehmann in Dresden 


Das Ziel der Untersuchungen. ist die Eingrenzung eines beliebig vorgeschriebenen ee linearer 
selbstadjungierter Probleme; dabei soll der Rechenaufwand für jede verlangte Genauigkeit möglichst gering 
bleiben. Die Lösung der :Aufgabe gelingt mit einer neuen Bigenwertdarstellung, jeder interessierende Eigen- 
wert erscheint dabei einmal als Minimal- und einmal als Maximalwert einer Integralform. Passende Ver- 
gleichsfunktionen ergeben daraus in einfachster Weise die ersten Eigenwertschranken. 


The investigalion aims at the limitalion of an eigenvalue of Linear self- adjoined problems, the order 
of which is arbitrarily preassigned. The amount of compuiation work shall remain as small as possible 
for any required exaciness. T’he problem is succesfully solved by a new representation of the eigenvalues: 9 
Each fixed eigenvalue vs represented first as minimum secondly as mawimum of an integral form. Suitably Yu 
chosen ‚‚neighbouring functions‘ yield the first boundaries of Ihe eigenvalue in a most simple way. 


Br a ne Di ze 


C’est le but.de cette-recherche de donner des limites pour une valeur propre, d’un ordre arbitraire, de 
probl&mes Tineaires identiques a leurs adjoints; pour atteindre toute approximation desiree on veut calculer 
le moins possible. On arrive a la solution de ce problöme par une representation nouvelle des valeurs propres: 
La valeur propre en question se pr&sente. une fois comme le minimum, l’autre fois comme le maximum d’une 
forme integrale. On y trouve, d’une maniere tres simple, da l’aide de "fonctions approwimalives convenables, 
les premieres limites des valeurs propres, 


Naunoe uecenoBaHne CTABUT cebe HEeNbE 'MATb ABYCTOPOHHWP OIEeHKy MIA JIOÖBIX 
COÖCTBEHHEIX YUCeAI HPOÖJeM, MubepeHmmasıbHoe BEIPAIKeHHE KOTOPEIX COBIANACT C „‚AnbIOH- 
THPOBAHHEM‘‘ nuhepeHunanbHeIm BbIpaskenuem. Ilpı 5ToMm KonnyecrBo pacuerHoä pabork 
aa 10004 Tpeöylmelch CTeIeHNH TOUHOCTH AOMKHO ÖBITb IIO BOSMOSKHOCTH MEHbIIe. ITA 
sanaya penmaerca IIPM IIOMOIIH HOBOTO BBIPaKeHHA MIA COÖCTBEHHEIX UHCEN! KAsKII0E HCKO- 
MO COÖTBEHHOE YHCHO IPEeNCTABIAeTCH CHAYANIA KAK MUHHMAJBHO®, & BATEM KAK MAKCH- 


1) Diese Arbeit stellt den Teil I der von der Techn. Hochschule Dresden 1948 genehmigten Dissertation 
des Verfassers dar (Referenten: Prof. Dr. Fr. A. Willers und Prof. Dr. O.-H. Keller). 
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342 Lehmann, Beiträge zur Lösung E 2 
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MaIBHO6 SHAUCHME HEKOTOPOrO HHTETPANBHOTO Burpamenns. Ilpm momomu Nomxons 
NOLYCTUMEIX DyHRImKa IOIyaalTca OTCINa HamÖONee IPOCTEIM O0PasoM HepBEle ABYCTOP 
OICHKEH IA COÖCTBEHHEIX YHCEN, nF, s 


Einleitung und Zielsetzung al, 

Die meisten in der Praxis auftauchenden linearen web { 

zwingen zu einer näherungsweisen Berechnung der Eigenwerte. Von der umfa ichen 
darüber sollen hier nur die zusammenfassenden Darstellungen von Hohe nemser?) und 
Collatz°)®) erwähnt werden, dort finden sich auch ausführliche Schrifttumsverzeichnisse. 
Es zeigt sich in diesen Arbeiten, daß es bei vielen Aufgaben verhältnismäßig leicht ist, ur 
Eigenwertschranken zu bestimmen — z.B. mit Hilfe von Vergleichssätzen®)*) oder den Ein- 
schließungssätzen von Temple®)#)5)®) und Kryloff?)®). Ebenso bekommt man mit 
erträglichem Rechenaufwand gute obere Schranken für die Beträge der absolut genommen 
kleinsten Eigenwerte; man nutzt hierbei die klassischen Extremaleigenschaften der Eigenwerte 
(evtl. zusammen mit Iterationsschritten) zur Bildung von Minimalfolgen aus (Ritz- Galer- 
kinsches Verfahren). Diese Methoden liefern auch obere Schranken für die höheren Eigen- 
werte, haben dann aber den Nachteil, daß die Rechenarbeit außerordentlich anwächst (ein 


' Überschlag ergibt, daß der Arbeitsaufwand auch bei günstigster Rechnungsführung wesentlich 


stärker als quadratisch mit der Ordnungszahl des interessierenden Eigenwertes ansteigt). Die 
größten Schwierigkeiten ergeben sich jedoch bei der Frage nach einer entsprechend brauchbaren 
unteren Schranke Die bisher bekannten unteren Eigenwertschranken finden sich schon fast 
alle im Eigenwertbericht von Collatz, die Situation hat sich seit 1939 kaum geändert (zu- 
mindest nicht nach der mir zugänglichen Literatur). Allen Anforderungen an Genauigkeit und 
Einfachheit genügt dabei nur eine Formel von Temple°), die für einen einfachen kleinsten 
Eigenwert eines positiv-definiten Problems eine untere Schranke bestimmt; bei Collatz, 
Eigenwertberechnung, wird diese Beziehung als ‚„„‚Hauptformel‘“ benutzt. (Eine von Col- 


‚lat z5)10%)11) gegebene Erweiterung dieser Schranken für höhere Eigenwerte ist zwar mathe- 


matisch interessant, aber praktisch kaum brauchbar.) Neu ist schließlich die Rellichsche 
Fehlerabschätzung für das Störungsverfahren !?); jedoch sind auch diese Schranken sehr grob, 
die Näherungswerte sind meist um mehrere Ziffernstellen genauer als die Abschätzungen ver- 
muten lassen. 

Mit den folgenden Arbeiten sollen diese Schwierigkeiten beseitigt werden. Dazu wird 
jeder Eigenwert einmal als Minimal- und einmal als Maximalwert (im Gegensatz zur Minimum- 
Maximum-Darstellang) einer quadratischen Integralform dargestellt. Daraus wird als erstes 
ein Einschließungssatz für die Eigenwerte entnommen, der oft eine überraschend gute Eigen- 
werteingrenzung ermöglicht. Nach einigen allgemeinen Bemerkungen und Überlegungen über 
die Handhabung und den Anwendungsbereich des Satzes wird seine Brauchbarkeit und Über- 
legenheit gegenüber bekannten Verfahren an einigen Beispielen vorgeführt. 


I. Eine neue Eigenwertdarstellung bei Integralgleichungen 
1. Die Aufgabe und Bezeichnungen 
Allen Untersuchungen wird zunächst die homogene Integralgleichung 
va-ilEadylade Ne. nee 


®) K. Hohenemser, Die Methoden zur angenäherten Lösung von Eigenwertproblemen in der 
Elastokinetik. Sammlung: Ergebn. d. Math. u. ihrer Grenzgebiete, Bd. I, Heft 4. Berlin: J. Springer 1932, 89 S. 
»)L. Collatz, Genäherte Berechnung von Eigenwerten. Z. angew. Math. Mech., Bd. 19 (1939), (ZB) 


8.224—249 u. 297318, (Im folgenden kurz zitiert als: Collatz, Eigenwertbericht.) 


*) L. Collatz, Eigenwertprobleme und ihre numerische Behandlung. Leipzig: Akad. Verlagsges. 
1945, 3388. (Kurz: Collatz, Eigenwertberechnung.) 

°»)L.Collatz, Schrittweise Näherungen bei Integralgleichungen und Eigenwertschranken. Math. Z. 
Bd. 46 (1940), S. 692—708, insbes. 8. 706—708. 

°%) L. Collatz, Einschließungssatz für die Eigenwerte von Integralgleichungen. Math. Z., Bd. 47 
(1941), 8. 395—398. 


”) Dieser Satz wird in der Literatur oft nach Weinstein benannt, man vgl. aber Collatz, Eigen- 
wertberechnung, 8. 208. 


Nun ®) Erdelyi, A. Eigenfrequenzen inhomogener Saiten. Z. angew. Math. Mech., Bd. 18 (1938), S. 177 
8 ; 
»), A. Temple, W. A. Bickl ey, OnRayleighs Principle, London 1933, 8. 21—23. 
... M) R.Iglish, Bemerkungen zu einigen von Herrn Oollatz angegebenen Eigenwertabschätzungen 
bei linearen Integralgleichungen. Math. Ann., Bd. 118 (1941), S. 263—275. 

‚“ı) Diese von Collatz gegebene Erweiterung der Templeschen unteren Schranke ist allerdings 
auch im Buche: OnRayleighs Principle (vgl. Fußnote °) auf S, 86—87 unter den Exercises angedeutet. 


A, vr Rellich, Störungstheorie der Spektralzerlegung. IV. Mitteilung. Math. Ann., Bd. 117 (1939), 


ü | ıb n reellen symmetr che: 

sche Integrabili ät und mittlere Stetigkeit er ERIC 
| Die Eigenverte h von GI. (1) sind reell ‚und können nach Ge DER are 
r % An We Ne 4% u FR ee SAL <A <o a <sh< % ” " = er Aue Su 
| geordnet und beziffert a (Abweichend vom Üblichen ist ER größte negative Eigenv 
mit A, bezeichnet.) Jeder Eigenwert wird seiner Vielfachheit entsprechend oft gezählt. 


Sa _ zugehörigen reellen rn: Y (%) werden als SENGESn Eh eHE und auf Ale Täng 12 
Re RN: normiert angenommen. BUN: 


2. le ze und Ra RR a aus der Theorie der. Bey. 
? Integralgleichungen BE 2 NONE 
Obwohl die Ne der Gl. (1) überall als bekannt vorausgesetzt wird 13), ‚sollen einige iX 
für nr anschließenden ee grundlegenden Sätze in Ger BeRoalee Fassung wiederholt = 
werden. N 


A. Die Lösung der zu Gl.(l) gehörigen inhomogenen Integer 


gleichung 
Fa)=y@)—A/K(ms)y(s)ds.. RN VOTEN He 


ist bei beliebiger reeller in @ stückweise stetiger Störungsfunktion fx)*) (kurz: für jede 
s-Funktion f(x)) eindeutig möglich, wenn der Parameter A nicht Eigenwert ist: A+A,. Bei 
4= 4, ist für die (zunächst nicht.eindeutige) Lösbarkeit notwendig en hinreichend, daR F ) 
auch zu allen zu A; gehörigen Eigenfunktionen orthogonal ist | 


E: ah . Sr) <yı(z)y)de—=0 oder kürzer fo) L <yi® (DDR NE 
- (<Y(a)> mit der spitzen Klammer vertritt hier und weiterhin immer alle möglichen Linear- 


kombinationen der zu A,— A; gehörigen Eigenfunktionen.) 
Mit der Hilfe der Resolventen 


Ras) ee nie) a ya Eee (5) a 


Bild die Lösung der Gl. (3) in der Form | 
yva=fÄa+ifRaszsa)fidds......- EN Sr 


gefunden; sofern Gl. (4) a ist, darf dabei auch A= 4; werden und man erhält dann die 


mit der Zusatzforderung | 
aaO re ee (7) 


auch für den Parameterwert A = A, eindeutig bestimmte Lösung der Gl. (3). Durch die Ortho- 
5 n. gonalitätsbedingung Gl. (4) wird in der Gl. (6) der Einfluß der Eigenwerte A, = A; völlig be- 
seitigt. 

M Den GI. (3) und (6) entnimmt man, daß sich bei festgehaltenem reellen Parameter A + 4, 
a Z: jede s-Funktion f(x) umkehrbar eindeutig durch eine andere s-Funktion y(x) in der Form der 
E Gl. (3) darstellen läßt; y(x) =0 liefert dabei f(x) = 0. _Durchläuft also y(x) den Bereich 
z 

} 


‚der auf die Länge 1 normierbaren s-Funktionen, erhält man mit f(x) aus Gl. (3) wieder alle Br 
Funktionen dieses Bereiches. Ist aber in Gl. (3) A= 4, so ergibt diese Darstellung nur die E 
_ Gesamtheit der zu (y;(x)) orthogonalen Funktionen: f(x) | <y; (x)), dabei genügt es schon "nl 
[man vgl. Gl. (6), (7)], wenn y(«) dem Bereich ae a y(x) | <y;(x)) entnommen 7 
wird; y(x) = <y;(x)) führt hier immer auf f(«) E 
-=B- Die Festlegung der er A, bzw. ihrer Reziprokwerte 1/), = %, 
(der charakteristischen Zahlen) durch ein Extremalprinzip nachCourent 1): 


22) Man vgl. z.B.: R. Courant, D. Hilbert, Methoden der mathematischen Physik I. Berlin: 
Julius Springer. 2. Aufl. 1931, 496 8. (kurz: Courant-Hilber t) oder G. Ha mel, Integralgleichungen, 
Einführung in Lehre und Gebrauch. Berlin: Julius Springer 1937, 166 S. oder G. Wiard a, Integralglei- 
chungen. Leipzig und Berlin: B. G. Teubner 1930, 183 8. 

14) Allgemeiner können hier und im folgenden an Stelle der stück weise stetigen Funktionen (vgl. Cou- 
rant-Hilbert S. 39) auch alle quadratisch summablen (d.h. im Sinne von Lebesgue integrierbaren) 
Funktionen in Betracht gezogen werden. Natürlich sind dann auch alle Integrale und Gleichheiten im Sinne 
von Lebesgue zu verstehen usf. 

15) Die Äquivalenz -ist im Sinnevon Hur witz benutzt. Vgl. z.B. Hamel, Integralgleichungen N 
4.2.0. 8.64 u. ff. 

16) Man vgl. z.B. Courant-Hilbert, 8.112, 113 a.a. O. 
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Es sei v,(2), vg(R), Ir ‚ %{x) oder zusammengefaßt {v}; ein System von 
s-Funktionen und ; | Be. 
M({v};) = Max [| K.(&, s) f(@) f(s) ded m({v})= Min [[K (m, s)f(a)f(s)dads  (SauK 

7 1 07 ; ; 0 1 1,0; 2 


das Maximum oder die obere Grenze bzw. das Minimum oder die untere Grenze der quac tische, 2 n 
Integralform [[ K(a, s) f(x) f(s) dx ds, wenn alle s-Funktionen f(x) konkurrieren, die im 
Grundgebiet auf die Länge 1 normiert und zum System {v}; orthogonal sind: Beh dc 
[f?(a) dad=1 und J[f(a)u(e)da=0 bzw. kürzer Lie: - (9b). 
ürr=1,2,..,5 zn ir. 
(Im Falle j = 0 fallen natürlich alle Orthogonalitätsbedingungen weg, das System {v},istleer.) 
Sucht man jetzt in Abhängigkeit von den Hilfsfunktionen z;(z) den Kleinstwert von , 
M ({v};): Min M ({v};) und ist dieser positiv, so gilt 
v 


Min M ({o}) = Min Max I Kto, )fa)f)drds=1/h4 3-0, Yon 
v v tıly ; 


344 ı.- Lehmann, Beiträge zur Lösung 
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analog dazu liefert der Größtwert von m ({v};): Max m ({v};), sofern er noch negativ ist: 
v 
Max m ({v},) = en SEK fa)fi)deds=14; j=dlh... .. .(10b). 
v vuTL UF) 


Die betrachteten Minimum-Maximum-Werte werden für v, = Y1, da = Ya» ..., y= 4 und 
f= Yj+ı angenommen; die Maximum-Minimum-Werte entsprechend bei %, = Yu, % = Y-ı, 
ee %—= Yı-j und = Y- i e i . 

Die beiden, positive und negative Eigenwerte unterscheidenden Eigenwertdarstellungen 
lassen sich formal zusammenfassen: 


Min Max (> d re 
N Koss a ily.. -..- Ad. 
Er ;>0 net 

es gehört bzw. zu ie: m Me 


Die Minimum-Maximum-Werte (Maximum-Minimum-Werte) können für 5=(,1,... auf- 
gesucht werden, solange die Integralform dabei noch positive (negative) Werte ergibt; das soll 
durch Min Max (> 0) (bzw. Max Min (<0)) angedeutet sein. (Im Falle 4; = 0, also betj= 1 
oder 0, fallen wegen {v}, = leer alle Nebenbedingungen entsprechend G]. (9b) weg, so daß ein 
einfaches Maximum- bzw. Minimum-Problem übrig bleibt.) Die zu den a,;-Orthogonalitäts- 
bedingungen / |_ {v}, ; gehörigen Extremalfunktionen sind bei den Bezeichnungen nach Gl. (11) 
fo) y(a). 

3. Ausnutzung der Resolventen zu einer neuen Eigenwert- 
darstellung. Zunächst eine Vereinbarung: Aus der Folge der Eigenwerte des Kernes 
K(z, s), die nach Gl]. (2) geordnet sein sollen, wird im folgenden immer wieder ein bestimmter 
Eigenwert A, herausgegrilfen; gehört der Eigenwert dabei zu einer Gruppe gleichgroßer (ist er 
also mehrfach), so wird in dieser Gruppe immer der mit größtem Index ins Auge gefaßt und 
unter A; verstanden. 

Weiter soll t; ein reeller Parameterwert zwischen A; und dem nächstgrößeren Eigenwert 
Ai+ı sein 

A<t <Arı En nn er (12), 
Ist einmal A;+1 der kleinste oder A, der größte Eigenwert des betrachteten Kernes, so fällt in 


. Gl. (12) die untere bzw. obere Grenze für t; weg — formal wird in diesen Fällen einfach A; > — 


bzw. A441 + w gesetzt. 

Die zum Parameter A = t; gehörige Resolvente R(z, 5; t;) des Kernes K (a, s) kann nach 
Gl. (5) wieder als Kern mit den Eigenwerten ,—t; = A437 —b >0 für » >ioderrj>0, 
<0 für» <ioder j SO und den zugehörigen Eigenfunktionen 4,(x) = Yi+j (2) aufgefaßt 
werden. Diese Eigenwerte gestatten mit Gl. (11) die Darstellung 


MinM >0) SD Dr 
Max min (<o) Rio 554) Fa) fie) da ds=1/ (44-1) Dre 13). 
v Ita, 3 an 


Die s-Funktionen f(x) müssen dabei nach (9a) auf die Länge 1 normiert sein. Die Eigenwerte 
werden bei f(x) = yi15(x) angenommen. 


Diese Eigenwertdarstellung ist für die Anwendung wichtig, weil jetzt 1/(A4ı — t,) als ein- 
facher Maximalwert, 1/(A, — &;) als einfacher Minimalwert einer quadratischen Integralform er- 


| Damit Sen es nöglich, das Figenwertspektrun 
s Wertes i; zu untersuchen. LE ER 
Die in GI. (13) auftretende unbekannte und A SeöTende Resolvente Tan sich leich % 
beseitigen. Nach den Betrachtungen im Abschnitt 2/A kann bei i; er die Gesamtheit der 
Be: normierbaren s-Funktionen fi) auch in der Form 


Be ya menuas re 


B; erhalten werden; 'y(«) hat wieder a Bereich der normierbaren s-Funktionen zu durchlaufen. % 

277: Beieinem Funktionenpaar f(x); y(x), das Gl. (14) genügt, le nach GI. (3) und (6) die Be- 

E ziehung i | > 
Be. [K: (X, s u N (25; Bu Br: ER (La) - ? 
Fi _ In Gl.(13) führt das A eh | 2 
© © MinMax(>0) x | EEE ar 
Max Min (<0) IK s) en Hedy (s) ds] y( a an 
Bee ya, EIN Er... 

r 3% Die Normierungs- und ie für y( x) folgen aus Gl. (9a, b):: 

Er, : u) —tu[K(a,s ALS 1 und y\ Wa >... 2 (16b,C). 

5 Als y(2) und für das System {v}a, —als0 0, +. ae — ist wieder die Gesamtheit der s-Funktionen B-: 
in Betracht zu ziehen. Es. mußte beachtet werden, daß die Bedingung EAN (a; nach GI. (Ybyrses = 
mit dem Ansatz Gl. (14) zunächst Bi 

(2) [y(@)— 6 [K (8, s) y(s)ds] da= [y® )Ter (2) —4 [KR (@; s) ve(s s)ds]de= Su a 0. er. ; 
für =1,2,...,0; | Br 
d.h: {O}a; ergibt; da aber nach Abschnitt 2/A %, mit v, gleichfalls den gesamten Bereich Me 


‘der a. Fimktigtien darstellt, kann dafür Gl. (16c) geschrieben werden. 


Die Extremalfunktionen yz(x) der Aufgabe Gl. (16) sind Eigenfunktionen des Kernes 
Ki, s), bei der Normierung Gl. (16b) wird y5 (2) = Ai; yırz (a)Aiyz; —t;). 

Die Eigenwertdarstellung nach (Gl. 16) soll in eine noch zweckmäßigere Gestalt gebracht 
werden. Bisher war in Gl. (16) die Gesamtheit der normierbaren s-Funktionen y(x) zur Kon- 
kurrenz zugelassen, die erhaltenen Minimum-Maximum- und Maximum-Minimum-Werte werden Be: 
aber nur verwendet, solange sie bzw. positiv oder negativ ausfallen. Von jetzt ab sollen im 


7 Minimum-Maximum- (im Maximum-Minimum-) Problem’von vornherein nur solche s-Funktionen Bi 

2 konkurrenzfähig sein, die in Gl. (16a) eine . positive (negative) Integralform ergeben; in ‚ 

E- diesem Sinne wird Gl. (16) weiterhin immer verstanden. Dabei bleiben die Extremwerte der 

4 Aufgabe offensichtlich ungeändert. [Denn die interessierenden Extremwerte werden bei be- ER 

stimmten Extremalfunktionen (nämlich den Eigenfunktionen cy;+;) angenommen und bleiben Ca 
erhalten, wenn man dem zugelassenen Funktionenbereich Zwangsbedingungen (Beschränkun- ie 


Br gen) auferlegt, denen die ursprünglichen Extremalfunktionen sowieso genügen !?).] 
Außerdem soll noch die Normierungsbedingung Gl. (16b) wegbleiben; an Stelle der qua- 
dratischen Integralform in Gl. (16a) muß dann der Quotient zweier Formen treten 1%): 


[IK @, 9 Ivo) —uS Kia, s) yis)ds] yi)dads _ 10 Ze 
[va —u[K(® s) y(s)ds)?dx TE: RR 

Der hierbei definierte Ausdruck 
un Ma Ze E@ a yiaddstade ei 


S[Ewyd—t[EKas)yidslyis)daads " ° 


soll wegen seiner großen Bedeutung Te mplescher Quotient genannt werden (kurz: T-Quo- 
tient). Temple?) erhielt diesen Ausdruck in etwas anderer Gestalt als Ergebnis einer Ab- 
sehätzung als untere Schranke für den kleinsten Eigenwert A, eines positiv definiten Eigen- 
wertproblems, dessen A, einfach sein mußte. Dabei blieb der unmittelbare Zusammenhang 
mit dem Eigenwertspektrum, der hier im Mittelpunkt der gesamten Überlegungen steht und 
„zufällig‘‘ auf diesen bekannten Quotienten führt, völlig unsichtbar. 

Mit dem Quotienten von Gl. (17) wird aus dem Extremalproblem Gl. (16) 


MinMax(>0)\ 1 en 152. 
MaxMin(<OQ)/)T,y)—u Ati )=9 1 2, ER IE 
v Be 


17) Man macht sich das leicht, an Hand der Gl.(8) bis (10) klar, oder man vgl. Courant-Hilbert 
8. 201 und 353 u. folg. 
18) Vol. z.B. auch Courant-Hilbert 8.346, 


“ durchgeführte Division durch den Zählerausdruck keine Schwierigkeiten bereitete. Die € tho- 
4. 
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Konkurrenzfähig sind im Minimum-Maximum- (Maximum-Minimum-) Problem alle s-Fu 
nen, für die der Quotient 1/(7, (y, 4) — 1) bzw. positiv (negativ) ist, so daß auch die mn 
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Eee Lehmann, 


gonalitätsbedingungen sind die gleichen wie bei Gl. (16), die Normierungsforderung © (16b) . 

ist weggefallen. j R 
Beachtet man in 1/(7, (y, t) — ii) einmal nur Funktionen y für die der Quotient positiv 

ist (oder: negativ ist), so entspricht im gleichen Funktionenbereich einem ee 


Max RUE An Min 2 
Min "Wert von Rear immer einem yay"Wert von T,yh)—u- & x 
Dieser Satz kann offensichtlich in Gl. (19) angewendet werden: Y 
Min Max (> 0) 1 2, 1 ". EN are 
Max Min (< 0 N \—L; : MaxMi 0 un —4)=9 —b.... 
a AR el rg Er 
v vida, 


Vergleich der Nenner und Weglassen der dabei belanglosen Konstanten ; führt auf die end- 
gültige Eigenwertdarstellung 


Max Min (>;)] | N Hr 
Min Max (<6) |  WFAH =, — 2...) ...... (20). 
f ® via, , 


Hier sind alle s-Funktionen zu berücksichtigen, für die 1/7,(9,t;) =0 gebildet werden kann 
(in ‚Gl. (19) war ja 1/(7,(,t) —t) = 0 auch ausgeschlossen); im Maximum-Minimum- und | 
im Minimum-Maximum-Problem werden beziehungsweise nur die Funktionen verwendet, j 
für die-[7, (6) —t] >0 und [7,(y, 4) —t;] <O ausfällt, wie das auch in Gl. (20) 
zum Ausdruck gebracht ist. Bei dieser letzten Aufteilung des zugelassenen Funktionen- 
bereiches wird immer noch die Gesamtheit der s-Funktionen erfaßt, die 1/7',(y, t;) = O ergeben. 
s-Funktionen, für die 7,(y, 6) = t; ist, können nicht vorkommen, da sonst nach Gl. (17) 
Ivy) —t [ K(®, s) y(s) ds?4&= 0 und daher t, im Gegensatz zur Voraussetzung ein Eigen- 
wert sein müßte. 

Die Eigenwertdarstellung nach Gl. (20) liefert den gesuchten Ansatzpunkt für die beab- 
sichtigte numerische Eigenwertberechnung. Bevor jedoch diese praktische Aufgabe in Angriff 
genommen werden kann, müssen noch einige Grenzfragen erledigt werden. 


4.Das Verhalten der Eigenwertdarstellung Gl. (20), wenn &; 
Eigenwert ist. Der in der Darstellung Gl. (20) auftretende Parameter t; war nach der 
Definitionsgl. (12) bisher nicht Eigenwert. Gelegentlich wird aber auch der Fall £;— 4; (oder 
ganz analog: t;— A;;ı) auftreten und berücksichtigt werden müssen. 

Da in der Gl. (20) die rechte Seite von t£; völlig unabhängig ist, kann damit auch für die 
linke Seite der Grenzwert 1; — A; oder — A;+1 definiert werden; dabei ist aber die Konvergenz 
für die Gesamtheit der konkurrenzfähigen y(x) nicht gleichmäßig. Setzt man von vornherein 
t; = A; fest (die Betrachtungen für t;—= A;;ı verlaufen ebenso), so zeigt sich, daß 7,(%, A,) 
für y= <y;(x)) eine isolierte Lückenstelle ergibt, der als Grenzwert 2;—4; der Wert A; zu- 
geschrieben werden kann; löst man die Aufgabe Gl. (20) mit #5; = A; , so fehlen in der Gesamt- 
heit der Extremwerte alle Eigenwerte A, = A,, es ist 


Max Min (>4,)] a Mr ART 
Min Max (<A) a Miu — (+1, — (ut di.. 
ES, a gr 

; k bij=—(;+k, k>0 


Hierbei bedeutet o, die Vielfachheit von A;. Die Funktionen y = <y;> erweisen sich als nicht 
konkurrenzfähig, weil 7, (<y;), A;) unbestimmt wird, 

Um diese Behauptungen zu beweisen, wird der Quotient 7,(y, A;) zuerst einmal für die 
Funktionengruppe <y;) und dann für die Gesamtheit der übrigen s-Funktionen betrachtet, Der 
erste Fall ergibt, gleich als Grenzwert 


. (21% 


tl 
Im Kylie N 
ni ı (<yi) t) Fi 1/7, —t/R A; a re (22); 


1%) Der Beweis für diese Beziehung kann auch direkt, ohne den Umweg über G1.(13) bis (16) geführt 
werden. Dabei müssen allerdings die Min-Max- und Max-Min-Eigenschaften der Form von Grund auf neu 
entwickelt werden, so daß dieser Weg etwas umständlicher wäre. 


4 bare nen y=+<yir 


| 3% dabei jöt PRO) infolge y+ AN eine nicht identisch es wegen H 


(RS i via)= Vu mit Kar en Yin a Y% @ 


Sm) yu@)da= ya) [y(e) —<y (al da= [yin) yula)az— [yıs) ae Nie. 


formel Gl. (20) auf die Forderung y | {y;>. Für den Funktionenbereich y j {y;) fehlen daher 
in der Darstellung nach Gl]. (20) die Eigenwerte A, = 4;; alle übrigen Extremwerte bleiben un- 


berührt, da auch die eingeengte Funktionenklasse die zugehörigen Extremalfunktionen (die 


Y + <yı») enthält. Weiter hat A; auch für den Parameterwert ti; seine Sonderstellung verloren, 
es darf £; = A; gesetzt werden, womit sich dann über Gl. (24) sofort Gl. (21) ergibt. Und daraus 


folgt schließlich, daß 7,(y,A;) tatsächlich keinen Punkt in der Umgebung von A; annimmt. 


Der Wertevorrat von T,(y, ti) zeigt folgende charakteristischen Merkmale: 


Er RT iR RAN | a: 
beih<u<a,, DE - JIIIITITEIL (III = Wertevorrat) 
Ai Aitı 
bei ;—1t LIT TITITEITITT 


dio, (A als Grenzwert einer Lücke) 


' Die Möglichkeit ,; =A;:ı erledigt sich entsprechend. 


Mit Gl. (21)/(22) ist das Verhalten von G]. (20) auch für die Grenzlagen von t; völlig ge- 
klärt: Werden bei Gl. (20) die Funktionen <y;> von der Konkurrenz ausgeschlossen, so darf 
im Maximum-Minimum-Teilproblem ohne Änderung der Extremwerte t;—= }; angenommen 

werden, im anderen Teilproblem fallen die Extremwerte },— 4; aus; beim Grenzübergang 
t; > A, können immer alle Eigenwerte gefunden werden. Für den Fall t; = ;;ı und y # (yi+1) 
bleibt umgekehrt Gar Minimum-Maximum-Teilproblem unberührt, im zweiten fehlen alle 


? »= Er Ai+ı- 


5.Der Templesche. Quotient und das Verfahren der ehrite 
weisen Näherungen?°)*)°).. Bei der Lösung der Integralgleichung (1) nach dem 
Iterationsverfahren geht man von irgendeiner s-Funktion y(z) = y'”(x) aus und bestimmt 
nach 
y9ı (x IK (0) yr m )de= [EM (a,s) yle)ds ı n=L2.... (05) 


iterativ eine ee ya), y’(a), yıılz),.... K®%(z, s) ist dabei der n-fach iterierte 


Kern. Damit-werden Schwarzsche Konstanten 


aly)= fyrmi(z RER REM ER ER (26) 


gebildet, sie sind von m unabhängig. Sobald @a,+1(4) + 0 vorausgesetzt wird, lassen sich auch 


die sog. Schwarzschen Quotienten 


N 25 DE (27). 


An+1!Y) 
anschreiben. Der Templesche Quotient Gl. (18) nimmt mit diesen Bezeichnungen die Ge- 
stalt an: 


HN—iHaM)_ „amt _, mM —ry] og 
rare Pr Kaly)—i Bir 


B: E | 
(wobei Yul 2) € <yıla)) jede zum Eigenwert A, — 4; gehörige Eigenfunktion sein kann) zu $: 
<yil a)> orthogonale s-Funktion. An Hand der GI. (18) stellt man nunmehr fest, daß = nCrENeH 


Tılyr )= T,(y.%) HR ee Le Er IE: EE en 


- Fir den Beweis von Gl. (21) brauchen also, da die Funktionen y= Zi) nicht Ka Kuren Ar; 
we fähig waren, nach Gl. (24) von vornherein nur alle zu (y;) orthogonalen Funktionen yL(y) 
N in Betracht gezogen zu werden. Wiederholt man jetzt die Überlegungen, die von Gl. (13) u 
Gl. (20) geführt haben, mit dem Funktionenbereich f | (y;), so erkennt man bereits beim An- 
satz Gl. (13) (s. auch die Bemerkung unter Gl. (7)), daß damit der Einfluß der Eigenwerte 
 %,= 4, völlig ausgeschaltet ist; die Bedingung f | (y;) führt nach Gl. (14) und (6) in der Schluß- . 


Br? nd 
d.h, falls a,(y) #0, sofern also 1/7,(y,t) nicht unabhängig von t verschwindet (mit a,(y % 


[+ 2 Br v Lupe ee a ra .-4 hg Rn ’ a 
=] AB ne ur . j Zee du v % 
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(Die Existenz der verwendeten Quotienten ist natüı Te 
dabei deutlich die Abhängigkeit 7,(y, 1) = fit) bei fest gewähltem y. Die S: 
Ungleichung liefert dazu noch 


arly) = LfyO la) yD(a) dar < [Iy”(a)Paz- [Ir®(al’de=aty)aty) - + 


En 


ist nach Gl. (29) auch a,(y) = 0!), bleibt 


u. — ala = usly) ul) —ueNM)S0 !,r. +... BO 
a yi Sobald also 7,{y,t) überhaupt von t ab- 
hängt, liegt nach Gl. (28) bei t= u,(Y) RE 
ein einfacher Pol und es ist immer mit 
t<u,(y) oder >u,(y) bzw. T,(y, 0) Zus >t 
oder <u,;<t, und umgekehrt folgt aus 
T,(y,t) >t oder <t bzw. m >t oder 
us <t. Tıy,t)=fil) stellt eine gleich- 
seitige Hyperbel mit den Asymptoten 
t= 4, T,=yu, und um 1350 gegen die 
t-Achse geneigter Hauptachse dar. Man 
vgl. Bild 1. Das gezeigte Verhalten des 
T-Quotienten wird noch oft benötigt. 

In die Templeschen Quotienten 
T,(y, t) können auch iterierte Funktionen 
y")(z) eingesetzt werden. Der Satz nach 
Gl. (20) bzw. (21) bleibt dann gültig, wenn 
die Gesamtheit der »-ma] iterierten Funk- 

Bild1. Verlauf der Funktion T,(y,t) = fl). tionen zur Konkurrenz herangezogen wird; 

diese Funktionen sind alle sogar stetig und 

es finden sich darunter wieder die Eigenfunktionen von X (z, s) und damit die alten Extremal- 
funktionen zur Gl. (20) und (21). Mit Gl. (26) bekommt man 


an (y) = | ym yMda= a, (y) ei N und | 


1 n=0,1,2... fürn >1 
an) = [ [Ka 5) ym (a) y ()drds= a, (y) 
und 


nl) a Nm warn (32). 
I m ra, lyM) al) ty) re 


Nachdem gezeigt ist, daß der Satz nach Gl. (20) für 72 +1 (%, £) gültig bleibt, taucht die Frage 
auf, ob er auch für den naheliegenden Quotienten 


T,(y®, = a, (y) —ta,(y®) a Aon+1ıly) —taanı2(Y) 


Bi = Tania) ce RC 
a,(ym) —taz(y) Aon+e(Yy) —tAgn+s(Y) +2 (4 8) ) 


benutzt werden kann. 


6. Der Templesche Quotient T,(y,d). 

Zunächst läßt sich sofort feststellen, daß die Gl. (20) für T,(%, ti) nicht mehr gelten kann, 
sobald der Kern K(z, s) der zugrundeliegenden Integralgleichung positive und negative Eigen- 
werte besitzt. Man überlegt das einfach beim Sonderfall t,—=0. Nach dem Entwick- 
lungssatz gilt für die einmal iterierte Filnktion yd (x) — | Ki«, s) y(s) ds= & c, y(a)/A, mit 
g= f y(x) y,(x) dx; für T,(y,0) nach Gl. (26) und (33) also ud 


7,40), Welse (34). 
v y 


Da positive und negative Eigenwerte vorkommen, erkennt man hierbei, daß der Quotient 
T,(y,0) in Abhängigkeit von y(x) (d.h. also den Cy), Im Gegensatz zu den mit Gl. (20) be- 
wiesenen und ausgenutzten Eigenschaften von T,(y,0) (vgl. Bemerkungen über den Werte- 
vorrat auf S. 347), jeden Wert zwischen A,...0.. .4, annehmen kann. 

Sobald nur positive Eigenwerte vorhanden sind, der Kern der vorgelegten Integralgleichung 
also positiv-definit ist, entfällt der soeben aufgezeigte Widerspruch und es läßt sich zeigen, daß 


dann Gl. (20) auch mit 7,(%,t) formuliert werden darf. Mit einer Hilfsfunktion Y(x) kann 


T,(y, t) jetzt immer in die Form T,(y, t) gebracht werden. Das zu y(x) gehörige %(x) wird 
in der Gestalt 


v=23ne [yawddafh = ylally.. 2... (35) 


2 Puhl on TR da es Eigener Ay wegen des. Ban “ 
<a ie sein. müssen, außerdem nach Rieß-Fische.r sicher summabel 14) u a 
[Kia s) y(«) y(s) deds = 2 c3/A, die. Summe der De der in Gl. 65) es! Ent “ BR 


wicklungskoeffizienten konvergiert. ran BG 3. 
Der Entwicklungssatz und Gl. (35) ergeben die ER E Ba 5 


am = BL (2, 5 N ()drds= 2 5, Y (x) deu) .: oo) 


Bund: 


& 2,(y) = 2a; - -Ktas) 72) die) dade- ud). = EER- a. 
Daraus folgt 


ne ta,(y) N) _ 
U en rg 
aly)—tazly) ya) 
7 Die a x) können dabei als s-Funktionen gelten, sie lassen sich nämlich durch N bbreehen der 
 Reihendarstellung Gl. (35) immer so durch sogar stetige Funktionen annähern, daß diein 7, (, t) 
mit diesen s-Approximationsfunktionen erhaltenen Integrale (vgl. (36)) von den formal mit 
Y (x) gebildeten beliebig wenig abweichen. Der mit y(x) nach Gl. (35) gegebene Funktionen- 
bereich enthält alle Eigenfunktionen des Kernes K(x, s), so daß T,(y,t,) = T,(y,t;) beim 
Durchlaufen der Gesamtheit der s-Funktionen Y( x) wie in Gl. (20) die Aurj als Extremwerte a 
darstellt. ee 
II. Ausnutzung der gezeigten Eigenwertdarstellungen für die numerische er. 
. Eigenwertberechnung 2 
7. Bin EinschlieBungssatz für die Eigenwerte des Kernes K(25): 
Kann um den Eigenwert A; ein Bereich abgegrenzt werden, der nur diesen einen Punkt- 
. eigenwert enthält — kennt man also eine untere Schranke };+1 für den nach A; nächstgrößeren a 
Eigenwert A;+ı und eine obere Schranke Li, für Ki-op dem nächsten Eigenwert unterhalb P 
des o;fachen 4 (A, <Lio, <A <Lbirı <A), so gestattet jede s-Funktion y(z), die Be: 
der Bedingung 


Pe Eee de RE DS 
B: genügt, die Eingrenzung i Se 
4 | 1719, Li—o,) > 4,> Ti l;+ı) AT (38b). I E 
2 = Ist einmal A; der kleinste‘ (oder größte) Eigenwert der Aufgabe, darf formal ı 4, 4 
ä km; >— x (kıı>+o) und damit u, ——o(li4ı> +&) gewählt werden, es B. 
32 ‚strebt dann T7,(y,t)—us(y). Wird eine der Ungleichungen in Gl. (38a) nicht erfüllt, so fällt I: 
| auch in Gl. (38b) die zugehörige Schranke aus. Die Eingrenzung mit Gl. (38b) wird im all- x 
Ä gemeinen um so besser, je kleiner L;_,, und je größer l;+1 gewählt werden kann, für eine Eigen- R 


e- funktion y; ist T,(y; t) = us(y;) = A und von t ganz unabhängig. 
k Zum Beweis setzt man 4; = 4< bi, = h<i4ı = A und <huı=b<in. 


Wegen Gl. (38a) ist nun u,(y) >t; bzw. u,(y) <t; und damit (vgl. Bild 1 und bei Gl. (30)) 
auch T,(y %) >us >t bzw. T,(y,t;) <u; <t. y(x) ist also im Minimum-Maximum- 
Problem von Gl. (20), bzw., wenn dort t, statt t; gesetzt wird, auch im Maximum-Minimum- 
- Problem konkurrenzfähig. Bedenkt man nun, daß das Minimum von T,(y,t) >t Arıı 
ist (nach Gl. (20) ohne Orthogonalitätsbildung), erhält man den ersten Teil von Satz (38b). 
Der zweite Teilfolgt analog; man vgl. auch die Betrachtungen über den Wertevorrat von 7, (y, t;) 
- am Schluß von Abschnitt 4. Ist eine der Bedingungen von Gl. (38a) nicht erfüllt, so entfällt 
lediglich der darauf bezügliche Teilbeweis. Die angegebene Abhängigkeit der Schranken von 
L;_., und I;;ı verfolgt man am einfachsten am Bild 1. 
Nach den Beweisen im 5. und 6. Abschnitt darf in Gl. (38) 7T,(y,t) immer durch 
Tantı(y, t) ersetzt werden; an Stelle von u,(y) muß dann u2n+2(Y) treten. Im Falle eines 
positiv-definiten Kernes kann sogar allgemein 7„(y,t) zusammen mit n+1(y) für n>1) 
Verwendung finden. 
Über das Verhalten von 7',(%, t) (oder Tan+ı(%, t)) bei n— oo orientiert man sich bequem 
mit Hilfe des Entwicklungssatzes. Da nach Gl. (25) bereits ie quellenmäßig dargestellt 
wird, gilt für n>1 


ym (2) = = I | ie 92 len mit = [ y(®) ©) Yol2) dw lodr 


v 


20) Vgl. z.B. Courant-Hilbert S. 93/94. 
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Ist neben A,, dem absolut genommenen kleinsten Eigenwert, bei ı 
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koeffizient e, nicht verschwindet, — A, kein Eigenwert, wird 
1l=1für,=) - 
i „| <1 für alle anderen 4, | 
und es konvergiert damit | 
my) Kyla)) fürn... ren 
((y:(x)) ist eine zu A, gehörige, allgemein unnormierte Eigenfunktion.) 


Daraus ergibt sich wieder unabhängig von be N | 
lim I, (Y; ) nn A er re (41). 
n oo 


Den bisher ausgeschlossenen Fall, in dem neben A, auch — 4, Eigenwert ist, erkennt man 
daran, daß die y®(x) für gerades und ungerades n nach voneinander linear unabhängigen 
Funktionen konvergieren. Auch dann läßt sich an Hand von’Gl. (39) überlegen, wie daraus 
gegen Eigenfunktionen konvergierende Folgen konstruiert werden; wegen der geringen prak- 
tischen Bedeutung erübrigt es sich, darauf genauer einzugehen. 


8.Der Einschließungssatz von Kryloff?). 

Die Eigenwertdarstellung nach Gl.(20) bietet auch einen interessanten und verall- 
gemeinerungsfähigen Zugang zum Einschließungssatz von Kryloff. N 

Die aus Gl. (20) hergeleiteten Betrachtungen über den Wertevorrat von T,(y,4) am | 
Schluß des 4. Abschnitts zeigen anschaulich, daß für jeden beliebigen Parameter 4—=t und 
jede s-Funktion y(z) zwischen £ und T7',(y, t) mindestens ein Eigenwert liegen muß — sofern 
auch 7T,(y,t) selbst dieser Eigenwert sein darf. Ob t einmal Eigenwert ist, bleibt belanglos. 
Damit liegt der Gedanke nahe, den Parameterwert t so zu bestimmen, daß die durch ?= ty 
und dem mit einer beliebig gewählten s-Funktion y(x) und iy gebildeten Quotienten T, (%, tx) 
gelieferten Eigenwertschranken möglichst eng werden: 


[T,(yt)—t®=Min d.h. 2[7, yy—1 [da T,(ytY)at—1j=0.... (2). 


Der erste Faktor der angeschriebenen notwendigen Bedingungsgleichung für das Minimum 
kann nur dann einmal verschwinden, wenn 7,(%,£) ein Eigenwert und daher von ? unabhängig 
wird. (Man vgl. die Bemerkungen unter Gl. (20), danach folgt aus T,(y,t) =t, daß t Eigen- 
wert, y Eigenfunktion ist.) Im allgemeinen muß zur Bestimmung des Minimums d T, (y, /dt=1 
herangezogen werden. Mit Gl. (28) folgt daraus für die gesuchten Parameterwerte ty und 


T,(9 tm): ; 

im, = Hey) + Varly) laıly)—azly)] und T,ly, ta, ) =usly)F Varly) ı)—asly)]) (43). 
Da zwischen im. und 7,(y, im, ,o) immer mindestens ein Eigenwert liegen muß, ist damit der 
Satz von Kryloff bewiesen 


NY arly) ln) —usly)) <A <usly) +VRely) ad) ae] - - - (44). 


Die Betrachtung von Tan+ı(Y, t), bzw. bei positiv-definitem Kern auch von Tan+2(%, f) 
zeigt, daß in Gl. (44) an Stelle der Schwarzschen Quotienten 4,(%), #2(y) noch allgemeiner 
Kan+ı(y) und uen+2(y) bzw. sogar n(y) und tn+1(Y) benutzt werden können. 

L,—ı AT Iszı 
„T, (%, Le-i) 


As+ı 


dam 


My) u (Y) 
Ma — Vf (in —ue) | — A — Hy + Yrz (u) —us) 


Bild 2. Schematische Darstellung der Ergebnisse nach Gl. (38), (44) und folgende. 


Sobald mit Hilfe einer s-Funktion y(x) für einen i-ten Ei { j tig: 
f genwert gleichzeitig Schranken 
aus Gl. (38) und (44) erhalten werden, dürfen die Schranken und Shhraplene 


K Ynlayıa] 

4; mad AR) l#ı (Y) —4:(Y)] und d= T,(y, Lu) — Tı(y: l+ı)- ae (45a) 
AR für Gl. (44) für Gl. (38) 

miteinander verglichen werden. Sie sollen. von Null verschieden sein, sonst geben beide Ver- 

fahren ja den gleichen Eigenwert exakt an. Da die mit Gl. (38) bestimmten re (sofern 

sie nicht selbst den Eigenwert darstellen) nach den dort angegebenen Bemerkungen um so 


9 N 
ON 


e kleiner L;_,, und je grö rä ‚kann, gi 
E nik er daß die Schranken von Gl. (88) a nd a ann günstiger si 
Er yı off schen en wenn Itı > dar, bzw. dns, <ta, bleibt. D.h. es ist 


i 
we a . .. a I NT, 
R > 


ae wenn Di, und l;ıı außerhalb der in Gl. (44) bestimmten Schranken Tiegen. 
 - Kommt es also auf enge Eingrenzung eines Eigenwertes (relativ zum Abstand von de 
Nachbareigenwerten) an, und sind halbwegs brauchbare L; ., und bıı bekannt, muß dahe: 
die Schranke nach Gl. (38) dem Kryloffschen Satz weit überlegen sein. Man vgl. hierzu 
die Bilder 1 und 2 sowie die Beispiele A (mit Bild 3) und B im 11. 3 beehuiie ” 


9. Bemerkungen zum Bereit Bone ien der DE ER EEN SAtzei 
Fi Die Anwendung der bewiesenen Sätze ist keineswegs auf Integralgleichungen von Typ 5: 
der Gl. (1) beschränkt; ohne auf irgendwelche Vollständigkeit zu achten, werden Bach einige 
für den Praktiker wichtige Anwendungsgebiete zusammengestellt. 
Zunächst braucht das Integral in Gl. (1) nicht wie bisher nur im Riemann schen ; 
Sinne verstanden zu werden. Die Schmidt-Hilbertsche az und: "42 
damit alle hier aufgetretenen Sätze bleiben auch für 


-A[Ka,)yli)dal)) -. - .- EEE Le Inte) 
7 gültig, wenn das Integral im N Stieltjes benutzt wird und die „Dichte“ o(«) 
g monoton mit & wächst ?!). Bei der Anwendung der gezeigten Beziehungen sind dann einfach 
alle vorkommenden Integrale als Stieltjes integrale mit dem Dichtedifferential d o(«) 
aufzufassen; quadratische Integrabilität z. B. ist immer in bezug Aue die Dichte a) zu Vver- 
stehen usf. 

Weiter soll nur erwähnt werden, daß sich auch Systeme von N a Wrat: 
einer einzigen vom Typ der Gl. (1) oder (46) zusammenfassen lassen und damit den hier auf- 
gestellten Sätzen zugänglich sind. Ebenso bieten Integralgleichungen mit allgemeineren her- 
miteschen Kernen keine Schwierigkeiten ®). Bei polaren Integralgleichungen, deren Kern 
positiv-definit, die Dichtefunktion o(x) nicht nur steigend ist, übernimmt der T-Quotient 
T,(y, it) die Rolle von. 7,(y, 0). ko s 

Wichtiger ist der Zusammenhang der Integral- und Differentialgleichungen. Das Ver- 
bindungsglied ist die Greensche Funktion *) im gewöhnlichen oder erweiterten Sinne. 

Es sei ein EN vom Typ 


(y)+Aoy=0, Ux{y) = 0 für =1,2,. ee Pe 


“ vorgelegt. Dabei ist Z ein Er homogener selbstadjungierter Ebenen oder partieller 
Differentialoperator, o(z) eine im Grundgebiet positive s-Funktion. U,(y) sind zugehörige 
lineare homogene Randbedingungen. Die Existenz einer symmetrischen quadratisch-integrablen A 
Greenschen Funktion mittlerer Stetigkeit @(x, s) (im gewöhnlichen oder, wenn A,=0 Eigen- N 
wert ist, im erweiterten Sinne) wird vorausgesetzt, so daß die reziproken Beziehungen ‚Va 


| L%W)=—u(2); U,()=0 entspricht v()= [E(&,)us)ds ...... (48) er. 
bestehen. .Als v(z) kann jede s-Funktion vorgegeben werden, lediglich bei der Greenschen Bi; 


Be, Funktion im erweiterten Sinne muß u(x) zusätzlich zu allen zum Eigenwert A, = 0 gehörigen 
I Eigenfunktionen der Differentialgl. (47) orthogonal sein: w_| <yp». Eine mit G(z, s) quellen- 


ig mäßig dargestellte Funktion v(x) bezeichnet man als ‚‚Vergleichsfunktion“; sie ist nach Gl. (48) | 3 
& dadurch gekennzeichnet, daß sie die Randbedingungen U,(v) =0 von Gl. (47) erfüllt und die 
- _ — Bildung Z(v) = (s-Funktion) zuläßt; falls in Gl. (47) A, = 0 Eigenwert ist, muß auch v | <yo> 
e- sein. 
er Dem Randwertproblem Gl. (47) ist nun die FABENEIUHR 
—=1/[6(x, s) y(s) o()ds=A [@(x, s) aje(d) ae RT 


21)K.neser, A., Die Integralgleichungen und ihre A, in der mathematischen Physik. 292 8. 
Braunschweig: Vieweg 1922. 2. Aufl. 

Die ‚‚belastete "Integralgleichung“ bei Kneser hat noch nicht die allgemeine ‚Gestalt von Gl. (46), 
jedoch bemerkt man sofort, daß die von ihm auf $. 119 als wesentlich aufgezählten Eigenschaften des be- 
lasteten Integrals auch beim Stiel tjes-Integral mit monoton wachsenden Dichtefunktionen erfüllt sind. 

22) Man vgl. hierzu z. B. die in Fußnote 2 genannte Arbeit von Hohenemser. 

23) Zur Theorie allgemeiner hermitescher Kerne vgl.z.B.Hoheisel, A. Integralgleichungen. 136 S. 
Sammilg. Göschen Nr. 1099. Berlin u. Leipzig: de Gruyter 1936. 

24) Vol. hierzuCollatz, Eigenwertberechnung, Courant-Hilbertund besonders Kamke,E., 
Differentialgleichungen, Lösungsmethoden und Lösungen I, 6668. Leipzig: Akad. Verlagsges. 1944; hier 

iA; finden sich auch weitere Literaturangaben. 
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- äquivalent. Während diese Aufgabe sonst durch Symmetrisierung behandelt wi: 
sie hier der Gl. (46) untergeordnet. Damit ist die Anwendung der für Integrz -hungen be 
wiesenen Sätze auch auf die Randwertprobleme der Gl. (47) möglich. Auf Einzelheiten soll 
hier und’ auch im folgenden nicht eingegangen werden, es kann auf die Literatur ?°) und e' t a MR 
die Beispiele im 11. Abschnitt dieser Arbeit verwiesen werden. Re as 
- Ebenso wie die Randwertaufgabe Gl. (47) lassen sich die nach Collatz Eingliedauf- 
gaben genannten eindimensionalen Randwertprobleme 
# L(y) HA 1 (gym)® = L(y) HAN (y) V.y)=0 .... (0a) 
9(x) >0 und n-mal differentiierbar im Grundgebiet ug 


(der Index n deutet bei dieser Gleichung einmal auf die Zahl der erforderlichen Differentiationen) 
auf eine Integralgleichung vom Typ der G]. (46) zurückführen, wenn unter den Randbedingungen 
Y)y=y()= = y"d()=gb)=ylb)= = yedb)=0 .. . (506) 
vorkommt ?®), 
Für die von Kamke?”) theoretisch, von Collatz) vom Standpunkt numerischer 
Eigenwertberechnung betrachteten allgemeineren Randwertprobleme L(y) +4N(y)=0 (und 
Randbedingungen), wo N (y) ein allgemeinerer homogener definiter selbstadjungierter Differen- 
tialoperator ist, gelingt die Zurückführung auf Integralgleichungen nicht allgemein ®®). Die An- 
wendung der abgeleiteten Sätze würde dadurch aber keineswegs problematisch: Der Beweis 
der Beziehung Gl. (20) läßt sich auch mit den schon bei Collatz: Eigenwertberechnung, ver- 
wendeten grundsätzlichen Hilfsmitteln durchführen; am einfachsten, wenn alle benutzten Funk- 
tionen als entwickelbar vorausgesetzt werden (für die Anwendung ist aber wichtig, daß man auf 
diese Annahme verzichten kann). 

Zum Abschluß sei nur noch bemerkt, daß mit Hilfe des Entwicklungssatzes (Spektral- 
zerlegung) die Anwendbarkeit der aufgezeigten Verfahren für das diskrete Eigenwertspektrum 
bei selbstadjungierten Randwertproblemen (Operatorengleichungen, Matrizen) auch dann ge- 
sichert werden kann, wenn das Problem keiner Integralgleichung der erwähnten Formen 
äquivalent ist (im Falle z.B. auch ein Streckenspektrum oder ähnliches vorkommt). Es wird 
beispielsweise an die Schrödinger- Gleichung der Wellenmechanik gedacht. Sofern sie 
rein reell selbstadjungiert ist, wird sie formal ganz ähnlich wie z.B. die Aufgabe Gl. (47) 
behandelt. 


u) 


u en Tec 


10. Bemerkungen zur Anwendung und Handhabung des Ein- 
schließungssatzes. 


Der Einschließungssatz nach Gl. (38) gestattet unter Ausnutzung einfacher bekannter 
Näherungen für die Eigenfunktionen und leicht erhältlicher grober Schranken für die Nachbar- 
eigenwerte sehr oft eine außerordentlich enge Eingrenzung des zwischenliegenden Eigenwertes; 
der nötige Rechenaufwand bleibt minimal. Die Näherungen für die Eigenfunktionen können 
dabei z. B. aus Ersatzsystemen mit endlich vielen Freiheitsgraden (Differenzenverfahren) und 
Interpolationsformeln, aus benachbarten exakt lösbaren Problemen, aus asymptotischen Eigen- 
funktionsdarstellungen oder direkt aus physikalischen oder anderen Erwägungen entnommen 
sein; natürlich kann man diese Methoden auch alle geeignet koppeln (indem man das benach- 
barte System vielleicht so wählt, daß die Eigenwerte dieses Systems wenigstens asymptotisch 
mit den Eigenwerten des vorgelegten übereinstimmen) ®°), Es ist für die Anwendung wesent- 
lich, daß der Einschließungssatz Gl. (38) über die Abweichungen der Näherung von der Eigen- 
funktion mittelt, und daher auch stark gestörte (im Sinne der Störungsrechnung) Aufgaben zur 
Bestimmung der erforderlichen Eigenfunktionsnäherungen herangezogen werden können. Wenn 
beispielsweise die Rellichschen !?2) Konvergenzkriterien der Störungsrechnung längst nicht 
mehr erfüllt sind, können die Eigenfunktionen des ungestörten Problems in (G1.38) noch aus- 
gezeichnete Resultate ergeben. 

Die in Gl. (38) notwendigen Schranken für die Nachbareigenwerte dürfen im allgemeinen 
sehr grob sein; sie werden etwa Vergleichssätzen entnommen oder mit Hilfe des Tem ple- 

25) Insbes. Collatz, Eigenwertberechnung, a. a. O. 

?6) Die bei Collatz, Eigenwertberechnung gestellte Bedingung 3n <2m — 1, wo 2 m die Ordnung 


an auch in der von Collatz gewählten Beweisführung wegbleiben. Siehe in Diplomarbeit des 
?) Kamke, E., Über die definiten selbstadjungierten Eigenw on bei öhnli i 

Bu an elglaichn en II u, III, Math. Z. Bd. 46 (1940), 8.231.286 a 
) Während bisher nur Sonderfälle auf Integralgleichungen zurüc eführt werden konnt B. bei 
Picone, E. Atti della Reale accademia dei Lincei, Rendiconti, Olanse die scienze eiche ee ? 
naturali Ne S ER, alles - Bukebe a Verf. inzwischen als Integrodiff.-Gl. allgemein gelöst 
en Überblick über diese praktischen i y i y 
ER @ s ragen gibt besonders Collatz, Eigenwertberechnung 


; werteingrenzungen ermöglichen)®) 2). 


_ Iteration bei der Ermittelung höherer Eigenwerte schädlich wirken. Ist die Eigenwertaufgabe. 
als Differentialgleichung der Form Gl. (47) gegeben, läßt sich leicht eine Iterationsvorschrift 
_ finden, die zu einer Eigenfunktion führt, welche zum nächsten Nachbareigenwert einer 


LW+roey=-I-Lytbeyta—mWey-iytroy Uy)=0 ..... on 
Zum Operator L gehören jetzt offensichtlich die Eigenwerte A = (A, — t;), deren absolut ge- 


‚eine Funktionenfolge y9, y®, ..., y, ... ‚die nach oben angeführten Überlegungen entsprechend 


“ wendung ergibt meist die Nebenbedingung Gl]. (38a); insbesondere lassen sich die geforderten 


Mit T,(y,t) bleibt die Rechenarbeit meist am Kleinsten. Nee bei der ae 
absolut genommen tiefsten Eigenwerte empfiehlt sich gelegentlich der Einbau eines Iter 
schrittes und Verwendung von T7,(y,t). Bei positiv-definitem Kern kann auch 7,(y, t) 
nutzt werden, gegenüber T,(y,t) bedeutet das Einschaltung eines halben Iterabionscehrifte 
(aus T,(y,t) entsteht nach Gr (32) und (37) T,(y,t) und T,(y, t), wenn y(«) = y®(«) bzw. Sea 
durch y)(z) oder nach Gl. (35) ersetzt wird. vgl. von Gl. (35) und die Darstellung von ah 
nach Gl. (39) geben die Berechtigung zu den verwendeten Bezeichnungen). Die Iteration ver-r 
stärkt, relativ zu den anderen Eigenfunktionen, die zum (absolut genomnien) kleinsten Eigen- 
wert gehörigen Komponenten der Ansatzfunktion (s. Gl. (39)). Daher kann diese gewöhnliche. a 


vorgegebenen Zahl t, gehört. Man zerlegt dazu Gl. (47) in 


nommen kleinster durch den t; nächstbenachbarten ursprünglichen Eigenwert Ay bestimmt ist; 
die zugehörigen Eigenfunktionen sind (yv(x)), die ja gerade berechnet werden sollen. Während 
das Iterationsverfahren entsprechend Gl. (25) bei der Bandwertaufgabe Gl. (47) mit Gl. (48) 
die Gestalt 

L(y) = — o yar—D( DE Ol een ER 
annimmt und auf eine zum absolut genommen kleinsten Eigenwert gehörige Eipenlunktion 
führt, gibt 


L(g®) = — eye (a) = L(y) + yo yM, vg) —0 ER. SO 


Gl. (40) bis auf einen Faktor (yx (2), ae irgendeiner Kombination der zu Ay gehörigen Eigen- 
funktionen, zustrebt. (Man beachte, daß die Zerlegung in Gl. (51) dem Übergang von K (z, s) 
zu R(z,.s;t,) reziprok ist!) 

Die Bedeufung des Einschließungssatzes nach Gl. (38) liegt in den sehr oft erstaunlich 
engen Schranken, die man damit bei geringer Rechenarbeit erhält; häufig können sogar für alle 
Eigenwerte gleichzeitig Eingrenzungen gewonnen werden. Man berücksichtige hierzu und all- 
gemein zu diesem Abschnitt die anschließend durchgerechneten Beispiele. Grenzen für die An- 


Schranken für die Nachbareigenwerte dann nicht finden, wenn der betrachtete Eigenwert selbst 
mehrfach ist ohne daß dies der Rechner erkennt. (Überhaupt wird die Methode dort ungünstiger 
arbeiten, wo die einzelnen Eigenwerte eng beieinander liegen.) Wie in einer weiteren Arbeit 
gezeigt wird, können die angeführten Schwierigkeiten zum größten Teil durch Verwendung 
mehrgliedriger Ansatzfunktionen beseitigt werden. 


11. Beis Diele 
Die folgenden einfachen Beispiele sollen einen Überblick über die Handhabung, den er- BR: > = 
forderlichen Rechenaufwand und die Leistungsfähigkeit des Einschließungssatzes nach Gl. (38) 
geben. So 
A. Als erstes Beispiel wird die Sturmsche Eigenwertaufgabe betrachtet. E 
Sie ist durch Transformation ?°) immer auf die Liouvillesche Form zu. bringen 
y"+(en?g (2) +2) y =0=L ( y)+Ay; Rkı,Y(0) + kız y’ (0) = 0 (53) 
kayl)tia.yt)=0 y 
Das Fundamentalsystem asymptotischer Lösungen (d.h. für A—oo) der Diiferential- 
gleichung (53) ist sin Ya, cos JA x; asymptotische Eigenwert- und Eigenfunktionsdarstellungen 
können durch Eintragen des Ansatzes y= (, sin YAx-+ O,cos YA x in die zugehörigen Rand- 
bedingungen erhalten werden. Die aeymptokiichen Eigenfunkuionen dürfen als Näherungs- 
funktionen in Gl. (38) verwendet werden und liefern meist sehr gut brauchbare Schranken für 
alle Eigenwerte. Genauer wird das für den Fall der ersten. Randwertaufgabe durchgeführt. Es 
sei jetzt q(x) der Einfachheit halber beschränkt und 
y()=y()=0 und daher „,-sinina=yW(a) i=lh2%....... (SA). 
(= = asymptotisch) 
23 


N nl Vergleichsfunktionen (s. be N 
gegeben als einmal iterierte Funktionen angesehen werden ; unter Beachtur re 
5 oder auch (B2a)ist | N 
in yo) [O3 y"(Del)ds oder Ziym)=— ey"; Urt)  K 


v 


und speziell 


Vu [yin]” + eng (a) yD=—yM bei yr =) =0. Re. 
0 Geht man von y(® aus, wird y{® durch einfache Differentiationen gefunden und damit ergeben 
E sich hier die Schwarzschen Konstanten Te 
I Ei ! 1 44 ! Du ‚ 
Bin = | (rdo= ZT —2etm dit 2er Bi; = | yP ymds=-—-— erh | 
En BIP Ä ö 7, (56). 
SEAN v % r h . 
Br A f (ydrdı= z mit A=[q() sin’izzdz und 2 Bf (2) sintinzdz 
2 n : { J 
% 1%; A; und B; sind also im wesentlichen Fourier- Koeffizienten der Funktionen g(z) und g°(2). 
ER  Sielassen sich in der Praxis zumeist durch Aufzeichnen von g(z), 9°(2) und Anwendung des 
Br; harmonischen Analysators bequem für genügend viele ö ausrechnen. Mit den Schwarz- 
Ir Ks . schen Konstanten hat man alle Bausteine für den Kryloffschen Einschließungs- 
30 . satz nach Gl. (44) [mit (27)], es wird: 
®. RT NER - (57) | 
RE oder bei A, auch: ), <m(1—284,) : 
Hierbei wurde gleich verwendet, daß diese Aufgabe höchstens endlich viele negative 


Eigenwerte haben kann und daher die Numerierung ohne Rücksicht auf das Vorzeichen der 
Eigenwerte beim kleinsten mit 1 beginnend erfolgen darf. Die Anwendung der verschiedenen 
Verfahren bleibt davon unberührt, alle Ergebnisse gelten auch bei negativen Eigenwerten. 
Ebenso wurde sofort benutzt, daß formal A, > — sound somit in Gl. (38) zweckmäßig L,— — 
gesetzt wird. Das gibt A, < u,, die zweite Zeile in Gl. (57). 

Wegen der beschränkten Störung jer”q(x)| < je; M ist nach bekannten. Vergleichs- 
sätzen?)®) |A;(e)— n2| < je| M und daher in Gl. (57), zumindest solange sich die einzelnen 
Schranken nicht überdecken, /, = 4, zu setzen. 

Mit Benutzung der meist sehr groben Schranken von Gl. (57) gestattet der Satz nach 
Gl. (38) eine in weitem e-Bereich sehr enge Eingrenzung sämtlicher Eigenwerte der Aufgabe. 
Durch Einsetzen der Zahlen aus Gl. (57) als Z;_, und l;+, in Gl. (38) bekommt man 


or 4&(B, — 4?) 

24 +1— 22 (Ay — 4) — 21e] Y Bi — uR 
4 &8(B,— A?) 

PR VERS TERID ER ER IP aa: Fir) 


Die Schranken in Gl. (58) sind wenigstens solange gesichert, alssich die verwendeten Kryloff- 
schen Schranken nicht überdecken; die Bedingung Gl. (38a) ist dann erfüllt. Natürlich liegt es 
hier nahe, für}, L,—— oo (also A, <4s), Lı = u, usf. auszunutzen; man bekommt durch 
rekursive Anwendung von Gl. (38) genauere Schranken als mit Gl. (58), es ist 


\ g2 ee, 
L=n? (2A ae (Bi E 
”—2e 4 —L_ı[m 

He Bi An) 


em (RP 224 — a) mit „=n"(®—2e A —2le Yß; | 


m? 2e4— 


h (98). 
2 (? —22.4+ 


) mit ,=n?(1 -28A,) | 
(59). 


[Die zur Berechnung unterer Schranken benötigte Ausgangsschranke l; (k beliebig) braucht 
selbstverständlich nicht unbedingt aus Gl. (57) entnommen zu werden.] 


An einem Zahlenbeispiel wird die Brauchbarkeit der gewonnenen Schranken deutlich. 
Es sei 4 (2) = «, dann ist A; = 1/4, 4 B,= 1/3 — 1/(2 27); es sei darauf hingewiesen, daß 
dies rechnerisch einfachste Beispiel nicht ausnahmsweise günstige Ergebnisse zeigt! In Bild 3 
wurden die für die ersten drei Eigenwerte nach den verschiedenen Methoden gefundenen Schran- 
ken in Abhängigkeit von e aufgetragen. Es gehört zu Gl. (57): 12; OZ, zu Gl. (58): 2; O,, 1; O 
I,; Ol und GI. (59): A; O,, 1; O, 15; ON (mit k= 4) a IT 


I 


8 = 4 ei 2 % er 
1 V aus Gl.(57). 


- 2 =70 
Unter el tung von lg. er REN | 
iese Werte folgen mit sin ir « als Vergleichsfunktionen aus 
gehören dazu’ die Schranken 17; DT NER N ÜBER ‘ 
end sei bemerkt, daß die Mittelwerte 8 und 6,5 auch von der Störungsrech 
\äheru gen, 1. Ordnung erhalten werden, jedoch erlaubt das Rellichsche Krit au 
bis e= 0,3 die Konvergenz von A; (e) zusichern EEE DE NE 
d den Fehler abzuschätzen.. on N a RER 
B. Es wirddie (auch bei Collatz, EBENE 
Eigenwertbericht diskutierte) Aufgabe \ \ Ir N ET 
VERLTON | REN SEEN 


dem Templeschen 


| 
b N HN BER 


ohne Transformation mit dem Ansatz 

 y®=sinia® behandelt. Hierbei läßt sich 
leicht ein Iterationsschritt einbauen, aus 
WOr=—U+D Hy M-0 
ergibt sich h 


2 ur OR (1-++-&)sinin x es 
N: U, ie 43708 i , az 


EAN NL REN Eee 

VER er a A N: 
a re ge bei Dedhd1 
an u (61). 


oo 


Mit a=— | 


_ ee 


7 
schen Quotienten den Wert 


B.: 3072 
N, Eee, IR 5% 62. 
E; de Ita 3 UTC, 11, INN 5 AR ( 2) 
er 7 42m? = 
e: 5 . Dahier ein positiv-definites Problem vorliegt, 


können beide 7-Quotienten 7, und T, ver- 
wendet werden. Für die ersten drei Eigen- 
werte ergeben sich die Schranken (in der fol- 
genden Tabelle geben die Zahlen der Reihe | = 
nach: obere, untere Schranke, Mittelwert der Bild 3. Schranken für die ersten Eigenwerte der Aufgabe Gl. (53)/(54) “ a 
verstärkten Werte mit Fehlerschranke in %): gr er 


x a Methode 3 | Schranken für A, Schranken für A, Schranken für A; 


= + Yan) *) | (6,58) 4,04 6,36 |31,5 21,5 26,3 70,6. 47,8 Bu 
Ha 4 Vs (U,—Hs) (6,50) 6,16 +3,1% | 31,0 18,4 +18% 68,9 44,8 +18% 
2 er) T(y, t) : 6,58 6,24 6,55 27,5 25,4 26,2 63,3 52,6 20 
{ T,(y,t) 6,56 6,54 +0,15% 126,9 23,1 + 3% 61,7 49,7 48% 
| {6) 
»##),Demplescher & 9,8 4,9 7,4 +33% | 39,5 19,7 29,6 33% 88,8 44,4 67 + 33% 
Bar » | v1 65684 5% |985 25,8 272 45% | 641 581 614 5% 
= *) Die in Klammer gesetzten Werte sind u, bzw. u; allein. u 
7, ' er) Für L, wurde u, bzw. u, benutzt, als /, ‘der aus ‘dem Vergleichsproblem y’+i2y=90, y() = 
5 y (l) = 0 entnommene.Wert I, = 8n?. Weiter wurde nach Gl. (59) rekursiv verfahren. 
Br, ***) Templescher Einschließungssatz°) ®)5)%) a) mit sinınz, 
Er b) mit V3,5— sin[(in In(1—x/3,5))/ln (5/7)] (s Collatz, Eigenwertbericht) als Vergleichsfunktion. 
z | 23% 
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356 Buchholtz, Reihenentwicklung einer zweireihigen Determinante a. 


An diesem Beispiel ist bemerkenswert, daß der grobe Ansatz für die Eigen 
rung — die nicht einmal asymptotisch mit den Eigenfunktionen übereinstimmt — m Gl. 
zu durchaus brauchbaren Eigenwertschranken führt. Weiter gibt die Aufgabe ein Beisf 
den Einfluß eines halben Iterationsschrittes; dabei wird die erste Eigenfuı mponente 
verstärkt (vgl. Abschnitt 10), so daß sich die 7-Quotienten in der Rich ung T 17 1 FA 
ändern. Auch bei diesem ungünstigen Beispiel zeigt sich die Überlegenheit von Gl. (38) gegen 
über den vorher bekannten Verfahren. Eine weitere Verbesserung würde sich bei Benutzung 
. der oben im Temple schen Einschließungssatz unter b) eingetragenen Vergleichsfunktion 

ergeben [vgl. die Ergebnisse a) und b)], allerdings erfordern dann die notwendigen Quadraturen 
größere Mühe. EA In 

C. Es wird die (von CollatzÖ) behandelte) Integralgleichung 


y(ia)=ıflae-siyi)ds ..... BE (63) 


auf den absolut genommen kleinsten Eigenwert untersucht. Mit der Ausgangsfunklion „P—=1 
hat Collatz die Iterationsfunktionen y bis y® und die zugehörigen Schwarzschen 
Quotienten u; bis 4, errechnet. Er kann damit in diesem einfachen Fall den absolut genommen 
kleinsten Eigenwert A, (hier positiv) eingrenzen: 


a 


unter Verwendung von Up...;4 A1= 23,865 + 0,5% 
Up. Me A 2,87804 + 0,017%. 


Dabei wird auch eine untere Schranke für den Betrag des nach A, absolut genommen nächst- 
größeren Eigenwertes A, (es kommt A, oder A,in Betracht) gebraucht. Mit der Trefftzschen | 
Schranke ®°) und u, 4,tindet Collatz |},| > 4,657. Nun kann auch Gl. (38) zur Anwendung | 
kommen. Aus |A,| 24,657 folgt A, <— 4,657 <A, <-+ 4,657 <A, (da ja A, der A, dem 
Betrage nach nächstbenachbarte Eigenwert sein sollte), die erforderlichen Schwarzschen 
Quotienten werden von oben übernommen. Dabei stellen sich gegenüber der von Collatz 
angegebenen Eingrenzungsmethode wesentliche Vorteile heraus: einfachste Handhabung, 
bereits mit zwei Schwarzschen Quotienten können Eigenwertschranken ermittelt werden, 
die Eingrenzung ist besser. (Man sieht, daß die Schranken nach GI. (38) auch bei wesentlich 
komplizierteren Kernen als |® — s| leicht berechnet werden können, während das beim Col- 
latzschen Verfahren fast ausgeschlossen ist.) 
Die Resultate lauten: 

mit u, tg vaRor Lan 

‚mit u, 4 also T,: A,)= 2,8749 + 0,16%, 

mit uu,tg also T,: A, = 2,878337 + 0,012%. 


(Da der Kern nicht positiv-definit ist, darf nur 7Tg„+1 benutzt werden.) 


Bisher wurde absichtlich nicht erkannt, daß hier zufällig A, der größte Eigenwert ist, sonst 
könnte man mit w,, ug sogar A, = 2,878462 + 0,001% behaupten. 


°) Trefftz, E., Fehlerabschätzung. bei Berechnung der Eigenwerte. Math. Ann. (108)/4 (1933), 
S. 595 —604. i e 3% z ER 
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Besondere Reihenentwicklungen für eine häufig vorkommende 
zweireihige Determinante. mit Zylinderfunktionen 
und ihre Nullstellen 


Von Herbert Buchholz in Darmstadt‘ 


Im ersten Teeil der Arbeit werden für die Ausdrücke (1) und (2) einige anscheinend bisher unbekannte 
Reihenentwicklungen aufgestellt. Sie schreiten entweder nach den Kugelfunktionen &% mit }= 0,1, 2... 
oder nach halbzahligen Besselschen Funktionen wachsender Ordnungszahl fort, oder es sind Neu- 
manmsche Reihen zweiter Art. Im zweiten Teil der Arbeit werden diese Entwicklungen dazu benutzt 
um die Nullstellen von (2) in bezug auf x oder die von (1) in bezug auf v angenähert darzustellen. Dabei 
wird eine Nullstelle von (2) berücksichtigt, die bisher anscheinend übersehen worden ist. 

In the first part of the paper, some series expansions for the expressions (1) cnd (2) are stated; they seem 
io be unknown hitherto. They proceed either with the Le gendre functions O wüh A=0,1,2,... 
or with the Bessel functions of the first kind and of the order n + us 
Neumann series of the second kind. In the second part, these ewpansions are used to represent ap- 


proximatively on the one hand the zeroes of (2) relating to x and on the other the zer: i 
088 1 2 h 
One zero of (2) is taken into account that has evidently been overlooked till now, Ve 


1/2, n being an integer, or they are 


Be z Dans la deuxitme Be 
ai ne approzimativement les 26ros % 2: par rapport. & z ou les zeros de (1) pa = 
ER: zero de (2) y est pris en consideration, & Tequel on semble, Jusqu’d resent, ne pas avoir fait a 


N 'B epBoü yacru padors BEIPAsKeHHS da) u (2) pasnaramrea # pansı,. KOTOpsBIe no 
= a op 6suım, HOBH/IHMOMY, HEHSBECTHEI. ITU PAnEI COCTABICHH HAH 13 dymkuni Jlemar 
0 Mei=0, 1,2... num us hyueumä Beccenm mopanka n + 12, HIU-ke OHM ABIAIOTCH pan 
A: Hebsama BTOpOro mHopanka. Bo BTopoä yacrı padorsı aTu PAABI CHOAEBYIOTCH 
IIPHÖNMKeHHOTO ONpeneneHns KopHeH BEIpaskeHnun (2) OTHOCHTEIBHO z H BhipazkeHnust 
 OTHOCHTEIBHO v. Ilpı 5ToMm yunTsBaeren onuH KopeHs ae a BR 10 
Hop ue Öpuı, Ber, oÖHAaPYy:KeH. & 


I: I | EL; Einleitung 
Vom „Verfasser wurde bereits an anderer Stelle [1] nachgewiesen, daß für ae auf der 
linken Seite von Gl. (1) auftretende zweireihige Determinante die darin auf der rechten 


J (R ( url 1 + In (Rjr) Ir Es A 
ET RE Be I En“ vn 
RAS T.eN url 12) —= e®».J (2 wg) dd... we ” e 
q \ — In (Ar) ° E 
enge RN RE Integraldarstellung mit der Erklärungegl. (18) für wobesteht.. 
i = RT DER HE DT NEN FE TREE da) 


v und 2 können in (1) SRH u oder komplexe Zahlen sein, während R N, r mit ER ir 
zwei reelle Parameter bedeuten sollen. | 
Differenziert man die Gl. (1) zunächst nach R und dann ein zweites Mal nach r, so gelangt. 
man nach einigen naheliegenden Umformungen des Integrals. rechter Hand, die wegen ihres rein 
formalen aakiegs hie» nicht interessieren, zu der Beziehung (2).. Löst man hierin Cof © in 
Ra Isa] ek Een 
2 r2 
3 var) Ya | 0 deren) 
Fi Das (R/r) i EN 
die Summe der beiden Exponentialfunktionen exp (+ ®) auf, so läßt sich jedes der entstehen- 
be -, den drei einzelnen Integrale vermöge (1) durch eine Determinante darstellen, und es entsteht 
© auf diese Weise die Gl. % In dieser Form kann sie auch bewiesen werden, indem man in - 


—&0|-a0. (2) 


y2 
Rr.-22 


ya 
Br lzrz) vn ) 
der Determinante (2) J und Y, mit Hilfe eines bekannten Rekursionsformelpaares durch 
die Summe oder die Differenz von J,+1ı und Y,+ ı ersetzt, die Determinante zerlegt und in 
denjenigen beiden von den vier entstehenden Determinanten, deren beide Spalten Zylinder- = 
funktionen mit verschiedenen Zeigern enthalten, wiederum J,_ı und J,+ı durch J, und J, 5 Be: 
ausdrückt. Be 
Kombinationen von Zylinderfunktionen nach Art von (1) oder 2) spielen eine große Rolle 
beim Studium der elektromagnetischen Wellen in kreiszylindrischen Hohlleitern mit einem 
axialen Mittelleiter, und es ist für diesen Zweck erwünscht, für sie sowohl über bequem bereehen- 
bare Reihenentwicklungen zu verfügen als auch die Lage ihrer Nullstellen etwa in bezug auf 2 
© zu kennen. Dabei ist » in der Regel gleich einer der ersten nicht negativen ganzen Zahlen, der 
Bi Parameter r = 0,05...0,5: R und z reell. Ein weiteres häufiges Vorkommen von (1) findet 
bei der Untersuchung der Wirbelstromerscheinungen in kreiszylindrischen Schirmleitern statt. 
r und R sind in diesem Fall in der Regel nur wenig voneinander verschieden, während dann z 
eine komplexe Zahl von der Form |z|-y7 ist. 
Be Im folgenden Abschnitt 2 werden wir uns zunächst mit einigen bisher nur wenig oder noch 
= gar nicht bekannten Reihenentwicklungen für das gegenüber (1) etwas verallgemeinerte Integral 
; (4) ner Im Abschnitt 3 wird in aller Kürze die Frage nach den Nullstellen des Deter- 
+ In (Rjr) +1n (R/r) 
| X (; Rn) = era a We) n=n. | ae 
E 2 ws)? . (2° ws)? 
3 re ö — In (R/r) 0 
minantenausdrucks (2) behandelt werden. Dabei wird das Hauptgewicht auf der Herleitung 
einer Potenzreihenentwicklung für eine Nullstelle von (2) liegen, die anscheinend bisher über-, 
sehen worden ist, obwohl ihr in den Anwendungen die größere Bedeutung zukommen dürfte, 


} ‚Be A,((Re)', =: Ayzı (Rz, rz) u. — ' Ay_ı(R2,r2) — 


.dd..(M 


i ‘ Betr eyes die Funktion, n. (2: we)/ | 
En konvergente Reihe zu ersetzen und. 


= 1. N ER 
- Um Integral a zu einer 


wurde bereits in der eingangs zitierten Arbeit [1] verfahren, Nach a sung in 
‚ führt man dann zweckmäßig wegen ) EAN 


Re+r—2Rr- or Eee) i Bi ; En 


die Sl Estiheton) e® — t durch. Dabei werden die neuen Integralgrenzen r/R und Rich ir 
dem Übergang von i zu a Variablen a der Gleichung ‚a 


jr) 


u 


gehen die Tofegralgtenzen für u in 0 und : über, und man hat dann ein Were 


‚Integral vor sich. Führt man in der zugehörigen hypergeometrischen Reihe die bekannten vu a 
'Kummerschen Transformationen durch, so gelangt man zu den vier Reihenentwicklungen (I, 


die erstmalig von P. Schafheitlin [2] auf dem Wege der Reihenumformung gewonnen 
worden sind. Sie gelten für alle Werte von z, solange (R? — r?)/r? oder R®< 1 ist. Wir führen 
sie explizite an, weil I von ihnen für die numerische Berechnung ihre besonderen kleinen 
Vorzüge hat. 


a! 
vr 2 BEER 
- 5)" 2 De ee MAR UAR i+1; 21.1) 


Et Be, hi | 
RvrrR—r <£ R: 2 Fa 

-(&)" 2P. RE Se a] )1 Platt, A413 224 2:15) | 
X, (2; R,r)= a) 

zR R2— r?\21 ke; 

Fr FR de. j77 ) “A! Yu 

-H Pa FA @A+ı] 'ıF, (atı4n +1; 2I+2; IR & 

zr R—2r2 R 


Herzen Kara Flatitr, +1; 22+2; ee 


Die in (1) auftretenden hypergeometrischen Funktionen können aber nach dem Vorbild 
der Transformationsformel’ 
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auch noch einer quadratischen Transformation mu + werden. Von den vier weiteren da- 

durch entstehenden Ausdrücken geben wir nur die beiden zu der Beziehung (2) vereinigten 

Formen an. Die darin vorkommenden hypergeometrischen Funktionen lassen sich aber durch 
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(die Kugelfunktion &% mit dem Argument (R?+ r)/(R?— r?) ausdrücken. Aus der zweiten 


Zeile von (2) entsteht dann die für jedes » gültige Reihenentwicklung (3). Sie konvergiert wie 
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die Reihen in (2) absolut und gleichmäßig für jedes beliebige zund alle R >r. Der Kugelfunktion 
in (3) liegt wie in der bekannten Formelsammlung von Magnus und Oberhettinger 
die Definition von Hobson [3] zugrunde. 

Ein besonderer Vorzug der angegebenen Reihenentwicklungen ist darin zu erblicken, daß 
die Größe 2, die in manchen Anwendungen eine komplexe Zahl ist, nicht auch als Argument 
der hypergeometrischen oder der Kugelfunktion vorkommt. Ist v»= n positiv und ganzzahlig, 
so besteht die Kugelfunktion ©, nur aus endlich vielen Gliedern. 


cn (5; R,r)= - (3) 


2.2. Eine nach halbzahligen Besselschen Funktionen fortschreitende 
Reihenentwickung für x , 


Mit Hilfe der durch die Gl. (la, b) definierten beiden Hilfsgrößen ö und n lassen sich w& 
und ® in Gl. (1.4) wie folgt schreiben 
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Mit n als neuer Integrationsvariablen entspricht der unteren Grenze ®= 0 in (1.4) der Wert _ 
n=0.und der oberen Grenze ®= In (R/r) der Wert n=n/2. Demnäch erhält man also 
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Wegen der absoluten und gleichmäßigen Konvergenz von (4) für ein 6 <I1 darf in en 
nach dem Einsetzen von (4) gliedweise integriert werden. Die dabei entstehenden ‚Integrale (5) 
lassen sich ihrerseits auswerten, indem man die Besselsche Funktion in eine unendliche 


Be u, #5 


EEE . j. 0 
- “ ae cs PY: i 5 = Si F 
” Fi 2% IT I a Ar = viry 2 % Bon 
“ ’ B m; ° r Sam ni. .) Q y . 4 EB. a F ® 
R 7 er ne aim Datreih an Det m | 1 
Reihenentwicklung einer zı nee 
» ad > F, NA re fe z_ 


Bi. 


BL (a + a) ne. Tea u 
— Dp+T. (ajaytr+ı2 rn ana part 


durch Gl. (5) angegebenen Formel geführt. Für p = 0 fällt darin die in der geschweiften lamme 
stehende endliche Summe fort. . 

Geht man mit diesem Ergebnis in die erste Entwicklung ein, so resultiert schließlich nach 
einigen rein formalen Umrechnungen die Reihendarstellung (6). In der endlichen Summe von (6), 
die für p = 0 gleichfalls fortzulassen ist, ließen sich noch die hypergeometrischen Funktionen 
als Ableitungen darstellen. 
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Ist in (6)v= m + 1/2 mit m=0,1,2..., so geht die Summe über A nur von 0... m 
Im besonderen folgt aus (6): 
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Auch die Entwicklung (6) ist für die numerische Berechnung von Be 


rechnung der Funktionswerte von Ji_p+172 erfordert dann aber wegen des komplexen Argu- 
mentes z(R — r) einen beträchtlichen Zeitaufwand, es sei denn, daß R— r sehr klein ist. Die 
Entwicklung (6) bietet aber in anderer Hinsicht sehr große Vorteile. Hier erwähnen wir nur 
ihre sehr bequeme Verwendung in solchen Fällen, wo das Verhalten von X für ein R,r > 
bei R—r= const untersucht werden soll, wenn also der kreisringförmige Bereich in einen 
Streifen konstanter Breite entartet. Gerade in diesem Falle ist sie wesentlich vorteilhafter als 
die etwa für = 0 aus (6a) dadurch entstehende asymptotische Entwicklung, daß man links 
für J, und Y,ihre asymptotischen Näherungsausdrücke einsetzt. 


2.3. Zwei Neumannsche Reihen zweiter Art für 4, 

Nach einer von Gegenbauer stammenden Verallgemeinerung des Additionstheorems 
der Zylinderfunktionen gilt bekanntlich auch für komplexe Werte von z die Gl. (1), in der die 
rechts stehende unendliche Reihe absolut und gleichmäßig konvergiert, solange exp (+ ®) < Rir 
ist. Geht man mit dieser Entwicklung in die Gl. (4) ein und integriert gliedweise, so erhält man 
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De Beziehung (2). Das darin auftretende Integral mit dm Gegenbauerschen Polynom 
03? Jäßt sich für beliebige Werte von p auswerten. Wir wollen uns aber hier damit begnügen, 
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‚solange z reell ist. Für komplexe Werte von z bleibt zwar die Konvergenz bestehen. Die Be- 


1 
1 
i 
h 


An 
ee 
Bet 


Reno nr SER ROM, Beikenenbrankrung einer olkeiiigen; Delekmnaal URL 


es nur für Dir 0 zu berechnen. Wegen I'(p):C CH” (Evi 8) — 2: Coj (A®)/A für : und 
1=1,2,3.!.sowie 0) (Cof 8) = 1 führt in diesem Fan die Gl. (2) zu der für alle z und v mit 
Ausnahme von v»—= A gültigen absolut konvergenten Neumannschen Reihe (2a) mit der 
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Definitionsgl. (2b) für die Hilfsgröße «. Für » = A tritt in (2a) für das betreffende Reihenglied 
der durch den Grenzübergang v—A aus (2a) DEE Ausdruck (3), der seinerseits für A= 0 
den Grenzwert «a hat. 
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Macht man in 2.2 (3) für p= 0 von der Fourier-Reihe (4) Gebrauch, so erhält man 
für X’ die Neumann sche Reihe zweiter Art (5), deren Koeffizienten A; sich aus der Gl. (5a) 


"Jı(zR) J,(zr) 


Een . 


Inl2.2-coso)= 2: Sn Jı(&) Ja (x) cos 5) N A (4) 
berechnen. Zu der für A, angegebenen Reihenentwicklung gelangt man durch die erlaubte 
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gliedweise Integration der Gl. 2.2(4). Ist r nur wenig kleiner als R, so ist natürlich die Ent- 
wicklung (5) wesentlich besser konvergent als die Reihe (2a), zumal bei einem kleinen Wert von 
R — r der beiden Reihen gemeinsame Nachteil, der in dem Auftreten des häufig komplexwertigen 
z im Argument der Besselschen Funktion liegt, nicht so schwer ins Gewicht fällt. 


3.1. Die Wurzeln der Gleichung 4,(h’, (R/r:h)') =0 
Im Hinblick auf die Gl. 1.(2) und (3) und auf die Gl. 2.2(6a) läßt sich schreiben: 
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Wir betrachten den Fall kleiner Werte von a. Sieht man dabei von dem Faktor vor dem Summen- 
zeichen ab, so verschwindet die linke Seite von Gl. (2) in erster Näherung, wenn entweder 
1. für alle v» >0 h=»-(1—.a)-!2 ist, weildann die eckige Klammer zu Null wird, 
oder 
2. für alle ah=nm mit m=1,2,3,.., weil dann J,,, (um) = (2/7 m)\?- sin (um) zu 
Null wird. 
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Wir haben demnach in der Funktion f(, t) von Gl. (2), die der für alle gleichmäßig un... 
genten unendlichen Reihe der Gl. 3.1(1) ae eine in x und t analytische Funktion vor 


ER N 
ER N, 


nn AR 


= 1 r ,. Jr bY2rO+ m) _ b „ga p) 
ZEN ar ne 


— 


% 


uns, die für {= 0 in @= 0 verschwindet. Nach einem allgemeinen Satz aus der Theorie der 


impliziten analytischen Funktionen muß demnach die Möglichkeit bestehen, die Wurzel der 
Gleichung f(z,t)=0 in bezug auf & zum mindesten in der unmittelbaren Nachbarschaft von’ 
t=0in eine EN nach ? zu en und dasselbe muß auch von jeder analytischen 
Funkton von x gelten. 
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so gilt nach einem bekannten Satz [4] für denjenigen Wert von Ah, der die Gl. 3.1(1) aus dem 
unter 1. erwähnten Grunde zu Null macht, im Hinblick auf (la) die durch Gl. (4) angegebene 
Entwicklung. 
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Auf Grund von (2) besteht nun für die Koeffizienten C,(x) in der Relation (5) die Deter- 
minantendarstellung (6). 
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Die nun noch erforderliche Berechnung der O, («) aus (6), (9) und (9a) ist lediglich eine 
rein formale Angelegenheit, deren Durchführung hier nicht besprochen zu werden braucht. 


Sie führt schließlich für 2r/h zu der nach geradzahligen Potenzen von (R— r)/2yr R fort- 
schreitenden Potenzreihenentwicklung (10). 


a 


RE 3. Die Entwicklung von 2 z/h nach geradzahligen Potenzen von (R—r) Rynd 
für he. nmla 


In ganz ähnlicher Weise kann man vorgehen, um eine Entwicklung für die zweite Wurzel- 


serie von 3.1(1) nach steigenden Potenzen von (R—r)/2/r R zu erhalten. Wir können uns aus 
diesem Grunde kurz fassen. Nach Einführung der beiden Variablen i{ und x von Gl. (la, b) 
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legt man jetzt zweckmäßig den weiteren Betrachtungen die Funktion f(x, t) von Gl. (2) zu- 
grunde. Die Bedeutung von b, (2) und b, (x) ist damit im vorliegenden Falle unmittelbar klar. 
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Mit dem Symbol C; (x ) von Gl. 3.2(5) gelangt man dann neuerdings zu der Beziehung (3). 
Die Auflösung der hierin auftretenden Integrale verlangt zwar ziemlich umfangreiche Rech- 
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nungen. Sie bieten aber keinerlei besondere Schwierigkeiten, so daß sie hier übergange 
en. Die schließlich entstehende Entwicklung (4) erscheint auch hier als e aA a 
zahligen Potenzen von (R—r)/2 Yr R fortschreitende Reihe. a 
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*" Von MacMahon [ö]sind für h, selbst schon vor langer Zeit ähnliche Formeln auf- 


estellt worden, die jedoch nicht als Potenzreihen nach dem Verhältnis (R—r)/2Yr R, sondern 
An asymptotische Entwicklungen für große Werte von m anzusprechen sind. Gemäß dem Auf- 
bau von (4) ist es aber möglich, beide Entwicklungsarten in einem gewissen Umfang miteinander 
in Übereinstimmung zu bringen, und dieser Umstand erklärt es auch, warum einige der Mac- 
Mahonschen Reihenentwicklungen selbst für den Wert m = 1 noch brauchbar sind, falls 


eben nur gleichzeitig das Verhältnis (R— r)/2Yr R genügend klein ist. 


4. Die Wurzeln von A4,(Rz, rz) =0 in bezug auf den Zeiger v 


Wesentlich schwieriger ist die Berechnung der Nullstellen von 1.(1) oder 1.(2) in bezug 
auf den Zeiger v. Mit einer Aufgabe dieser Art bekommt man es z. B. zu tun, wenn die Fort- 
pflanzungskonstante elektromagnetischer Wellen [6] in einem gleichmäßig gekrümmten un- 
begrenzt langen Hohlleiter rechteckigen Querschnitts bestimmt werden soll, dessen beide Kanten 
von der Länge 2a in der Ebene der Krümmung liegen, während die anderen beiden Kanten von 
- der Länge 2b senkrecht dazu stehen. Die Größe z in unseren bisherigen Gleichungen steht dann 
zur Abkürzung für den Wurzelausdruck (k? — (nr n/2 b)?)!/? mit k als Wellenzahl. Er ist neben 
Rund r mit R—r = 2a als gegeben anzusehen, und es gehört nun gerade zur vollständigen 
Lösung der zugehörigen Randwertaufgabe, in Rücksicht auf die vorgeschriebenen Randwert- 
bedingungen den Zeiger », der zugleich in dem Exponentialfaktor exp (? vg) auftritt und der 


umgekehrt proportional der Winkelgeschwindigkeit der Welle ist, so zu bestimmen, daß für den. 


einen Wellentypus die linke Seite von.Gl. 1.(1) und für den anderen Wellentypus die linke Seite 
der Gl. 1.(2) verschwindet. 

Einen ersten Einblick in das Verhalten der Lösung gewährt der Grenzfall des gerad aus- 
gestreckten Hohlleiters. Für R„= (Rr)!2 — x geht vg in das Produkt y -z aus der Fort- 
pflanzungskonstanten y und der Koordinate z in Richtung der Längsachse des geraden Hohl- 
leiters über und 9 R„ in z selbst. Mit zunehmendem Krümmungsradius R,, des gleichmäßig 
gekrümmten Hohlleiters muß sich also » mehr und mehr dem Wert y- R, nähern und mithin 


4? y?(R—r)? 
gehen und somit endlich bleiben. 


4.1. Eine weitere nach halbzahligen Besselschen Funktionen und Potenzen von ö fort- 
schreitende Reihenentwicklung für X 
Unter diesen Umständen wird es sich empfehlen, die Gl. 2.2(6a) vor der Berechnung der 
Nullstellen in bezug auf » so umzuformen, daß darin » nur im Produkt mit Ö auftritt und danach 
die Berechnung der Wurzel nicht für » selbst, sondern für » ö vorzunehmen. 


Man stellt nun leicht fest, daß für die Entwicklung von (», A) nach wachsenden Potenzen 
von » die Beziehung (1) gilt, in der die Zahlen De) durch die Gl, (2a) oder (2b) definiert sind. 
Die 29 sind darin Polynome p-ten Grades in n. Ihr gesetzmäßiger Aufbau ist bis zu ?o=6 
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in dem Buch von N.E. Nörlund [7] zu finden. Geht man dann mit (1) in die Gl. 2.2 (6a) 


ein, so gelangt man nach einer einfachen Reihenumformung zu der Formel (3) mit der Erklärungs- 
gleichung (3a) für die innere Summe S$,, die für alle Werte von o und 2 »2ö? absolut konvergiert. 


© P) . 
xy (5; B,r)=—n- 4, (Ra rn) 28 mn. 2 e=(R—r)'z.. (3) 
vu x 
1 
(arsasn 7 (n+r+ 2 


5,282; 0)=2P- = = ae Dee "Antprı2 (0) A(2)=J,(2)]e” (3a) 
ee) u 
2 


# 


Da sich diese Summe über unendlich viele Glieder erstreckt, so wäre ihre direkte Berech- 
nung sehr umständlich und mit der bisherigen Umformung nichts. gewonnen, falls es nicht 
gelingt, den Reihenwert in endlicher Form zu bestimmen. Das wird nun aber dadurch möglich, 
daß man.in (3a) zunächst die Zahlen D durch ihre Integraldarstellung (2) ersetzt und die hier 
erlaubte Vertauschung der Reihenfolge von Summation und Integration vornimmt. In der auf 
diese Weise entstehenden Beziehung (4) machen wir für die Summierung der in der geschweiften 
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Klammer stehenden unendlichen Reihe von der Integraldarstellung (5) für die Funktion A, (z) 
Gebrauch und erhalten dann für diese Summe den Ausdruck 
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der sich unter Zuhilfenahme einer bekannten Transformationsformel für die Ku mmer sche 
Funktion ‚f, auch in der Gestalt 
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schreiben läßt. Nun ist aber offenbar 
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so daß wir nunmehr mit Hilfe der vorübergehenden Abkürzung ß (p, A, r) von Gl. (6) für S, (, 0) 
die Darstellung (7) erhalten. 
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RR SE EERN Lie, 
re -ap!T(p ar PL (a) >“ E Bip Ar) Anttrrn Ve 2a). (m, 


=Ul _|). r=0 
a) 


Die innere Summe von (7) erstreckt sich nur noch scheinbar von r = 0. RE denn man 
erkennt leicht aus (6), daß die Hilfsgröße ßfür aler=p+ 1,p + 2...identisch verschwindet, 


ip u 
ER ie h 


a ausdrücken. Ih Ki 
0 er ve; Br 
4 va Macht man zunächst allein von der Se Tatsache: Gebrauch, so Jap ser RT Er en 
N “ ö naheliegenden Reihenumformung EN f | Bi: Bi: B 3 Is 
R- 5, (2, 0)= 


"A (» +- 5) 20 (— z)r Ph 21 

SEE Ba gm: . Ansaeil e— 22 Br en 
u ae Ka r(;) um) Bu 

| Ersetzt man nun hierin in der inneren Summe ß (p, A, n — 4) durch den a (8) un Kai kehrt 
wiederum die Reihenfolge der Summationen um, so läßt sich die neu entstehende innere 


2% Summe durch die hypergeometrische Funktion 5, (- 2, — 6, 1) zu dem Summenwert 


BEE r(; +?+ s) Ir ri5 + p)r (5 + N) summieren. Nach geeigneten Definitionen entsteht auf 


r diese Weise schließlich für die Funktion 8, (2 „2; 0) von Gl.(3) die endgültige Darstellung (9) 
e mit der Erklärungsgleichung (9a) für die Koeffizienten B,(n). Die Zahlenwerte für die ersten 
neun dieser Koeffizienten sind in der Gl. (9b) angegeben. 


2 0)= I mon. m) Ans (VE ERE) a . 
Wet 
Bun) = 2". or Be an Ey Ve Ma) ; 
5 | 

B, (0) = | 
B,0)=— B, = —2 a--5 |. | 
B,0)=27 Bl)=116 B,@Q)=62  B,ß) — B,()=2 | | 
Ist o= 2» öund also z= »/YRr, so folgt aus (3) und (9) im besonderen die Gl. (10). Mit 


EN Sun (dew je: mer en 


Pin+p+: 


den Gleichungen (3) und (9) für X) ist das erstrebte Ziel erreicht und in (9) für die Größe S, 
eine jedenfalls in den ersten Glie dern bequem berechenbare Formel gefunden. 


4.2. Die Entwicklung einer beliebigen Funktion der Wurzeln von A,=0 in bezug auf » 
nach geradzahligen Potenzen von (R—r)/2 YRr 
Da die Herleitung einer solchen Entw icklung auf dieselbe Weise bewerkstelligt wird wie 


in den Abschnitten 3.1...3.3, so können wir uns hier entsprechend k ar 
urz fas r 
rechnung ist es am bequemsten, von der Funktion i gen. EDLER 


2» 


oo ip \ 
Sy ale) mit Rp = I Bot) (E-ay" Antosınla) =) (1) 


vu wa 
auszugehen. Da für = 0 


so(2) =bu(2) = An (a) = (2j/m)t8- Sin 
2 


a 
Man= 
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ist, so liegen für genügend kleine Werte von i die Wurzeln & = x, (t) von f(x,t) =D in der 
Nachbarschait von 2 =nm mit m=1,2,3... Anstatt nun lediglich die Entwicklung von 
X (t) nach i selbst anzugeben, bestimmen wir an dieser Stelle etwas allgemeiner die Entwick- 
lung einer beliebigen Funktion F (x, (t)) dieser Wurzeln. Dabei soll # (x, (t)) eine im Kon- 
vergenzbereich von (1) reguläre analytische Funktion von 2, (£) sein. 

Bedeuten dann die O0, (x) dieselben Größen wie im Abschnitt 3.2, so ist wieder nach 
2 un/desl 0.1 14] 


1 (am-+) y 
ne 7 
la): a Fe vi )Cy'a)- BRRueı.: 


Die Berechnung der Integrale in (2) erfordert selbst in den ersten Gliedern bereits einen ziem- 
lich großen Rechenaufwand, im übrigen bietet sie aber keinerlei Bessere Schwierigkeiten. 
Man findet schließlich den folgenden Ausdruck: 


“ ı 


Fin) = 


F(2,(8°2))— F(am) — 6°: Bam 4 


m 


P"(gm) [25 1 1 
ie EL IE TEL ee Be eu 22 —. 03,.(3232 . 772 m? —1,2,3.:. 
Eu mE la + 30m (429 — 287 mw) +, 0% (3232 — 303 rn? m?) (m=1,2,3.:.) ir 
ne (3315— — 3257? m?4 It mt) + 04 (8315 — 39052®m® + Brtms)| +0 (89) 


On (08 arm?) [me m8, e=2(R—r, ö= (R—r)/2YRr 


Setzt man in (3) F(x) = x und also F’(n m) = 1, so geben in (3) die von ö abhängenden 
Glieder die an ” m anzubringende Korrektur an, die infolge des Auftretens der über p=0 
hinausgehenden Entwicklungsglieder in (1) erforderlich werden, falls die Gleichung f (z, it) = 0 
bis zur Ordnung ö? noch erfüllt sein soll. Wählt man F(x) = 0% — 28, so ist F’(nm) = — 2m 
und die rechte Seite von (3) stellt dann die Größe 4 »? ö?2 in Abhängigkeit von 6? dar. 

Mit Hilfe der Beziehung 1.(3) kann man in ganz ähnlicher Weise die Lösung der Aufgabe 
bewerkstelligen, die Nullstellen der Gleichung 4, ((R 2)’, (r 2)') = O in bezug auf v zu bestimmen. 
Aus Raummangel führen wir von diesen Rechnungen nur die der Gl. (3) entsprechenden Schluß- 
formel (4) an, die in der Umgebung jeder der einfachen Nullstellen = nm für 6 — 0 mit 

4 2 992 
212.10), Flam) u 62: en fl on ee ci | 


Fam) [68 1 


Fe ki 


a LA Om: (699 — 462 m?) 08, (1537 — 139 0° m) (ms) 
(mm 


3 2 = (4) 
4 x om (4585 — 475 70° m®— 2 et m?) + : Om (4389 — 490 n? m — 3m mi) +0 (68) 
5 


Om = (0 — 2 m?) In? m? R=2lh-r) 6= (R—r)[ARr 


= 


m=1,2,3...gilt. In der Umgebung der Stelle x = 0, die hier für d = 0 ebenfalls eine Null- 
stelle darstellt, besteht eine besondere Entwicklung, auf deren Wiedergabe verzichtet werden 
möge. Über die Bedeutung von F(&) in dieser Formel gelten die gleichen Bemerkungen wie 
zu Gl. (3). 
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Herleitung der Näherungsformel von Laplace für < 
Binomialverteilung, ohne Grenzübergang 

ta Von Paul Lorenz in Berlin ni 
Eine für die Anwendungen wesentliche Bigenschaft der L B, 27 5 schen Näherungsformel ER 


auch für Binomialverteilungen nichtgroßen Umfangs brau i 
Kerleung über ins Unendliche wachsenden Reihenumfang ist daher unerwünscht, und ist auch vermeidbar, 


Er 


’ 
wie durch geometrische Betrachtungen und mit Hilfe der Differenzenrechnung gezeigt wird. Die Methode 
liefert DR die Näherungsformel für beliebig viele Dimensionen. Die Auswertung des auftreienden. viel- > 


fachen Integrals gelingt durch drei Reihen von Substitutionen. 


The approximative formula by Laplace has a property that is essential to applications: It yields — 


still useful approximative values for binomial distributions of rather small extent. The usual de 

based on series of infinitely increasing ewtent is Iherefore not desired and can be avoided, as it is shown 
by geometrical considerations and by use of the difference calculus. The method yields the approzimative 
formula for an arbitrary number of dimensions loo. The computation-of the forth-coming multi-dimen- 
sional integral is carried through by ihree sequences of substitulions. 


Pour les applications de la formule approzimative de Laplace importe son propriäit& de donner 


des valeurs approchees utiles encore dans le cas de distributions binömes non tres ötendues. Sa derivation 


habituelle @ l’aide de series d’une &tendue croissante A l’infini n’est pas desiree et se peut Eviter: l’auteur 
le montre par des considerations geometriques et s’aidant du calcul de differences. Aussi on arrive A la 
formule approximative pour un nombre quelconque de dimensions. L’integral multiple, ce qu’on y trouve, 
est &value par trois sequences de substitutions. 

Ilpnönnxennas dopmyaa Jlanıaca nMeer crenyimee cyIlecTBeHHoe AA IPHAIOMeHnH 
CBOÜÄCTBO: OHA NAeT NOCTATOYHO TOYHLIE IIPHÖNKSKEHHBIE 3HAYeHHA H JUIA ÖHHOMHAJAIBHOTO 
pacıpepereHng HesHayHTeIbHOTO IIPOTAKeHuA. ÜÖsIyHBIÄa BEIBON, OCHOBEIBAWINHÄCH HA 
ÖecKOHeYHOM BOSPACTAHHM IPOTAIKeHHA PA]NOB, IIO9TOMy HeskesratejieH; ETO MO3KHO 
usÖe:KATb, KAK 9TO TIOKA3AHO IIPH IIOMOIIM TEOMETPHYeCKUX PAaccyKkleHHÄä H IOCPeACTBOM 
yPasHeHnü B KOHeYHBIX PasHocTax. IlpnBerennsrä sech MeToA NaeT IPHÖNHSKeHHLIe POPMyY-EL 
H MIA 1R000BOTO YHcCNAa USMepeHHn. Ilossısormmüca MHOTOKPATHBHÄ HHTErTpa.I yaaeTca pemHTb 
HpH MOMOoIIM TPpex PAAOB HONCTAHOBOK. 


Unsere Aufgabe lautet, für das allgemeine Glied der Entwicklung 


N n! n 
u Kr Sa Er > 
N,=U) 


=) 
wobei p, p, und nfest, n, und n, veränderlich, und die Bedingungen , +9, =lunddn +, —=n 
erfüllt sein sollen, einen handlicheren Näherungsausdruck zu finden. 
Wir betrachten dazu das geometrische Gebilde, das durch die Gleichung 


n! 
A engen ch ee Aa 


in einem rechtwinkligen kartesischen Koordinatensystem &,; 2,; y definiert ist, wenn die Be- 
dingung n, +n,=n fallen gelassen wird, und n, —=np, + &, N=Nnp; + 2, gesetzt wird. Es 
ist, da für n, und n, nur ganze positive Werte in Betracht kommen, ein Netz von diskreten Punk- 
ten, die alle in einer glatten krummen Fläche liegen. Ich werde für dieses Netz der Anschau- 
lichkeit und Kürze halber das Wort Fläche‘ im übertragenen Sinne gebrauchen. Ihre Pro- 
jektion auf die &,; z,-Ebene ist ein quadratisches Gitter, dessen Gitterlinien zu den Koordinaten- 
achsen x; und x, parallel sind. Partielle endliche Differenzierung von y (x; %,) nach z, ergibt 


. 1 
A N ERROR TIL AL =. Pı 
1% (25 %) eg ji en + +1) Inp, —n,Inp, =In en 
1 
in 3) 
pı N +1 1 (« 22) (8 a ( 
—= In ———- = 11 — =—Inn— 14) — [2 - 
mr +1 nes l np 2\ np 3\npı a: j 
nPpı 
Die Entwicklung ist konvergent für —1< A s1. Wir gewinnen aus ihr die gesuchte 


nPı 
Näherungsgleichung auf einfachste Weise, nämlich dadurch, daß wir alle höheren Potenzen 
weglassen, 


& +1 
N 


412 (05%) =—Inn— 
Hieraus folgt durch unbestimmte Kumulation nach x, 
(2 + 1)® 


dr —alan- 4m) 27 


Pa 


a 
N 
| 


pr r 
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und wegen der Gleichheit des Aufbaus der Funktion yin n, und n, bzw. x, und 2, sowie in 9, 
und p, steht nichts im Wege, zu definieren 


(2 +1®% (& 719 
2npı 2nps 

Das ist, so lange man für (€ die Möglichkeit der Veränderung offen läßt, die Gleichung 

einer Schar paralleler Flächen. Zwischen ihnen liegt, ihnen nicht unähnlich, die durch die 


Gleichung (2) beschriebene ‚‚Fläche“. Um das Gebilde zu untersuchen; schneiden wir es auf 
durch den ebenen Schnitt nn, + nz =n, das ist auch + ,=0. 


Die Schnittfigur der ‚Fläche‘ ist eine endliche Folge diskreter Punkte, die, bzw. deren 
seitliche Projektionen, den log nat der Verteilung (p, + Ps)” abbilden. 


Die Schnittfigur der Flächenschar ist eine Kurvenschar. Für ihre Projektion in die 


2 (25 %)=— (4, + 9%) Inn — 


A 1Or ARE San 


2; z-Ebene gilt nach der Substition &, =— 2, = die Gleichung | 
(2 a 
(H1)®, 841) 2 
= — — ErEn — [6 =— ee ern Y 
2npı 2nPe 2nPı Pa “ (9 
und entlogarithmiert 
(« a) 
2 

RUE de RR N (8). 


e®—ge 


Aus der Kurvenschar (8) greifen wir die Kurve heraus, welche ein ebenso großes Ebenen- 
stück überdacht wie die Punktfolge, wobei wir uns dieses aus gleichbreiten Rechtecken mit den 
Höhen W zusammengesetzt denken. Wir erhalten, indem wir die seitlichen Projektionen. der 
Ebenenstücke gleichsetzen, für g die Bestimmungsgleichung 


D—P:\ 2 
A ze) 
| ER RERRN BL (9); 
—00 
woraus in bekannter Weise folgt 
1| 


e Vn2 nPpı Pr’ 


so daß die Gleichung der Projektion der herausgegriffenen Kurve lautet 


1 SD eg 


REP TER ER en Re er ET AR Se ZORUNe 
a2 npı pa KR 
Das ist, wenn man noch das Glied EZ im Exponenten wegläßt, die Laplacesche 


Näherungsformel, die wir somit ohne Grenzübergang gefunden haben. Für nicht sehr große 


PıTzPe bei, denn es verbessert die An- 


Reihenumfänge n behält man jedoch besser das Glied — 5 


näherung nicht unwesentlich. 


Es muß nun noch gezeigt werden, daß die Formel (10) innerhalb eines nicht zu kleinen 
Bereichs um den Scheitel wirklich als Näherungsformel für die Projektion der Punktiolge 1 in 
die &,; z-Ebene brauchbar ist. 


Zunächst kann man in bekannter Weise leicht zeigen, daß die Streuungen der beiden Ge- 
bilde gleich sind, und ferner mit Hilfe der Stirlingschen Näherungsformel für die Fakultät, 
daß ihre maximalen Höhen näherungsweise gleich sind. Wir wollen weiter zeigen, daß der senk- 
rechte Abstand der Laplaceschen Näherungskurve von der Punktfolge innerhalb eines gewissen, 
praktisch bedeutenden Bereichs um den Scheitel in bestimmten, im voraus festsetzbaren Grenzen 
bleibt. 


Man gelangt von einem Punkte x,; ©; y der Punktfolge zu einem Nachbarpunkt dadurch, 
daß man eine der unabhängigen Koordinaten um 1 vergrößert, die andere um 1 verkleinert. 
Die Nachbarpunkte haben mithin die unabhängigen Koordinaten &, +1; 2, —1 und 2, —1; 
2, + 1. Wir wollen die Zunahme des senkrechten Abstands A zwischen der Punktfolge und der 
herausgegriffenen Kurve von der Stelle &,; &, bis zu der Stelle x, + 1; 2, — l ermitteln. Es sei 


24 


= 


en An) AA +) 


und speziell für „=-—, = 


Br Bari ee. ek -]- - ( 2 
4 A(a,— a) -[-5( np “r np ae 1 3 ve D 5 


' den zweiten mit II, so gilt 


nen min ng Er a 
aaa; )=Amı +1; BD Aain) = | 
tale; Bel lan Hm Amt im 


=[4,2(9; 9) — 4, y(25 ))—[A2n +1; gr 
1a +T Ya ar = Zu 
-; 2) Br) et; ©. 


Mithin gilt für den Abstand der beiden Schnittfiguren über der Stelle 2; —a 


Au) —-AG—d-— K "E - aa +—.: | | ; a | 


Re) . iy Tre Be 


wobei 0<®9 < 1 sein spll. 9; —d sind dann die Koordinaten des höchsten Punktes, bzw. eines 
der beiden höchsten Punkte, der Verteilung (p} + p,)” oder eines Nachbarpunktes, und 
|4(#; — d)| ist ein kleiner Wert, der meist unbeachtet bleiben darf. * 
Man kann den Wert des Ausdrucks (13) leicht in explizit anzugebende Grenzen einschlie- 
Ben. Der Spielraum wird um so enger, je mehr Glieder man in den Kumulanden berücksichtigt, 
jedoch werden die Ausdrücke für die Grenzen gleichzeitig um so umständlicher. Wir beschränken 
uns deshalb auf eine .einfache Schätzung. Bezeichnen wir den ersten Kumulanden mit I und 


a 


ie 2 “; Wr = k 
PER FRAME. Fr J 


- | 
> £ R 
aa 


al ee | - ) 1 (Er, 

in u ee re FE Sa 7 

2\npı 2\npı ER nPı „e+1= = nPpı 
Bery* 

Betrachtet man nur solche Werte von £, für die 


2 AL ; eh 
Be 2 1, so gilt 3m sIs-\——) ..... (14). 


Entsprechend gilt für 


SaE” 


ya 


I I En Ar 

0 ee ee er 

se s— 3, ) sr <( -) U > (15). 
Aus (13), (14) und (15) folgt 


rer!) +3.R (2) SA na; — 9) Ss 52 a (2) +&(l) (0) | 


e22; s\nPpı 2-0 \nPp, np, = 
das ist 
1 # E A © 
pm: EE-MRE-N—PElEHN)(2EH+ 1) sA (2; — a) — A(9; —Ö) 
f a7). 
S EnannmbRrerZ IR mEe+NDA@E+HN], | 
oder D 
1 2 2 2 2 > 
mp melt a Ze A En ae MP), s Au; —n) — 49; —9) 
1 B (18). 
SHE, 
Für © <0 findet man eine entsprechende Eingrenzung von Al; — 2) — Ad; —9), 


indem man in (17) p, mit 9, vertauscht und & durch — x sowie ® durch 1— ® ersetzt. 


Beispiel: 


3 
BT Be n=144. 


P fr r >: “ r 
P Hai) . 4 
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Die Verteilung hat ihren größten Wert bei n, 2 np, = 36; n, = np, —= 108. Er ist 0,0766. Es 
gilt 2, —n, — 36 und 2, =n, — 108; ferner 9 —=0 und y(0;0)=In 0,0766 — — 2,5692. 
Die Näherungskurve hat die Gleichung 


Een 
108 

Die Eingrenzung lautet für positive @ 

— Googr SE + 1522 +40) + 0,0012 < Als; — a) < 


und für negative & 


a — 2,5680 ren mr. 


1 


a er x) + 0,0012 


1 
Zagar (82° — 152? + 4x) + 0,0012 < Az; — a) 


69984 = 102976 (26 # — 42 2? +13 2) + 0,0012. 
Diesen Eingrenzungen entsprechen die äußeren Kurven des Bildes, und zwar für die Abszissen 
von =—12bis 2 —=15. Die Werte sind für die bildliche Darstellung entlogarithmiert. Die . 


innere Kurve zeigt das Wachsen des Abstandes der Näherungskurve von der Punktfolge, beim 


SRG (1 ae). Peer a ee DE Et UNEe DE He aut a a / Baer PA 


Entfernen vom Scheitel, gemessen an den Ordinaten der Punktfolge. Der Abstand beträgt an 
den Grenzen der Zeichnung 234 und —295°/,0.. Gemessen an der Scheitelordinate ist er 14,8 
und — 5,7%) g0- 

Allgemein lehren die Eingrenzungen (17) und (18), daß der Abstand zwischen der La- 
placeschen Kurve und der Punktfolge in der Nachbarschaft des Scheitels absolut klein ist 
und beim Entfernen vom Scheitel unter kleinen Schwankungen absolut wächst, jedoch bis zu 
einer bestimmten Abszisse unter einem festen Werte bleibt, so daß die Laplacekurve in einem 
je nach dem geforderten Grade der Näherung abzugrenzenden Bereiche als Näherungskurve 
der Verteilung (9, + 95)” dienen kann. 

Die Betrachtung gewinnt an Anschaulichkeit, wenn man eine neue unabhängige Variable 


4 


% -=- einführt, wobei o =] np, p, die Standardabweichung der Punktfolge und also auch der 


Laplaceschen Kurve sein soll. Die Eingrenzung (18) geht dadurch über in 


1 2 2 BEN 2 2 & 2 au t, ’ 1 n 
2,2 mp) o Vz u u RE ne er u 
1 Er ee... w re) 
sn, BRPT Br ar Dil, 


Durch diese Substitution beschränken wir den Spielraum der unabhängigen Variabeln, 
der für die Anwendungen allein von Interesse ist, auf die Werte von —4 bis 4. Anden Grenzen 
dieses Spielraums sind die Ordinaten der Näherungskurve nämlich bereits kleiner als 0,40/,, des 
Größtwertes und werden mit Zunahme der Entfernung schnell noch kleiner, praktisch also null. 
“ Die Punktfolgen und ihre Näherungskurven, gleichgültig welchen Reihenumfangs, n, werden 
für festes p, und p, alle „‚gleich breit“. Macht man sie auch noch gleich hoch, etwa, indem 
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Näherungskurven der Gleichung 


x \ \ y g | } 
2 > | Balz ) er 
N u=he | k 


x 


: und ihre Bilder unterscheiden sich nur ein wenig durch die Lage. Die 


2° sich bereits für.mäßig große n ihrer Gestalt nach nur wenig, aber die Folge der Punkte wird mit 
vuee wachsendem n immer dichter. Etwas anschaulicher wird die Vorstellung, wenn man sagt, daß 
Br. die Ordinaten mit wachsendem n immer mehr zusammenrücken. (Vergleiche dazu v. M 5, 
Bi: Wahrscheinlichkeitsrechnung, $. 135.) Wählt man für p, und p, ein anderes Wertepaar, 


werden die Gebilde alle entweder breiter oder schmäler. { 4 

a Die Eingrenzung (19) lehrt sehr anschaulich, daß beide Grenzen für den Abstand zwischen 

Fa der Laplace kurve und der Punktfolge bei festem p, und p,, aber veränderlichem n, und also 

Be) auch o, an einer bestimmten Stelle @’ absolut um so kleiner werden, je größer n wird, woraus 

5 “folgt, daß die Laplacekurve für jedes Wertepaar p, und p, als Näherungskurve benutzt 
werden kann, wenn es nur möglich ist, n genügend groß zu machen, auch für extreme Werte 
von p, und ?%. 

"hs een angebracht, diese Behauptung ‘durch ein Beispiel zu belegen, in anbetracht 
der oft in der Literatur zu findenden Meinung, daß in solchen Fällen allein die Poissonsche 
Näherungsformel brauchbar wäre, Ir 

Ist 9, = 0,05, 9, = 0,95, n = 8421, soist o=20, und für den höchsten Punkt der Punkt- 
folge gilt n, = 421, 2 = m —npı = — 0,05. Wir setzen 9 = 0,95, so daß 9’ = 0,0475 wird, 
und überzeugen uns leicht, daß die von © abhängigen Größen in der Eingrenzung (19) unbe- 
achtet bleiben können außer für absolut ganz kleine &, für die eine Eingrenzung ohne Interesse | 
ist. Die Eingrenzung lautet dann, mit kleinen Abrundungen, für positive =’ 


—0,000009375 (800 #’3 + 60 2’? +2’) < A(#’; — x’) < — 0,000006215 (800 2’? + 60 2’? + =) 
und für negative x | 
0,000009375 (800 #3 — 60 #2? + x’) < Ala’;— x’) < 0,000018785 (800 z’? — 60 «+ =’). 


Hieraus ergeben sich die folgenden Zahlenwerte für die Eingrenzung der verhältnismäßigen 
Abweichung der Näherungskurve von der Punktfolge, in Promille 


I. der Ordinaten der Punktfolge II. der Scheitelordinate | ” 
2 —=1;—1 —8 und; +7 und-+ 14 —4,9 und —3,2; +4,2 und + 8,2 | 
u =2; —2 —60 und —40; +59 und +123 —8,1 und —5,4; +8,0 und +16,6 
& =3;—3 —188und —129; -+218 und +485 —2,1und —1,4; +2,4und + 5,4 
X —4;,—4 —0O,1und —0,1; +0,2und + 0,6. 


Man darf schließen, daß die Wiedergabe der Verteilung (0,05 + 0,95)%2 und ähnlicher 
Verteilungen durch die Laplacesche Näherungsformel bei den Anwendungen in vielen 
Fällen brauchbar sein wird. 


Mehrdimensionale Verteilungen 
Die im vorstehenden dargestellte Methode der Untersuchung von geometrischen Gebilden, 
die zunächst nur für diskrete Punkte definiert sind, erweist sich als besonders vorteilhaft, wenn 
man sie auf mehrdimensionale Gebilde anwendet. Wir wollen sie auf das geometrische Gebilde 
anwenden, welches durch das allgemeine Glied der Entwicklung von 


MR: Fer (1), 
das ist 

= REBEL N n nn 
mim emgihi Pr’. » . Py® De Te - (2) 


beschrieben wird, unter der Voraussetzung, daß n p, P2. . . px fest, aber N Ng... N, unabhängig 
veränderlich Sind, und die Bindung 9, +9 + "+ pr = 1 besteht. Da die n; nur ganze 
positive Werte annehmen können, besteht das Gebilde aus lauter diskreten Punkten und ist 
nach allen Seiten begrenzt. Die Projektion der Punkte in den Grundraum Ni5Ng5...Nzistein 
rechtwinkliges Gitter, dessen Gitterlinien zu den Koordinatenachsen parallel sind, und man 
kann folglich von einem Punkte. des Gitters zu jedem andern kommen, indem man sich aus- 
schließlich parallel zu Koordinatenachsen bewegt. Der Bewegung längs einer einzigen Gitter- 
linie, etwa parallel zu der Koordinatenachse n, entspricht analytisch partielle endliche Diffe- 
renzierung nach n; (wofür noch ein bequemerer Ausdruck gefunden werden muß) und ihre Um- 
kehrung, die Kumulierung, längs dieser Linie. 
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Wir betrachten wieder, wie bei der Herleitung der Laplaceschen Formel für eine Di- 
mension, In W, setzen ihn gleich y, und differenzieren y endlich nach »,. Wir erhalten genau 


wie dort, 
Ay-mn at fact Imst burn = SeMalte ; .. 8) 


en ap AP, 


wobei u =np; + x; sein soll. 


Die Entwicklung ist konvergent für —1< 2) <1. Wir gewinnen aus ihr die ge- 
suchte Näherungsgleichung auf einfachste Weise, nämlich dadurch, daß wir alle höheren Po- 
tenzen weglassen, 21 

% 


A255 EM AN NE pp RR 
N Pı 
Hieraus folgt durch unbestimmte Kumulation nach x, 
(& +)@ 
nenn, a en ET 5) 


und wegen der Gleichheit des Aufbaus der Funktion y in den n; bzw. x; sowie in den 7; steht 


nichts im Wege, zu definieren 
Ds +9 HDD, + N“ 
= —- (+++ )an— Inpı REN DT HE 0 Ab): 


Das ist, solange man für C die Möglichkeit der Veränderung offen läßt, die Gleichung 
einer Schar paralleler krummer k-dimensionaler Räume. Zwischen ihnen liegt, ihnen nicht un- 
ähnlich, der durch die Gleichung (2)beschriebene, aus diskreten Punkten bestehende k-dimensio- 
nale „Raum“, Um das Gebilde zu untersuchen, schneiden wir es auf durch den ebenen 
Schnitt , + n+ "+ m, =n, dasist auch ++ '' + =0. 

Die Schnittfigur des ‚Raumes‘ ist eine (k—1)fache Mannigfaltigkeit diskreter Punkte, 
die, er deren ‚‚seitliche‘‘ Projektionen, das heißt die anne in die ebenen Räume %; 
%35 "2%; W und so weiter, den In W abbilden. 

"Die Schnittfigur der Raumschar ist eine aeunehdiiche Schar krummer (k — 1)-dimen- 
sionaler Räume. Ihre, aus (6) hervorgehende Gleichung lautet nach einer geringen Umformung 


4) er. 


2= n u olng a er 


2npı 2nP 2nPr 


k (si4 2)? 


w=d=ge 4-1 


und entlogarithmiert 


Aus der Raumschar (8) greifen wir denjenigen Raum heraus, welcher ein ebenso großes 
Stück ebenen k-dimensionalen Raumes überdacht wie die Punktmannigfaltigkeit, wobei wir uns 
dieses aus ‚‚k-dimensionalen quadratischen Säulen‘ mit den Höhen W zusammengesetzt denken. 
Wir erhalten dann für g ” a 


EB J 1 2 er | 
en 2 „3 2N Ds 2np 
d.xg 2" de dl ml 


Das Integral ist nur (k — 1)fach, weil die Zahl der unabhängigen Veränderlichen zufolge 
der Bindung 5, +++ =0 nur gleich k — 1 ist. Zur Vereinfachung der Schreib- 
weise habe ich die bisher a hltene Reihenfolge der x; umgekehrt, und betrachte das dabei 
an die letzte Stelle kommende x, als durch die Bindung bestimmt. Ferner ist klar, daß man 


die Glieder + z im Exponenten weglassen kann, denn eine Parallelverschiebung ändert offenbar 


nichts an dem Problem. Schließlich habe ich bei den Integralzeichen die unteren Grenzen, die 
alle —®, und die oberen, die alle & sind, weggelassen. 

Die Integration gelingt durch drei Reihen von Substitutionen, von denen die ersten beiden 
lauten 


mt%=L, Ppı t+P2 = Pr 
Grm de P,-+ 93=P; 
Se (10). 
Aa + = K-fdı Pıa + mMi=Prı 


Kıta=Xm0 Pat m= P=]1 


ee 


Indem mit dem Integral 9 identischen Integral 11 


z al. ® Sr S.j® eig un "di, | N ee: 
x Erreicht der, nur a einer Veränderlichen abhängige Integrand unter dem letzten nte | 
zeichen ein Maximum an der Stelle &, — Xp; n=ÄL, B Diese Stelle mache Ich ZURE ESSEN B en i ö 


einer neuen Veränderlichen t, gemäß der Substitution ee 
5 “ 5 Sr 
1=—htäp und a eh 
Dann wird 
RE BELEN Re TER | 2 + En pn) 
anp .2nPpı 2nmppıl’ \P 
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und folglich das letzte Integral von (11) 


X. Dres x. Ver yr 
vr =, ER nn I TER TE 
8 | 3 zn,pıP dt, je unPp,; F, 


Dadurch wird die Zahl der Integrationsvariabeln um eine vermindert, aber das Gesamtintegral 


n2npıp RR ee ar, 
| een ar 1 In pr 2nP nf day: -daz 2. a (12) 
behält seinen Aufbau, denn es gilt , + "+9 #P,;=1und u +" +2+ Be er 
ya n2np, Pa | 
Folglich darf das Verfahren wiederholt werden. Dabei tritt nun der Faktor Bag | 


3 
vor die Integralzeichen und vereinigt sich mit dem bereits dort befindlichen Faktor zu 
(mn)? PpıPs Pa 
er,D 1 


In dieser Weise kann man fortfahren, bis das Gesamtintegral vollständig 


ausgewertet ist, und findet schließlich als Gleichungen der entlogarithmierten Näherungsformel 


Sets) 
Leine Be (13). 


Mr 2 net pıPp2... Pr 
: P+mtttm-]l Hart +rm=0. 
Die Formel bleibt gültig, wenn man sie auf die beiden miteinander verbundenen Variabeln x, 
und z, beschränkt, und geht dann durch die Substitution & =— 2, =zin die Laplace- 
sche Formel über. y 

Wie für diese kann man auch für die allgemeine Formel (13) leicht zeigen, daß ihre Streu- 
ung gleich der Streuung der zugrundeliegenden Punktmannigfaltigkeit ist, und ferner mit Hilfe 
der Stirlingschen. Näherungsformel für die Fakultät, daß die maximalen „Höhen“ der 
beiden Gebilde, d.h. Ymax und Zmax näherungsweise gleich sind. Weiter folgt, wenn man z.B. 
%35 %45 -» .» % festhält, genau auf dem im ersten Abschnitt beschriebenen Wege, daß in einem 
gewissen, praktisch bedeutenden Bereiche um den ‚‚Scheitel‘‘ der Abstand zwischen der Punkt- 
mannigfaltigkeit und dem Näherungsgebilde in bestimmten, im voraus festsetzbaren Grenzen 
bleibt. Daraus folgt aber, daß in einem ganzen, (k—I)fachen, praktisch bedeutungsvollen Be- 
reich um den „‚Scheitel‘‘ die Wiedergabe der Punktmannigfaltigkeit durch das Näherungs- 
gebilde gewissen im voraus festsetzbaren Bedingungen genügt, denn man braucht, um vom Schei- 
tel oder von einem gewissen Punkte 9; ®,;...®;, in der Nähe des Scheitels zu irgendeinem 
Punkte des den beiden Gebilden zugrunde liegenden Gitters zu gelangen, höchstens k—1 Teil- 
strecken zurückzulegen. Längs jeder Teilstrecke ändern sich der Abstand :und seine formel- 
mäßigen angebbaren Eingrenzungen. Der gesuchte Bereich ist der, innerhalb dessen die Sum- 
men aus dem Anfangsabstand und den Änderungen die vorgeschriebenen Höchstbeträge nicht 
überschreiten. Es folgt aus den Entwicklungen des ersten Abschnitts, daß bei festen ?1; 
Ps; * Pr die Wiedergabe um so besser ist, je größer » ist. 


U Ca 


Eingegangen am 2. 3. 1949. 
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Genauigkeitssteigerung belm Beilschneiden- 
planimeter. 


Die Anwendung des von Prytz!)2)®) erfundenen 


ilschneidenplanimeters zur Ausmessung von Flächen - 


ist seitlangem bekannt. Einer allgemeinen Verwendung 
dieses sehr einfachen Instruments stehen allerdings 
einige ganz erhebliche Nachteile entgegen. Es liefert 
nämlich den Flächeninhalt nur in angenäherter Größe. 
Der Fehler wächst um so mehr an, je größer der Durch- 
messer der zu planimetrierenden Figur im Verhältnis 
zur Fahrarmlänge des Planimeters ist. Eine praktische 
Regel besagt, daß der größte Abstand eines Peripherie- 


“ punktes vom Schwerpunkt der Figur nicht größer als 


etwa 20% der Fahrarmlänge sein soll, wenn die untere 
Grenze des Meßfehlers 1—2%, des wahren Wertes er- 
reichen soll. Diese Beschränkung ist aber nicht in jedem 
Falle gegeben, zumal der Fehler auch noch stark von 
der Form der zu planimetrierenden Figur abhängig ist. 
Im folgenden soll daher ein Verfahren angegeben wer- 
den, das für beliebige Verhältnisse sehr genaue Resul- 
tate zeitigt. 

Das Gerät, welches mit dem Fahrarmstift die Punkte 
P der gegebenen Kurve y = f(x) durchläuft (Bild 1), 
beschreibt mit seiner dem Fahrarm parallelen Schneide 
die Schleppkurve y* = f* (x*), so daß in den Punkten 
P* die Richtung des Fahrarms gleichzeitig die Richtung 
der Schleppkurventangente darstellt. Esgilt also, wenn 
die Fahrarmlänge ! als Maßeinheit gewählt wird, für 
den Bogen P P’inkartesischen bzw. in Polarkoordinaten 


dpo=sinpde —cospdy bzw. 1 
dp = sin (p — 9) dr —rcos(p — d) dd \ N) 


ey 
und für die Umfahrung des Elementarsektors S,P,P’ die 
Differentialgleichung *) 5) 


do —d) ___ = 
N UL TE [4A + Bcos ( —®)], 


(2) 
wo A=koir—r ©inr; B=rbofr —Einr 


Diese Differentialgleichung ist nur für einige Sonder- 
fälle elementar geschlossen lösbar. Insbesondere ergibt 
sich für den Kreis, mit Einschluß der Geraden, bei der 
Umfahrung S-P-Bogen-P-S das Integral 


ayl—r?= ern Ä 
nl Po 
VeV0 & 
d.h. tg =tg— = 
? : (3) 


tgay1l—r? 
A-+ Becos 

2 ( Po) Yin 
ae er 
ee ber: 
yı-r: 


Im allgemeinen Fall läßt sich (2) auf eine Differen- 


‚ tialgleichung vom Riccatischen Typus zurückführen 


Eine Reihenentwicklung führt auf folgenden Zusam- 
menhang 


2n 2n 2 
ee RE I BE 07% 
=; | + 204 [ra0+ 
) i ) ) da 
7 


)> 


B| 1 
En +3 | Ptn”+ ..|0007 00 
0) 


1) Prytz: Engineering 57 (1894), 687 u. 813. 

2) Willers: Mathematische Instrumente. München 1943. 

3) Meyer zur Capellen: Mathematische Instrumente, 
2. Aufl., Leipzig 1948. 

°s) Runge: Z.f. Verm. Wes.24 (1895), 321/331. 

6) Schreiber: Z.f. Verm. Wes,37 (1908), 689/702. 


d.h. mit R=r/ und F als Flächeninhalt 
Pa=li= tm | Mary 


en RtdF+..: (4b). 


2 ae 
a + pr + \eoswar 


Bild 1 


Hierin ist y der Winkel zwischen Anfangs: und End- 
lage der Fahrarmrichtung. Die Entwicklungen (4) zei- 
gen zugleich, daß der Flächeninhalt in den auf die Fahr- 
armlänge bezogenen Einheiten nur dann einigermaßen 
richtig wiedergegeben wird, wenn die Reihe infolge der 
Kleinheit der Glieder nach dem ersten abgebrochen 
werden darf. Das ist jedoch nur dann gegeben, wenn 
stets R< 1 bleibt. In diesem Fall wird aber auch ß 
klein, so daß schon kleine Fehler bei der Abmessung des 
Bogens ß erheblich ins Gewicht fallen können. Außer- 
dem müssen in (4) aber auch die y enthaltenden Glieder 
zum, Verschwinden gebracht werden. Das läßt sich 
allerdings in bekannter Weise durch Umfahrung der 
Figur mit den Anfangslagen 9, und 9, + r) des Fahr- 
armes und anschließender Mittelbildung erreichen: 
2B = Pw+ Bm +n. — Es sind verschiedene Wege vor- 
geschlagen worden, um die Genauigkeit des Beilschnei- . 
denplanimeters zu steigern, u.a. von Berroth®)in 
der Weise, daß der vorstehende Prozeß mit zwei ver- 
schiedenen Fahrarmlängen /, und /, vorzunehmen sei, 
und daß durch Extrapolation auf | > 00 aus den ab- 
gebrochenen Reihen 


3 1 
= Ftzp | war, | 


= 1 
,ß2= F+ zu | war] 
auf den Flächeninhalt 
ri, _UFR—LR, er (da) 
2 —1} 2 —15 
geschlossen würde. Theorie und Praxis zeigen jedoch, 
daß auch hierbei noch erhebliche Abweichungen ent- 


stehen können, da. der Näherungswert F des Flächen- 
inhalts aus der Form 


en M R® 1 7 R, 
r= ra= 18 = [s(plar=5 |) GRZ (6) 
N) 


der Gleichung (4) folgt. Ermitteln wir jedoch in (6) 
für R einen konstanten Mittelwert R, so folgt 


® 
Fe, (Gem d.h.a=g(Am) -» . (N, 


(5) 


*#) Berroth; Schweiz, Z, f. Verm. Wes. 32 (1934), 1/11. 
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« 2 Hi w N 
Wo m= "5" den auf den Kreis mit dem Radius R. 


Sry Bogen im Einheitskreis bedeutet. Die Be- 
ziehung & = 9 (&m) ist aber für den Kreis sofort zu be- 
stimmen. Nach (3) haben für die Anfangslagen , und 
tm: 


wu. AtBmemde, | 
3 ‚1— Bsino,'@ 


& Pe RE) 
rn . (A—Beos 9) @ 
Br 1+Bsino,'o’ 
tgryl—r? 
wo ee 
ur er) 


BR +8 (2 A+ B2sin2Q,) 0 
somitigda= tg "3 = 771 (48 B!con2p,)e% 
Hiermit ist also & für den Kreis als Funktion des 


Radius mn und damit als Funktion des Wertes 


sı BR? 


DEF eg 
nV 


ben. Diese Werte sind in nachfolgender Tabelle 

für, = 0 berechnet worden. — Im praktischen Falle 

werden wir daher so vorgehen, daß wir nach den Um- 

fahrungen mitg, = 0 undp, = aus den gemessenen 
1 


Bögen ß, und ß, (s. u.) den Wert = 31 (Bı+ß,) er- 


mitteln und hierzu aus der Tabelle das zugehörige dm 


bestimmen. Mit F = I? „erhalten wir dann auchschon 
bei manchen ungünstigeren Fahrarmverhältnissen recht 
brauchbare Werte für den Flächeninhalt. Eine in 
allen Fällen befriedigende Genauigkeit erreichen wir 
jedoch durch Wiederholung des Meßvorganges mit einer 
zweiten Fahrarmlänge, z. B. der halben, und Bestim- 
mung der Werte &,, für jede Umfahrung. 


Tabellel. 
BR 12 1’ - &u-tda 

r 7 Om 5 Di 5 Dr 5 
0,00 0,0000 0,0000 0,0000 0,0000 
0,10 0,0314 0,0158 0,0157 0,0315 
0,20 0,1257 0,0637 0,0634 0,1271 
0,30 0,2828 0,1461 0,1435 0,2896 
0,40 0,5027 | 0,2682 | 0,2559 | 0.5241 
0,50 0,7854 0,4396 0,4028 0,8424 
0,60 1,1310 0,6731 0,5847 1,2578 
0,70 1,5394 0,9780 0,8158 1,7938 
0,80 2,0106 | 1,3342 | 1,1580 | 2.4922 
0,85 2,2698 1,5140 1,4464 2,9604 
0,866= 1 2,3562 1,5708 1,5708 3,1416 
0,90 2,5447 | 1,6873 | 1,9232 | 3.6105 
0,95 | 2,8353 | 1,8412 | 2,5876 | 4.4288 
1,00 3,1416 | 1,9708 | 31416 | 5.1124 
1,10 3,8013 2,2037 3,6715 5,8752 
1,30 5,3090 2,5004 4,0607 6,5611 
1,60 8,0425 | 2,7440 |: 4,3142 | 7.0582 
1,90 11,3413 | 2,8787 | 4,4501 | 7,3288 
2,20 15,2052 2,9626 4,5334 7,4960 


Man wird natürlich nicht immer erwarten können, 
daß die ermittelten Flächenwerte F, = I? %n, und 


F,= 13%, genau übereinstimmen. Extrapoliert man 
aber jetzt mit den (5) entsprechenden Formeln 


PEN : e i 
Hom=P,=F+pö Hom=F,=F+70 


Er b BER? sone Pre nr, 


EU wenn 2: r=th—H 


falls mft inagensmat, 2 Umsfahfungen aus, nänlinS3ENENEE 
un ra, oderp=n. Das vor- 


nde Verfahren also keinesfalls eine wesent- 
er Amen mit sie zumal man natürlich durch 
entsprechend bezifferte Maßstäbe unmittelbar die &- 


Werte ablesen könnte; doch gehen wir auf diese Dinge 
hier nicht ein. 


Du allaee befnälichen Krsisbogena (EM 
fangs- und e £ indli un (Radius 1) 
kann man sich der folgenden rungskons 
von Lambacherjbedienen’). Sind in Bild2 44’ 
und BB’ die beiden senkrechten Durchmesser 21 des 
gegebenen Kreises 
um M und AP der 
gesuchte Bogen, so 
schlage man um A 
den is mit dem 
Radius AB, der A4’ 
in O trifft und dann 
um O0 den Kreis mit 
demRadius 0A. Die 

Verbindungsgerade 
OP trifft diesen 
Kreis in @. Dann ist 
die Entfernung AQ 
sehrgenaugleichdem 
Bogen AP. Bei Win- 
keln 0< a<{ 100° 
übersteigt der Feh- 
ler nirgends die 
Größe von 0,05%. 


Bevor wir die zweckmäßige Ausgestaltung des Beil- 
schneidenplanimeters für das vorstehende Verfahren 
beschreiben, sei an einem praktischen Beispiel seine 
Leistungsfähigkeit nachgewiesen. Es wird dies hier mit 
allen Einzelheiten und Möglichkeiten behandelt. Im 
konkreten Falle wird natürlich immer nur je ein Wert 
aus jeder der nachfolgenden Gruppen benutzt. Gegeben 
sei ein Rechteck mit den Seiten a=5cemundb= 30cm, 
welches einmal mit der Fahrarmlänge !ı = 20 cm, 
zum anderen mit /ıı = 10cm umfahren werde. Im 
ersten Falle mißt man fürp,= 0 bzw. @,= 7 die Bögen 


Brı = 8,33 cm, Bra = 8,19 cm und daraus Pı= 
—Bıı + Bıa) = 8,26 cm, d.h. Fıı = 166,6 cm?, Fıs 


Bild 2 


= 163,8cm?, Fı = 165,2cm? (Fehler + 10%). Die 
entsprechenden Werte für lmı= 10cm lauten Br 
= 24,0lcm, Prıa = 22,6l cm, Bu = 23,3lem, Fin 
= 240,1 cm?, Fıra= 226,1 cm?, Fır = 233,1 cem?(Fehler 
+ 55,2%). Extrapolation auf Grund von (5) liefert F 
= 142,6 cm? mit einem Fehler von —5%. Unter Be- 
nutzung von Tabelle 1 erhalten wir 
&lı == 0,4165, d. Hr Am, = 0,3947 und Fjı = 157,9 cm? 
%12 = 0,4095, ,„ Am,= 0,4045 „, Fra=161,8 ,, 
öı = 0,4130, „ Am = 0,3990 „ Fı =1596 „, 
ferner 
AIm= 2,401, och, Am, = 1,833 und Filı = 133,3cm?. 
alle 2,261 ,„, .am,= 1,973 „ Fır = 197,3 ,, 
811 = 2331 „, om =1,901 „ Fır =190,1 , 


‘) Lambacher:; Z.f.math. u. nat. Unterr.1. 1948, Bd.2, 
3.94, 
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Extrapolation mit Hilfe von (9) ergibt 

aus Fr, und Fir, den Wert ?=149,4 cm? (Fehler -0,4%), 
>> Fig 3 FTIa ”> ” F=150,0 Er) ( ” 0,0%), 
”’ Fi 3 II Er} Er] F=149,5 Er] ( > -0,3%). 


Wegen der gestreckten Form desRechtecks(@a:b=1:6) 

ist dieses Beispiel für die Planimeteruntersuchung 

besonders ungünstig. Wenn trotzdem die Fehler über- 

raschend klein werden, so wird hiermit die besondere 
Leistungsfähigkeit dieses Verfahrens beleuchtet. 

Das Beilschneidenplanimeter in seiner 

gewöhnlichen Form hat eine feste Fahr- 

'armlänge und ist daher für den Gebrauch 

des vorliegenden Verfahrens nicht ohne 

weiteres geeignet. Um die beiden Fahr- 

armlängen darzustellen, läßt sich eine 

Reihe von Vorrich- 

tungen ersinnen. Am 

einfachsten wäre es, 

die Schneide oder 

den Fahrstift am 

Fahrarm mit Hilfe 

| einer Muffein Längs- 

richtung beweglich 

zu machen (Bild 3a). 

Wird. die feste Beil- 

schneide durch eine 

bewegliche Rolle er- 

setzt, so wäre es 

möglich, die heraus- 

nehmbare Rolle in 

Vertiefungen, die im, 

Fahrarm mit u-för- 

migem Querschnitt 

anzubringen ' wären, 

zu versetzen. (Bild 

3b). Neben weiteren 

anderen Lösungen er- 

scheint die folgende 

auch theoretisch sehr interessant. Nehmen wir nämlich 

an, daß die Planimeterschneide um den Winkel egegen 

die Planimeterachse gedreht ist, so gilt für die Schlepp- 

kurve (Bild 4) 


Bild3a Bild 3b 


dk 

ı* = + 1cosp; y*=y+1sino; Tr =tgp+ &). 

Nach kurzer Umformung folgt aus diesen Bezie- 
hungen 

rt =tg(P+ o) 

«—Isinpdo 

lcos edp=sin(p+ e) de— cos (P+ e) dy: (10), 

mit@ + .= y also 


oder 


lcosedy=sinyde —cosydy .. . (11). 


Das ist aber dieselbe Differentialgleichung wie (1), 
mit.dem Unterschied jedoch, daß die Fahrarmlänge im 


Bild 3c 


Bild 4 


Verhältnis cos & : 1 kleiner geworden scheint. Mit den 
Anfangslängen y,=0, d.h. 9, =—e und „=n, 
d.h.9%=n— släßt sich bei dieser neuen (mathema- 
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tischen) Fahrarmlänge I! cos & auch die Tabelle I in der 
beschriebenen Weise benutzen. Ist z.B. &= 60°, so 
wirkt sich das aus wie die Verwendung der halbenFahr- 
armlänge. Montieren wir also für diesen Fall die 
Schneide auf einen gleichseitigen Dreiecksfuß, so kann 
man durch Umstecken dieses Fußes in der gleichge- 
formten Vertiefung des Fahrarmstabes sofort die 
Drehung um 60° herstellen (Abb. 3c). 

Eine weitere Maßnahme, um mit einfachen Mitteln 
genaue Ergeknisse zu erhalten, besteht in einer Um- 
zeichnung der Berandung der zu, planimetrierenden 
Fläche in einen verzerrten Maßstab, Die obige, für den 
Kreis berechnete Tabelle 1 kann auch so aufgefaßt wer- 
den, daß die Werte & als wahre Flächeninhalte betrach- 
tet werden. Das kann man aber nur dann tun, wenn 
die Randkurve der zu bestimmenden Fläche vorher in 
einen kleineren Maßstab umgezeichnet worden ist. 
Offenbar ist das veränderliche Verkleinerungsverhältnis 


.der Radienvektoren ]/ 2. Damit geht dann das gege- 


& 
bene r über in das verkleinerte 0. Nach Umfahrung die- 
ser neuen verzerrten Fläche mit 9,=0undg,=n kann 
dann unmittelbar aus dem erhaltenen & mittels F=1} 
der richtige Flächeninhalt ermittelt werden. Tabelle 2 
gibt den Zusammenhang 0 = o(r). 


Tabelle 2. 
a Maar % 
0,000 0,000 0,400 0,390 0,800 | 0,736 
0,100 0,099 0,500 0,482 0,900 | 0,804 
0,200 | 0.198 | 0.600 | 0.572 | 1.000 | 0,864 
0,300 0,295: 0,700 0,658 1,200 | 0,957 


Das. Schneidenplanimeter ist herstellungsmäßig und 
in der Handhabung ein außerordentlich einfaches Ge- 
rät. Im praktischen Gebrauch ist es von erstaunlich 
hoher innerer Genauigkeit, d.h. die Streuungen der er- 
haltenen Werte bei Wiederholungsmessungen bewegen 
sich in sehr kleinen Grenzen ®). Aus diesem Grunde 
kann es bei Verwendung eines genauen Auswertungs- 
verfahrens wie des vorstehenden an vielen Stellen die 
komplizierteren und teureren Polarplanimeter ersetzen. 
Das ist um so mehr überall dort möglich, wo bei empi- 
rischen Werten der zu planimetrierenden Flächen über- 
triebene Anforderungen an die Meßgenauigkeit nicht 
am Platze sind. 


Elmshorn/Hamburg. H.Athen. 


8) Vgl. Schreiber: Z.£f. Verm. Wes.38 (1909), 8.401/14. 
Über das Quadratwurzelziehen aus kleinen 
Zahlen auf der Rechenmaschine. 

Das Studium von Gauß’ numerischen Rechen- 
methoden (Werke X2: Philipp Maennchen, 


' Gauß als Zahlenrechner) hat mir folgende Anregung 


zur Berechnung von Quadratwurzeln aus kleinen Zah- 
len gegeben. 

Gauß berechnet diese Wurzeln aus ihrer Ketten- 
bruchentwicklung. Er geht dabei in der Reihe der Nä- 
herungswerte so weit, daß er den Nenner noch in Prim- 
faktoren unter 1000 zerlegen kann, stellt dann den 
Näherungswert als eine Summe von Brüchen mit diesen 
Nennern dar und entnimmt jeden Bruch aus einer von 
ihm selbst angefertigten Tabelle, mit der er überhaupt 
dauernd arbeitet. Aber sowohl zur Zerlegung einer 
großen Zahl in Primfaktoren als zur Zerlegung eines 
Bruches in eine Summe von Brüchen der genannten 
Art gehören die Erfahrung und die seltene Kombina- 
tion tiefgehender zahlentheoretischer Erkenntnisse mit 
rechnerischer Fähigkeit, wie sie außer G a u ß wohl nur 
Euler besessen hat. 

Zudem wäre in unserer Zeit der Rechenmaschine das 
Gaußsche Verfahren sowieso überholt. Wir haben 
keine Scheu mehr vor Divisionen durch mehrstellige 
Zahlen, falls nur diese Zahlen nicht m ehr Stellen be- 
sitzen als das Quotientenwerk der Maschine. 


werden, die sich für die Zwecke der Rechen in 
besonders eignet. Die bekannte Methode, um aus dem 


Näherungswert a der Wurzel von x einen genaueren zu. 


finden, nämlich x durch a zu nr und aus dem 


-Quotienten x/@ und a das arithmetis 
den, also. 


Wer safe ren a rekdh 


liefert ein Resultat, das nur bis auf etwa doppelt so 

viele Ziffern genau ist wie a war, und kann daher nur 
"doppelt so viele richtige Ziffern liefern, wie das Quo- 

tientenwerk der Maschine faßt. ; Y 

Will man weiter gehen, so hat man 'von besseren 
Näherungswerten auszugehen als von einem Dezimal- 
bruch «a, dessen Genauigkeit ja schon mit der letzten 
Ziffer aufhört. Dazu bieten sich die Ketten- 
brüche. Man geht aber hier besser.nicht wie Ga uB 

 — wenigstens in den vonMaennchen herangezo- 
genen Beispielen — tut, vor, sondern trifft folgende 
erungen. 

1. Man benutzt nicht die regelmäßige Kettenbruch- 
entwieklung — d.h. die mit den Zählern 1 —, weil 
deren Periode schwierig zu berechnen ist, obwohl sie 
für kleine Zahlen (und nur darauf bezieht sich ja unser 
Verfahren) tabuliert ist. Man benutzt vielmehr besser 
die unregelmäßige Entwicklung: 


b 
Be y 2 Se Ale 
Yr = a Ir b a Ar D) a Si b 
2a-+b .% 
20 +... 
In Wahrheit braucht man diese Entwicklung nicht 
explizit aufzuschreiben, sie dient nur zur Ableitung und 
Begründung der Methode.) 

2. Nunist es aber nicht notwendig — wie es sonst ist 
—, daß man die Näherungswerte sukzessive berechnet, 
sondern vermöge der quadratischen Irrationalität be- 
steht ein einfacher Zusammenhang zwischen dem i-ten 
Wert z;/n; und dem 2i-ten z3;/ng;, nämlich, wieman 
einfach zeigt: 


2i=2it un, Ni=2ziMmi. . . «+ (2) 


Dies ist ein spezieller Fall des allgemeineren: aus dem 
i-ten Wert z;/n; und dem k-ten Wert zx/nz den (i+k)- 
ten Wert z;;z/ni+z zu finden, nämlich: 


Rirk — Zizk + CNiNK, Num= ink + 2eni (3). 


Mittels (2) kann man also auf einfache Weise nach- 
einander die Werte der Ordnung 1, 2,4, 8, 16, 32,... 
finden. Gegenüber der gewöhnlichen Methode ist dies 
eine bedeutende Arbeitsersparnis. 

Aus (2) sieht man, daß, wenn z; zu n; prim ist, der 
folgende Wert zsi/nzinur durch einen Faktor von x 
eventuell gekürzt werden kann. Da dasselbe aber auch 


für 2;/ni gelten müßte, so ist auch z,,/n,: irreduzibel, 


wenn z;/n; es ist. Ähnliches gilt für (3). 
3. Für manche Zwecke ist es vorteilhaft, (2) umzu- 
ormen: Es ist 
Zu= 2, + en; = (@d;— and) + 2}; 
— (»2 ey?‘ ‚m2; 
= (2? —ın?)? + 2cnd. 
Aus dieser Beziehung 
a — any; = (a —aen?) 
folgt, wenn 
2 /n, = al, 2,[n, = a + b/2a = 
= 4-+ b,/a, mit b, prim zu a, 
gesetzt wird: 


[4 
2,51 (du)? + 2en, Eee  (A)R 


4. Ist man durch Kombination von (2) mit (3) oder 
mit (2) allein zu einem Näherungswert Z/N gekommen, 
der noch gerade in das Quotientenwerk der Maschine 
paßt, so wird Z/N auf der Maschine in einen Dezimal- 


KEN oIl daher i folgenden eine Methode entwickelt 
Es sol r im fo | 


Mittel zu bil- 


8 und besitzt 
A eh 
so viel ul 
5. Einen 


keit, ei isch ag u Mr , 
man nun gemäß Formel (l), wenn 
dividiert. So bekommt man zN/Z. DaZ so 
wurde, daß es indas Quotientenwerk paßt, so 
die Division auf der Maschine ohne weiteres a | 
Es un auf dieselbe Zahl von Ziffern gerechnet wie 
bei f 

6. Jetzt bildet man nach (1) das arithmetische Mittel 
M, welches jetzt auf t so viele Dezimalen genau 
ist wie Z/N. > 


7. Auch M kann aber noch verbessert werden. Da 


M = (Z|N + zN/Z)j2 N 


ist, so ist 


zem—la......) 


wo 

d=ZIN—ıN]Z= (Z2—ıN2/ZN=C/ZN. 
Hier ist C eine leicht nach (4) berechenbare ganze Zahl. 
Somit ist 
Ye=M-—a@l8Yz +... = M—-OJZN8 Ye +... 
oder, wenn man für Vx den Wert ZIN einsetzt: 


Ye=M— BZPN +... .. . (ba). 


Man bekommt also zu M nochmals so. viele genaue 
Ziffern hinzu wie Z/N hat. 

Bevor wir die Entwicklung noch weiter treiben, 
wollen wir ein Beispiel rechnen. 


Es liege zunächst eine Rechenmaschine mit 6stel- 


ligemQuotientenwerk vor. /6soll auf mög- 
lichst viele Dezimalen berechnet werden. 


Hier ist « = 2, also 
zıln, = all, z,/n, = 5/2. 
Wegen 23 — 6n3 = 1 wird nach (4): 
4=1+2-2=-9,,-%, 
2: = 1-+ 12: 20? = 4801, n, = 1960. 
Hier sind wir vorläufig am Ende, da z,, und n,, zu 
groß werden. Wir wenden aber, da andererseits 2, und 


n,; uns zu klein sind, Formel (3) an und zwar für? = 8 
und k = 4 und bekommen: 


21. = 470.449, n,, = 192 060. 


Obwohl z,, zu groß wird, kann man doch hier noch 


einen Schritt weiter gehen, da z,, durch 2 teilbar ist: 


Man bekommt nach (3) für?! = 12, k=1: 

23 = 2:1046629, n,,= 854569. 
Nun ist doch z,, noch gerade zu groß geworden, denn 
der zweite Faktor hat 7 Stellen und läßt sich auch nicht 
schnell zerlegen. Wir bleiben daher bei z,,/n,, und fin- 


den, da 1/n,, = 5:'10%, also 1/n?, =2-104, auf 
35 Dezimalen: 


ZIN = 2,44948 97427 88711 86087 68093 
30417 57784 — 21, N». 


Ähnlich hat man 6 N/Z = 1152360/470449, also 


6N/Z= 2,44948 97427 77644 33551 77713 20589 
79412. 


Das arithmetische Mittel beider gibt: 


M = 2,44948 97427 83178 09819 72903 25503 
68598 


| 
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und ist auf doppelt so viel Stellen etwa genau wie z/n, 
also auf etwa 22 Stellen. Die anzubringende Korrektur 
ist —1/8Z®N: 


— 1/8Z?N = — 62 50797 63722 - 10>28, 
so daß: f 


Ve = re 97427 83178 09819 72840 74706 


Bei einer Maschine mit Sstelligem Quoti- 
entenwerk bekommt man: 


216 —=1+ 12 = 22 = 46099201 
N = 225g = 18819920 


Zsl%ıs = 2,44948 97427 83178 67451 08374 
53081 62840 22461 30695 56087 

6 Nyolzıs = 2,44948 97427 83177 52188 37306 
96330 28997 61755 08985 50281 

Mittelwert = 2,44948 97427 83178 _ 09819 72840 
| 74705 95918 92108 19840 53184 
— 1/8 23,0, = 6779 72448 72359 85922 
V6= 2,44948 97427 83178 09819 72840 

74705 89139 19659 47480 67262. 


Genau dieselbe Methode — wenn auch ohne Korrek- 
tion —d?/8x, welche nochmals so viele Stellen liefert 
wie die erste Approximation — hat nun auch Gauß 
in späterem Alter angewandt. In Kapsel 44 seines 
Nachlasses findet sich das Beispiel der Wurzel aus 2. 
- Er gibt dafür den Näherungswert z/n, wo 


2 = 215807? + 2: 195025, n = 2 195025 - 275807. 


Dies ist in der Tat genau gemäß unserer Formel (2). 
Es müßte also 275807/195025 ein Näherungswert der 
Kettenbruchentwicklung von Wurzel 2 sein. In der 
Tat ist er der l14te solche Wert. Dann aber ist z/n 
der 28te. = 

Maennchen führt in seinem genannten Kom- 
mentar auch dieses Beispiel (auf p. 14/15) an. Doch 
trifft sein Kommentar nicht den Kern.der Sache. Dieser 
ist vielmehr in unseren obigen Entwicklungen zu sehen, 
die sicherlich Ga uß auch schon bekannt waren. Ver- 
mutlich war der Ga ußsche Gedankengang dabei fol- 
gender: 

Da 275807/195025 ein Näherungswert von Wurzel 2 
ist, so ist auch der Quotient von 2 und diesem Nähe- 
rungswert ein solcher, also 2 - 195025/275807. Das 
arithmetische Mittel beider gibt dann obigen Wert z/n 
und ist auf doppelt so viele Dezimalen genau. 

Ob Gauß nun wußte, daß dieses arithmetische 
Mittel der 283te Näherungswert ist, läßt sich nicht 
nachweisen. Wahrscheinlich war es so, denn der Zu- 
sammenhang ist ja ganz einfach. Wenn er aber dies 
wußte, so wußte er wahrscheinlich auch noch mehr. 
Denn es scheint mir so, als ob die Zähler aller jener 
Kettenbruch-Näherungswerte Primzahlen sind. Das 
hat zur Folge, daß man bei der Reihe der Näherungs- 
werte zu einer Grenze kommt, wo der Nenner sich nicht 
mehr in Faktoren von — sagen wir — 8 stelligen Zahlen 
zerlegen läßt, so daß die Rechenmaschine von da an 
nicht mehr gebraucht werden kann. 

In unserem Beispiel der Wurzel aus 6 wäre dies dann 
der nächste Zähler 


2, = 1 + 12n2, = 425027 26656 76801. 


Er wäre nach meiner Vermutung eine Primzahl. In 
der Tat habe ich große Anstrengungen gemacht, um 
ihn zu zerlegen, was mir nicht gelang. 

Damit sind wir zugleich zueiner Erweiterung 
unserer Methode gekommen, die gestattet, die 
Näherungswerte höherer Ordnung zu berechnen, auch 
wenn die dabei auftretenden Zähler und Nenner das 
Quotientenwerk der Maschine weit übersteigen. 

Um die Entwicklung jetzt weiter zu treiben, be- 
zeichnen wir den Näherungswert der Kettenbruch- 
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Entwicklung von /6 mit zn. Dann ist der Näherungs- 
wert doppelter Ordnung gleich } 


a9, = (1 + 12n2)/2zn = zn — dj2, 
d=zm—6nle=1jen. 


(Hiermit ist zugleich der Fehler des Näherungswertes 
z|n bestimmt zu fast genau 1/zn). Offenbar kann ich 
nun fortsetzen: 


2,|n, = 2m, — ıl2; 
und in unendlicher Reihe: 
6 = zu./n. = en - 42 —d,/2—d,j2—d,]2—... (6) 
Nun ist 


wo 


we d,=1kın 


d,/d = 1/22, = 1/24n2. 


eh z=2nV6, also d 2 1/(m? Y6) ist aber dies 
gleich 

d,/&=2d/]10, also d,=0,1d. 
Ähnlich d,=0,1d?. Das heißt die Reihe () 
konvergiert quadratisch. 


Ist nun 
z|n = 46099201/18819920 , 
soist d= 10, daher d, = 10 "dl!° —10 255, Bricht 
man also Reihe (6) hinter d, ab, so ist der erhaltene 
Ausdruck bis auf 255 Dezimalen genau. 

Für die praktische Berechnung handelt es sich jetzt 
noch darum, die Reihe so umzugestalten, daß man nur 
8 stellige Divisoren bekommt. Bezeichne ich nun mit 
A,, As, A,,.. . . stets das vorhergehende Glied der Reihe 
absolut genommen, so ist 
d, = 1/z,n, = 1/2zn(1l + 120?) 

= dj24n? — A,[12n°+ A,j12m2 — A,[12n + .... 

— 4,,/12n2. 


al 


Hier werden alle Glieder auf 255 (und noch etwas mehr) 
Dezimalen berechnet. Mit dem 15ten Gliede bricht 
man ab, da das 16-te Glied zu klein.wird. Ferner ist 

d, = d,[24n? — A,/12n? + A,/]12n?—.... 

= d, [96 2?n? — A,/482?n? + A,]482n? —... 

— A,[48 2202. 

Schließlich ist noch 

d, = d,/96 2?n? — A,[482!n? +.... 

= FH Her Wo 

Q = d,/96z?n} = d,]55296 2?n® — A,/6n? + A,/8n? 
— 4,/9n? + 4A,[9,6n? — 4,/10n? 
+ 74,/72n? — 4,[10,5n? j 

Q, = Q]48 zini = Q/27648 2? n® — A,/6n? + A,|8n? 
— 4,/9n? + 4,/9, 6n? 

Q, = Q,[27648 220° = 15 - 1029, 
Damit ist die Rechnung beendet. 
Den Haag. E.Bodewig. 
Allgemeine Differenzenausdrücke für die Ablei- 

tungen einer Funktion y (x). 

Für die m-te Ableitung y'”° (x,) einer Funktion y (®) 
an der Stelle x, wird im folgenden zu beliebig ver- 
teilten Interpolationsstellen x, ein Differenzenausdruck 
(‚‚finiter Ausdruck‘“) samt Fehlerglied in geschlossener 
Form angegeben. Finite Ausdrücke finden bekannt- 
lich Verwendung beispielsweise bei der numerischen 
Differentiation und bei Auflösung linearer Differential- 
gleichungen. 

Das Ergebnis werde vorweggenommen. 

a) Die Funktion y (x) sei in einem z-Intervalle, das 
die sämtlich voneinander verschieden vorausgesetzten 
aber sonst beliebig gelegenen (n+ 1) Stellen x, 


20 +&oh(e = 1,2,...n) umfaßt, mindestens (n + 1)- 
mal stetig nach x differenzierbar; es bedeute 77 
d’y(®) 
yr) (®) I d 2 ; 
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4 


& 35/24 15 "Pi 
x ie _ 381120 — 7/36 Es 
— 7/6 | 21/20 
— 49/36 


p=o = —_ —_ 1/6 5/12 17/24 49/48 | 967/720 
En — —_ — 1/2 — 5/6 — 25/24 | — 7/6 | — 889/720 
2 = 1/2 0 —5/12 | —35/48 | — 77/80 
31 — 16 5/6 , | 5/12 0 — 49/144 
4+| — 52 25/24 | 35/48 49/144 
5 | — 17/24 6%: 1.20 71j80 2=3 
6 | — 49/48 | 889/720 
7 | — 967/720 ‚ 
1 
p=0 — — — — 1/24 7/18 
1 = _ = — 1/6 — 3/8 2 
2 —_ —_ 1/4 1/4 7/48 
3 2 — 1/6 1/4 7/18 7/18 
4 1/24 — 3/8 7/48 7/18 
5 1/8 — 7/12 0 
6 35/144 — 7/9 
7 7/18. 
p=0 = m — 1/120 7/240 | 161/2520 
1 — — — 1/24 | — 7/60 | — 77/360 
2. — 1/12 7/48 21/120 
3 — 1/12 0 7/72 
4 1/24 — 7/48 — 7/72 
5 |— 1/120 7/60 — 21/120 sie 
6 | — 7/240 | 77/360 
7 | — 1161/2520 


= 
1 


IOoUPODH-O 


=} 
l 
1.1 IE 


SS[OUTPRWOD-O 


— 1/5040 
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Ts Ar a a A a m Fe Eu 


In (1), (2) ist 
Fim,n,0)=1— 09 +02 —...+ 
+ Elan: a Sm2 + (1! es Sm-1 ; 


(2<m<n) 
Fi, Nn,%)=1. 


Dabei bedeutet S$, die Summe der Kombinationen 


N 1 il 
sämtlicher Elemente EUR A-ten Klasse ohne 
1 N 


N 
Wiederholung. I (%&—%,) ist (bei festem oe mit 
v 


1<eo<sn) zu bilden für v=1,2,...n mit Aus- 
nahme von » = 0. Anstatt der unterklammerten Aus- 
drücke n + 1,» + 1— m, n— m können dann (aber 
nur dann) die um 1 vergrößerten Zahlenn +2,n +2 
— m, 9 + 1— m genommen werden, wenn m gerade 
ist und überdies die » Punkte ©9+ &h symmetrisch 
zu &, liegen. 

b) Bira+tl=p+g+1 äquidistanten 
Interpolationsstellen 2, +oh mit e=0,1,2,...q, 
—1,—2,...—P(mq4>0 ganzzahlig, p+g >|) 
folgt speziell aus (1), (2) fürl<m<p+g : 


(-1)m m! 2, 
hm 


D(m,p,q, 0) 
2 Bla oil reh + Rm,»,a... (3) 


miptgt Sr (= De, 


yr)(z,) = "y(2) - hm om 


ep 


mit 
q 
m!p!g! PHraril—m De 
= — m, — 10» 
Bm, p,a FRRPSITK | 1) 
Ir 
ptatl=—m  B(m,9,94,0) : (P+ı+D 
un, 1 DE Del x 
® Pro (at? Ranch) 


(01-1, &l<)... 2; 


Dabei sei vorausgesetzt, daß y (x) in einem Intervalle, 
das sämtliche Interpolationsstellen enthält, mindestens 
2-+g-+ 1) mal stetig nach x differenzierbar sei. 
Er FAREEN 


In (3), (4) bedeuten 
(m, p,g; 0) == i—B2, = 0° 2, au ee 
+ (— 1)m-2 om-2 Zn E (— 1)m-1 om-1 a 3 
2<m<sp+tg 
®(1,2,9,0=]1. i 


Weiter sei &, die Summe der Kombinationen zur A-ten 
Klasse ohne Wiederholung der (p + q) Elemente 


1 1 ö 1 EN: d 
er m, l,—7,... bei 
Bas..7 > - PR el 
der q bzw. p Elemente 
1 1 1 1 
a am = ginn 
beiip=0bzw.gq=0 
q 
In 2” ist zu summieren über o=—p, — p+J],... 
e=-P 


—1,1,2,...g mit Ausnahme von ge = 0. Diein (4) 
und anschließend auftretenden unterklammerten 
Zahlen p+g+1,p+g+ 1-— m können um 1 er- 
höht werden bei p =g und geradem m. 


Der gesuchte finite Ausdruck für y®) (x,) wird dar- 
gestellt durch die rechten Seiten von (1) bzw. (3) bei 
Vernachlässigung der Fehlerglieder Rm,n bzw. Kim, p, q. 
Er ergibt sich, wenn man durch die (n + 1) Punkte 
[2 Y(ao), [üo+ &e h, Y(a, + ae h)] das (eindeutig 
bestimmte) Interpolationspolynom höchstens n-ten 
Grades in x legt und dessen m-te Ableitung nach z an 
der Stelle x, als Näherung für y(m) (x,) nimmt. Aus- 
gehend von der Lagrangeschen Interpolations- 


formel läßt sich die Beziehung (1) in der gleichen Weise 
herleiten wie in dem (für allgemeines m) bereits ander- 
weitig 1) betrachteten speziellen Falle symmetrischer, 
äquidistanter Interpolationsstellen. Die Darstellung 
(2) des Restgliedes kann man, wie man sich leicht über- 
legt, sofort aus dem finiten Ausdruck (1) erhalten, wenn 
KR Ko letzterem y(2,) durch 0 und y(x, + %e h) 
ure 


lan! 
(nr. 1j1 


bzw. bei geräadem m und hinsichtlich «, symmetrisch 
ee Interpolationsstellen 2, + &oh y(z, + &oh) 
urch 


ya do'&oh), 


(Keh)" * 
(n + 2)! 
(bei natürlich jeweils voneinander verschiedenen Zah- 
len |9e| < 1) ersetzt. Auf weitere Ausführungen zur 
Herleitung von (1), (2) darf hier verzichtet werden. 
Bei symmetrisch zu z, verteilten Interpolations- 
punkten lassen sich ersichtlich die Differenzenaus- 
drücke in (1) und (3) noch etwas vereinfachen. Für den 
Fall äquidistanter symmetrischer Stellen 2, +oh 
sind diese vereinfachten Darstellungen in den in der 
Fußnote 1 angeführten Veröffentlichungen angegeben. 
Zur Erleichterung der numerischen Berechnung 
spezieller finiter Ausdrücke?) nach (3) werden nach- 
stehend einige Werte 2, zahlenmäßig angeführt. 


Stuttgart. E. Pflanz. 


2) Allgemeine Formeln für y(m) (z,) beizu x, symmetrischen 
äquidistanten Interpolationspunkten sind für m =1}; 2 her- 
geleitet bei L. Collatz, Schriften des Math. Sem. u. d., 
Inst.f. angew. Math. d. Univ. Berlin, Bd. 3 (1935) S. 8S—9 und 
(auf andere Weise)fürallgemeines mbeiE.Pflanz,Z.angew. 
Math. Mech. 17 (1937) 8.296. { 

2) Zusammenstellungen von üblicherweise benutzten finiten 
Ausdrückenfindensich beiL.Collatz, Eigenwertprobleme 
und ihre numerische Behandlung. Leipzig 1945, 8. 261if. — 
E. Kamke, Differentialgleichungen: Lösungsmethoden und 
Lösungen, Bd.I. Gewöhnliche Differerentialgleichungen, 
Leipzig 1943, S. 2251f. — G. Schulz, Formelsammlung zur 
praktischen Mathematik. Sammlung. Göschen, Bd. 1110. 
Berlin und Leipzig 1937, 8. 126ff. 
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Zwei Sätze über die Nullstellen der Bessel- 
Funktionen. 

Wie im folgenden gezeigt wird, lassen sich aus der 
Darstellung des Temperaturverlaufes beim Abschrecken 
einer n-dimensionalen Kugel zwei Sätze über die Null- 
stellen der Bessel-Funktionen gewinnen, die bisher 
anscheinend noch nicht bekannt waren. ‘ 

Da der Temperaturverlauf beim Abschrecken Kugel- 
symmetrie aufweist, wird der Temperaturverlauf beim 
Abschrecken der n-dimensionalen Kugel beherrscht 
durch die Differentialgleichung: 

oT * x» [®T n—107 
a dmtTT er A) 


Hierin ist « die Wärmeleitfähigkeit, c die spezifisch- 
Wärme und o die Dichte der Kugel. h 
Nehmen wir an, daß die ,‚Oberfläche‘‘ O der n-dimen- 
sionalen Kugel mit dem Außenradius r, und dem ‚,‚In- 
halt‘‘ V und der gleichmäßigen Anfangstemperatur 7’, 
im Zeitpunktt= 0 auf die Temperatur T = 0 gebracht 
und dauernd auf dieser Temperatur gehalten wird, so 
wird die Gl.(1) unter Berücksichtigung dieser Anfangs- 
und Randbedingungen befriedigt durch den Ansatz; 


> ir ee "el sol) 2) 


To 


Diese Summe ist zu erstrecken über sämtliche 
Wurzeln r = 7% der Gleichung: 


ee Jp()=0. 


e (2) Ko : ARE 
a Tr 
EP Inland a 
Ar = To — ae 
JE ETae 
Wir erhalten: SR 2-7, 
| ax Inr1(08) 


und damit ergibt sich für den Verlauf der Temperatur 7 


während des Abschreckvorganges die Gleichung: 


u En g 
ea mr? () (3) 
r 
Ina) 
- In+1 (ar) 
Kir wol, pn 5 geht diese Gleichung in die be- 


kannte Fourier-Reihe für den Temperaturverlauf 
beim Abschrecken einer ebenen Platte mit der Dicke 
2. r, über. 

Der durch Gleichung (3) gegebene Temperatur- 
verlauf muß folgenden Bedingungen genügen: 

a) Die gesamte Wärmemenge, die aus der „„Ober- 
tläche‘‘Oder, ‚Kugel‘ herausströmt, muß gleich demge- 
samten Wärmeinhalt der ,„Kugel‘‘ mit dem, , Volumen“ 
V zu Beginn des Wärmeausgleiches, also für i=0, sein. 

b) Die Temperatur im Mittelpunkt der Kugel muß 
zu Beginn des Ausgleichsvorganges, also für t=0, 
gleich 7, sein. 

Aus diesen Forderungen lassen sich nun Aussagen 
über die Nullstellen der Bessel-Funktionen gewinnen. 

Die Forderung a) wird befriedigt, wenn die Be- 

ziehung: - 


oT ) 
.0/:T,=—x:0- — di 
IR s r & | | or 1=r, 

0 

erfüllt ist. Nun gilt für die n-dimensionale ‚Kugel‘ 

die Beziehung: war, | 
OR 

Mit oo aM. 

£ (= ) e'c ; 
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kl 
gelangen wir dann zu folgender Aussage über die Null- 
stellen der Bessel-Funktionen: 


oo 


1 TE PN 
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n=ı 


Jp(&r) 0, 


oder: 

Die Summe der reziproken Quadrate der 
Wurzeln der Bessel-Funktion mit dem Para- 
l 1 


2n Mp+l)' 
Für die verschiedenen ganzzahligen Werte von n 
führt diese Gleichung (4) zu folgenden Reihen: 
aan=1, -— 2: ebene Platte mit der 


Dicke 2 r;: 


eure 


ai (art. 3 3 5 


meter p =; —1 ist gleich 


2 


(überführbar in die bekannte Reihe: 
ala |) 


1 1 
a a 


b)n=2,p=0: Kreiszylinder mit dem Halb- 
messer rg: 


oe; p= + 4: Kugel mit dem 


y er 


1} . i 6 bs 

unse NE 7 2 
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Er 

“ er 
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DEN TE 
Zum mtamt ie; 
(wel.da)) J) 


| d)n=4,p=1: vierdimensionale Kugel mit 
Ir hi 


dem Radius rg: 


Ne aa) 
Du, für Ja) =0,.>0 7 


usw. *) h 
Die Forderung b), daß die Anfangstemperatur im 
Zentrum der Inge] gleich T, sein soll, ergibt aus 


Gleichung (3), indem wir t = 0 und r = 0 setzen, die 
Beziehung: \ 32 
1 1 Sa 1 a} 


für er -Ip(a,)=0: 


Die hierdurch gegebene alternierende Reihe weist 
lediglich für Werte n< 3 Glieder mit abnehmendem 
Absolutbetrag auf, so daß wir uns darauf beschränken, 
ihren Wert für nr = 1 und für » = 2 aufzuführen. 


re ebene Platte mit der 


a n=1,9P =— 
Dicke 2 r;: 


(bekannt als Leibnizsche Reihe in der Form: 
nt 1 1 1 . ) 
he 


b)n=2,p = 0:Kreiszylinder mit dem Halb- 
messer rg: 


VIE 1 

=: ee I 
Diese Identität scheint nicht bekannt zu sein; Jahnke- 
Emde (III. Aufl., S. 166) bringen lediglich die ähn- 
liche Formel: 


wo ı wo 1 
3 ee E20 2 Eee © 
Er PATER: 
für J,(Pr) = 0. 


für J,(ar) =0 (öb). 


München 23. W. Harries. 


*) Herr Professor Magnus macht mich darauf aufmerksam, 
daß man die Gleichung (4) auch als Kongruenz aus der Formel für 
die Produktzerlegung der Bessel-Funktionen gewinnen kann. 


Es ist: 1,9? © 2 
pie) = AR -II (' Serß 
k=1 Dr 
wobei 2=& -1,—2,—3,..., aber sonst beliebig ist, und x 
die %-te von Null verschiedene Nullstelle von J„(z) ist. Aus dieser 


Entwicklung folgt Gleichung (4) als Kongruenz und ähnliche For- 
meln für die Summe der reziproken 4-ten Potenzen usw. der Null- 
stellen durch Vergleich der Koeffizienten der Potenzen von z auf 
beiden Seiten nach derselben Schlußweise, mit der man aus 


©o 
snzo=ng: 17, (1-2) 


k* 
die Formelh 


ableitet. 


Diese Ableitung der Gleichung (4) zeigt, daß diese Gleichung 
für PENERIBR Klaade bzw. ngilt, mitalleiniger Ausnahme der Werte 
Re En a 

Fbenso hatten auch wir in unserer Ableitung in keiner Weise 
Gebrauch von der Ganzzahligkeit von n gemacht, wenn es auch 
schwer fällt, sich Wärmeausgleichsvorgänge in einer Kugel vor- 
zustellen, deren Dimension keine ganze. Zahl ist, 
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y "Es ist mehr eine Sammlung von Rohmaterialien 
zum Buchtitel als ein Buch über den durch den Titel 
bezeichneten Gegenstand. Um den Leser vor einer 


 * Enttäuschung über den bekannten. Verfasser zu be- 


wahren, möchte bemerkt werden, daß die im ersten 


Kt + Kapitel dargestellte Logik die veraltete formale 
Logik der Abstraktions- u. Inhalts-Umfangstheorie 


ist, daß als „Grundlegung der Analysis‘ nıchts ge- 


‘ braeht wird, was nicht in jedem Lehrbuch der In- 
finitesimalrechnung stände, und daß das Schluß- 


 kapitel ‚Mathematik und Erkenntnislehre“ für den 


genau da aufhört, wo die Probleme anfangen. 

Für wen ist das Buch geschrieben? Laut Vorwort 
' für den Lehrer und Studierenden, ja auch für den 
interessierten Schüler! Dazu nur eine Bemerkung, 


weil das Vorwort „moderne Stoffpläne‘ erwähnt, dıe 


„Rückblicke auf die Mathematik unter philosophi- 
schen Gesichtspunkten‘ vorschreiben. Dem Primaner 
möchte diese Kost anders und in anderer Zubereitung 
geboten werden. Eine Vorlesung vor künftigen Leh- 
rern dürfte sich erst: recht: nicht mit: Rohmaterialien 
zufrieden geben, sondern hat diesen Stoff in klarer 
philosophischer Durchdringung und mit eingehender 
Würdigung der verschiedenen Mathematikauffas- 
sungen zu bringen. 
> Dresden. 


1 


M. Draeger. | 


Th. Pösehl (o. Prof. a. d. Techn. Hochschule Karls- 
ruhe), Einführung in die analytische 
Mechanik. (Reihe: Mathematische Methoden in 
Naturwissenschaft und Technik.) 166 S. m. 37 Abb. 
Karlsruhe 1949. G. Braun-Verlag. Preis: Halbl. geb. 
10,— DM. 

Infolge der großen Fortschritte, welche im Laufe 
der letzten Zeit auf dem Gebiete der theoretischen 
und angewandten Mechanik, der Atommechanik und 
der Himmelsmechanik erzielt wurden, knüpft nicht 
nur die Fachwelt, sondern auch die breitere Öffentlich- 
keit an jede Neuerscheinung innerhalb dieses Wissens- 
zweiges besondere Erwartungen. Dies gilt um so mehr, 
als erst: vor wenigen Jahren das Buch Mechanik 
eines der größten deutschen theoretischen Physiker, 
‚Prof. Dr. Arnold Sommerfeld (als erster 

Band seiner Vorlesungen über theoretische Physik) 
erschienen war, und es liegt nahe, die vorliegende 
Einführung in die analytische Me- 
ehanik des durch eine große Zahl wertvoller Bei- 
träge auf dem Gebiete der Mechanik bekannten 
Prof. Dr. Th. Pöschl mit dem thematisch so 
stark verwandten Buche von Sommerfeld zu 
vergleichen. Hierbei zeigt sich zunächst —, wie nicht 
anders zu erwarten — weitgehende Übereinstimmung 
in der Wiedergabe der grundlegenden, auf La- 
grange, Hamilton, Jakobi u.a. zurück- 
gehenden dynamischen Theorien materieller Systeme, 


- "die man heute als das klassische Fundament der 


theoretischen Mechanik ansieht. Hinsichtlich der 
weiteren Behandlung des sehr vielseitigen Stoffes 
weisen beide Bücher jedoch wesentliche Unterschiede 
auf. Im Pöschlschen Buche, das laut Vorwort 
eine Wiedergabe der in den letzten Semestern an der 
Technischen Hochschule- Karlsruhe gehaltenen. Vor- 
lesungen des Verfassers darstellt und entsprechend dem 
gewählten Titel mehr der mathematischen Durch- 
dringung des Stoffgebietes gewidmet: ist, überwiegen 
die der Theorie besonders zugänglichen Abschnitte der 
Mechanik, was bei dem an sich geringen Umfang des 
Buches naturgemäß etwas auf Kosten der Reichhaltig- 


v gege 


: sche Integralgleichung, Anwendungsb 


Uneingeweihten die Dinge eher verwirrt als klärt und 


geben) und führt den Leser geschie 
ıathematisch-formalistischen Gipfelpı 


stellung und Lösung der Hamil on 
schen Differentialgleichung, einfache Fälle d 
tentheorie, verallgemeinerte Koordinaten, jedo 
Bezugnahme auf die Relativitätstheorie). Demge 
über wird in dr Sommerfeldschen Mecha 
der theoretische Untergrund aus pädagogis 
Rücksichten mehr in wohlabgewogener Knapp 
dennoch vielseitig dargeboten und durch mannigf 
eingestreute Hinweise auf besondere Ergebniss 
Experimentalphysik belebt, so daß der Leser etwas 
unmittelbarer angeregt wird. Dessen ungeachtet fin- 
det auch der’ Leser der Sommerfeldschen 
Mechanik, der sich noch nicht in die weiteren Werke 
der theoretischen Physik vertiefen möchte, bei 
Pöschl wertvolle theoretische Ergänzungeninkurzer, 
sprachlich eleganter Form, die er sich sonst mühsam 
aus der Literatur heraussuchen müßte (z. B. dieHin- B 
weise zur Atom- und Quantenmechanik im letzten 
Teil des dritten Abschnitte). Das Pöschlsche 
Buch, das sich auch durch sorgfältig ausgeführte Ab- 
bildungen auszeichnet (im Gegensatz zu Lagrange, 
der in seiner m&canique analytique Wert darauf legte, 
ohne Abbildungen auszukommen), ist daher als neue 
Ergänzung der Fachliteratur sehr zu begrüßen, umso 
mehr, als in der gegenwärtigen Situation die Stan- 
dardwerke der rein-analytischen Mechanik (z.B. 
Whittaker) vielfach schwer zugänglich sind. 


Dresden. ‚H. Neuber 


Dr.-Ing. W. de Beauclair, Untersuchungen 
überdieFouriersynthesederLadungs- 
verteilungin Kristallen. Band I. Ver- 
fahrenund Gerätezurmehrdimensio- 
nalen Fouriersynthese. VII-+71 S. mit 
38 Abb. -Berlin 1949. Akademie-Verlag. Preis 
15,— DM. 

Das Buch gibt eine Schilderung des augenblick- 
liehen Standes der Arbeiten zur mehr-, insbesondere 
zweidimensionalen Fouriersynthese, wie sie vor allem 
in der Kristallforschung gebraucht wird. Da diese 
Verfahren vielfach auf die eindimensionalen zurück- 
greifen, werden zunächst die verschiedenen Methoden . 
hierfür dargestellt. Sodann werden rechnerische und Er 
instrumentelle Verfahren zur zweidimensionalen Syn- 
these geschildert und, soweit sie dem Verfasser brauch- 
bar zu sein scheinen, so ausgearbeitet, daß sie auchvon Fe 
Hilfskräften angewandt werden können. Insbesondere 
wird ein Verfahren, das Phasentafeln benutzt, be- 
schrieben. Derartige Tafeln sind im zweiten Band 
gegeben; sie dienen der Bestimmung von 323 A7x cos 
2rhx- cos Arnky und geben die Werte von cos 2rrhx 
-c0o82rky für ganze Vielfache von 1/24 der Orts- 
koordinaten. Auf diese Blätter werden Transparent- 
blätter mit den Koeffizienten Anr gelegt, so daß die 
zu multiplizierenden Größen neben- oder übereinander 
stehen. Die Multiplikationen und fortlaufenden 
Additionen werden mit der Rechenmaschine aus- 
geführt. Auch eine Reihe mechanischer, elektrischer 


‚und optischer Verfahren wird angegeben. Neben dem 


Meßrollengerät erscheint dem Verfasser die Über- 
lagerung mittels elektrischer Frequenzmodulation die 
beste Methode. Zum Schluß wird über die bisher gelei- 
steten geringen Vorarbeiten zur dreidimensionalen Syn- 
these berichtet, die wenigstens zur Zeit noch der zwei. 
dimensionalen wesentlich an Wichtigkeit nachsteht, 


Sein Ziel, ein brauchbares numerisches Verfahren 
für die zweidimensionale Synthese zu schaffen und 


% 


darüber hinaus eine’ Basis zu geben, anf der weiter. 
7 aus ei s zu geben, 1te 


gebaut werden kann, um zu anderen bra 
numerischen oder instrumentellen Verfahren zu kom- 
men, hat der Verfasser mit diesem Buche voll erreicht. 


Dresden, Willers. 


Dr. Fritz Oberhettinger (Dozent. f. Math. a, d. Univ. 
Mainz) und Dr, Wilhelm Magnus (Prof. d. Math. a. d. 
Univ. Göttingen, Anwendung der ellip- 
tischen kit enen in Physik und 
Technik. (Die Grundlagen der ‚mathematischen 
Wissenschaften in Einzeldarstellungen mit besonderer 
Berücksichtigung der Anwendungsgebiete, heraus- 
‚gegeben von W. Blaschke, R. Grammel, E. Hopf, 
F. K. Schmidt und B. L. van der Waerden.) VII -+ 
126 S. mit: 54 Abb. Berlin, Göttingen, Heidelberg 1949, 
Springer-Verlag. Preis geh. 15,60 DM, geb. 18,30. DM, 
Die Verf. haben in einem kürzlich in zweiter Auflage 
erschienenen Buch Formeln und Sätze für die speziellen 
Funktionen der mathematischen Physik zusammen- 
gestellt. In dem voliegenden Buche zeigen sie nun, 
wie die im siebenten Kapitel des ersten Buches be- 
handelten elliptischen Funktionen auf Probleme der 
Physik und Technik angewandt werden, Dazu haben 
sie eine Reihe von Beispielen aus den verschiedensten 
Gebieten ausgesucht. Behandelt werden Fragen aus 
dem Gebiet der konformen Abbildungen und der 
Greenschen Funktion, aus der Elektrostatik, der 
Hydro- und Aerodynamik und aus anderen Gebieten 
(Mathematisches Pendel, Potential eines Ellipsoides, 
gedrückter Stab, Tschebyscheffsche Approximation). 
Da dem Physiker und Techniker die hierfür erforder- 
lichen Formeln und Beziehungen im allgemeinen. nicht 
gegenwärtig sein werden, sind sie im ersten Kapitel 
ähnlich wie in Kapitel sieben des erwähnten Buches 
zusammengestellt. Am Schluß sind dann ‚einige 
Zahlentafeln gegeben, die die numerische Auswertung 
der Ergebnisse gestatten. Die Beispiele des Buches 
sind sehr geschickt gewählt. Sie zeigen, daß die 
elliptischen Funktionen in den mannigfachsten Ge- 
bieten von Physik und Technik eine Rolle spielen. 
Die Darstellung gibt eine gute Anleitung für die Be- 
handlung solcher Probleme. Möge das Buch in die 
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Dr.-Ing. habil, Walter Kaufmann, 
werke. .Tei 
1.Bd., dritte ergänzte und 
VIIL + 314 S. mit 364 Abb. Berlin, Göttingen, - 
berg 1949. Springer-Verlag. Preis brosch. 25,50 DM. 
Das. Buch ist mit den beiden ersten Auflagen für 
zahlreiche Studierende technischer 
Männer der Praxis zum wertvollen Führer beim 
Studium der Baustatik geworden. Dies wird die neue 
Auflage ebenfalls für sich in Anspruch nehmen können. 
Sie behandelt wiederum ausschließlich die Statik der 
Stab- und Fachwerke, ist aber durch neue Beiträge 
ergänzt worden, die sich für den konstruktiven In- 
genieurbau in der jühgsten Zeit als zweckmäßig und 
wertvoll erwiesen haben. Im übrigen ist der bewährte 
Aufbau des Buches unverändert geblieben. Er be- 
handelt den Spannungs- und Verschiebungszustand 
im allgemeinen und enthält ausführliche Rechen- 
vorschriften für alle diejenigen Tragwerke, die besondere 
Bedeutung für den konstruktiven i ge- 
wonnen haben. Das sind die durchlaufenden Träger, 
die Bogenträger und Rahmen. Ihre statische Be- 
rechnung wird klar und zweckbestimmt Ba ie 
so daß jeder Leser das Buch voll Dank an den Ver- 
fasser aus der Hand legen und immer wieder gern zu 
Rate ziehen wird. 
Dresden. 


® 


K. Beyer. 


Götz Quarg, Wider den technischen 
Kulturpessimismus. 55 S. Düsseldorf 
1949. Deutscher Ingenieur-Verlag G.m.b.H. Preis 
brosch. 1,70 DM. 

Die kleine Schrift wendet sich gegen eine bestimmte 
Art unberechtigter Angriffe auf Naturwissenschaft 
und Technik wie sie in den letzten Jahren vielfach, 
besonders in Deutschland, erhoben sind. 

Dresden. Willers. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen 
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Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten) 


Dr.-Ing. Bruno Eck, Technische Strö- 
mungslehre, Dritte verbesserte und erweiterte 
Auflage. X + 398 S. mit 372 Abb. Berlin-Göttingen- 
see 1949. Springer-Verlag. Preis brosch. 
24,— DM. 


Prof. Dr. Werner Schaub, Weltraumflu 
(astronomische und physikalische Grundlagen). 93 8. 
mit 20 Abb. Bonn 1949. Ferd. Dümmlers Verlag. 
Preis brosch. 6,35 DM. 


Prof. Wilhelm Staufer und Oberreichsbahnrat 
Werner Splett, Durchhang und Zug- 
spannungen von Starkstrom-Frei- 
leitungen. IV-+10 S.mit 32 Zahlentafeln und 9 
Kurventafeln. München 1949. Richard-Pflaum-Ver- 
lag. Preis kart. 14..— DM. 


Dr.-Ing, habil Otto Luetkens, Die Methoden 
der Rahmenstatik, Aufbau, Zu- 
sammenfassung und Kritik VI-+ 
281 S. mit 38 Abb. und 9 Zahlentafeln. Berlin- 
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Berlin-Göttingen-Heidelberg 1949. Springer-Verlag. 
Preis brosch. 36,— DM, geb. 39,— DM. 


Dr. R. Rothe (ehem. Prof. a. d. Techn. Hochschule 
Berlin, Höhere Mathematikfür Mathe. 
matiker, Physiker, Ingenieure, Teil 
III: Flächen im Raume. Linienintegrale und mehr- 
fache Integrale. Gewöhnliche und partielle Diffe- 
rentialgleichungen nebst Anwendungen. Vierte Auf- 
lage. 236 S. mit 167 Abb. Leipzig 1949. B. G. Teub- 
ner-Verlag. Preis karton. 7,— DM. — TeilIV. Heft 
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